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Introducere 


Disciplina „Matematici aplicate în economie” este una dintre disciplinele generale 
pe care studenţii le studiază pentru ca pregătirea lor să fie completă. Ea se alătură altor discipline 
cu același caracter (limbile străine, dreptul, economia etc.) fără de care un specialist în economie 
nu poate să își exercite profesia. Trebuie precizat că în alte țări din Europa sau din afara 
continentului nostru, volumul de cunoștințe matematice predat studenţilor este mult mai mare şi 
ocupă programe de studii pe mai multe semestre. Ceea ce a fost introdus în această programă și 
s-a concretizat în prezentul curs reprezintă un minim şi un strict necesar. 

Disciplina este pe de altă parte una obligatorie. Ea este una dintre acele discipline 
cu credite şi fără acestea promovarea nu este posibilă. Numărul de credite în cazul disciplinei de 
faţă este egal cu cel al celorlalte discipline obligatorii (adică 5 credite). Este o nouă și importantă 
dovadă a locului care i se acordă în planul de învățământ al facultăţii noastre. 

Ar fi ideal dacă toţi studenții noştrii ar avea cunoştinţe de matematici elementare 
echivalente cu cele ale unui absolvent de liceu cu profil economic. Cum acest lucru nu este 
posibil, mulţi studenți provenind din licee cu alte profile unde matematica are o pondere mai 
mică, ne-am străduit să facem cât mai puţine referiri la cunoştinţele anterioare de matematică ale 
studenţilor. În consecință am cuprins în manual aproape tot ceea ce este necesar pentru 
parcurgerea şi înţelegerea textului. Elementele care sunt preluate pot fi aprofundate şi clarificate 
în activitățile tutoriale. 

Obiectivele cursului sunt multiple. Pe de o parte ne propunem să familiarizăm 
studenţii cu modalităţile uzuale de modelare matematică a unor fenomene economice. În acest 
scop vom obişnui studenții să extragă elementele caracteristice esențiale ale unui fenomen 
economic şi să le transpună în relații matematice. Acestei prime etape de modelate îi va urma 
etapa studiului modelului matematic rezultat şi, în final, interpretarea economică a rezultatelor. 

Pe de altă parte vom face cunoscute studenților două mijloace puternice de 
optimizare folosite în economie (metode algebrice şi metode bazate pe derivată) care au permis 
rezolvarea unor probleme clasice ale economiei (problema gestionării resurselor, problema 


transportului, problema repartizarii forţei de muncă etc.). 
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Nu în ultimul rând, studenților li se formează obişnuinţa de a gândi şi a acționa 
respectând regulile stricte ale matematicii. Nu ne propunem să împingem la extrem cerințele 
privind exactitatea calculului, dar nu trebuie să uităm că avem de-a face cu viitori economişti 
care vor utiliza mijloace materiale şi vor purta răspunderea pentru eficienţa şi corectitudinea 
muncii lor. 

Coroborând cunoștințele de matematică cu cele de informatică, absolvenţii noştri 
vor putea crea aplicaţii pentru computer care vor completa soft-ul deja existent în domeniul 
economic. Ei vor deveni nu doar utilizatori ci şi creatori de soft economic, ceea ce reprezintă o 
treaptă calitativă net superioară. 

Materialul prezentat este structurat în nouă unități de studiu. Această împărțire, 
împreuna cu planul calendaristic al activităţilor tutoriale, permite gruparea în trei părți, fiecăreia 
corespunzându-i o lucrare de evaluare. 

Conţinutul acestor lucrări va privi aspectul teoretic tratat în curs, dar şi întrebări 
de control sub formă de grilă la care studenţii vor trebui să răspundă. Întrebările de control la 
autoevaluare se află în materialul existent pe Blackboard şi este accesibil eşalonat studenților, 
etapă cu etapă. 

Pentru verificarea cunoștințelor şi evaluarea acestora s-au prevăzut trei teme de 
control. 

Prima temă de control privește elementele de algebră pe care studenţii le parcurg 
în primele trei unităţi de studiu. Ea se intitulează astfel „Elemente de algebră liniară”. 

A doua temă de control priveşte metodele algebrice de obținere a soluţiilor optime 
ale problemelor de economie. Această temă se intitulează “Programare liniară şi algoritmi 
asociați” şi utilizează material parcursă în unităţile de studio numerele 4, 5 şi 6. 

Ultima temă de control priveşte material parcursă în unitățile de studiu numerele 
7, 8 şi 9 şi este intitulată "Elemente de analiză matematică aplicate în economie”. Studenţii vor 
primi temele de control în cadrul activităţilor tutoriale şi le vor depune la termenele indicate la 
secretarialul ID. 

Pentru a asigura o asimilare uşoară a materialului teoretic, cursul şi activitățile 
tutoriale cuprind nenumărate exemple. 

Nu este de dorit şi nu încurajăm memorarea fără aprofundare a unor rezultate 


teoretice. 


Nu sunt necesare mijloace auxiliare (calculatoare de buzunar, tabele de valori, 
indrumare de calcul etc.), ci numai materialul editat. 

Evaluarea se face pe baza temelor de control (60% şi a evaluării finale 40%) cu 
condiția ca ambele să aibă măcar nota 5 (cinci). Dacă notăm cu M media celor trei lucrari de 


control 
Na + N2 + Na 
ai 3 


şi cu EvF nota la evaluarea finală (M > 5, EvF > 5), atunci nota finală este 


N; = 0,6 M + 0,4 EvF. 


Cuprins 


T AERO EEE sc zapada A aluat ta led drd a tata 8 a al aaa i a a 2 
2 CUPIS seu da iaca a ar d a a a a dia a a i ia 5 
3. USI Trâånsformări elementare; sui actul nd apa aa apoape 8 dag dat tao ba 6 
4; > 4052 Spatii haté ss aa e aaa ea a d aaa aa 0 ai 20 
Si USI Forme Pale ati GE acetic ata aaa a a aaa atata Bat 35 
6. US4 Programare liniară. Considerații teoretice... nenea nenea nn aan na 40 
Ti -USS Metoda SIMPLEX seci pct pe cae cata bad aa at da d aaa n lol a a d aa 50 
8, -US6-Přobleme de: transport setei ata a a ao d ac EA la 65 
9: -157:Seru:de numere reales sie eaină ca eta ad enter aa eo oda ia ea 80 
10. US8 Serii de puteri. Dezvoltarea unei funcții în serie de puteri... nenea 107 
111059 Puneţii:demai multe variabile sa si onest pana oa ta aaa aaa tea pla aaa iata 115 
12 BIbII0giă ÎNC sate ae al ou n ua ne aaa La a dai gt a 134 


USI Transformări elementare 


Nenumărate probleme din matematică utilizează noţiunea de matrice iar 
transformările elementare de linii aplicate acestora pot înlocui calculele lungi şi chinuitoare pe 
care unele metode algebrice le implică. Ne propunem în această unitate de studiu să prezentăm 
transformările elementare de linii într-o manieră practică şi unitară, şi insistăm asupra aplicaţiilor 
lor în determinarea rangului, în calculul determinantului şi a inversei unei matrici. De asemenea 
aplicăm transformările elementare pentru discuția compatibilității şi determinarea soluţiilor 
sistemelor de ecuații liniare. Putem astfel să evităm rezultate clasice (teoremele lui Rouché, 
Kronecker — Capelli şi Cramer) şi metode complicate care implică prea multe calcule. 

Studentul trebuie să se obișnuiască mai întâi să lucreze cu liniile unei matrici, 
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după care trebuie să îşi stabiliească obiectivul: formarea pe coloana „i? a coloanei „k” din 
matricea unitate I„. În continuare trebuie stabilite operaţiile necesare scopului propus şi 
efectuarea calculelor implicate de acestea. 

Ulterior, studentul aplică efectiv metoda însușşită pentru determinarea rangului sau 
inverse unei matrice sau pentru discuţia şi rezolvarea unui sistem de ecuaţii algebrice liniare. 

Timpul necesar însuşirii materialului conținut este de aproximativ 6 ore (2 ore 
pentru învățarea metodei, 2 ore pentru determinarea rangului şi a inversei, şi două ore pentru 
aplicaţii la discuţia sistemelor). 

Metoda este prezentată pe larg în materialul nostru dar şi în lista bibliografică la 


poziţiile [3], [4], [5] şi [6]. 


x x 


Studiul matricelor este legat de discuția şi rezolvarea sistemelor de ecuații algebrice 
liniare şi constituie elementul de control al aproape întregii programe analitice a disciplinei 
„Matematici aplicate în economie”. Noutatea pe care o propunem se numeşte transformare 
elementară efectuată asupra liniilor unei matrice, obţinând din aceasta o matrice echivalentă din 


punctul de vedere al rangului. Departe de a fi o abstractizare inutilă, folosirea transformărilor 


elementare poate duce, fără calcule dificile, la determinarea rangului, la obţinerea inversei, sau, 
la obţinerea soluţiei unui sistem algebric liniar. 

Se cunosc dificultățile pe care le întâmpinăm la discuţia sau rezolvarea acestui tip de 
sisteme. Metodele bazate pe teoremele lui Cramer, Kroneker-Capelli sau Rouche sunt greoaie şi 
folosesc determinanţii, al căror calcul este laborios. Prin transformarea matricei lărgite a unui 
astfel de sistem se obţine forma explicită din care se pot trage concluziile privind tipul şi soluţiile 
sale. O categorie specială de soluţii ale unui sistem liniar compatibil nedeterminat o formează 


soluţiile de bază care vor juca un rol deosebit în rezolvarea problemelor de programare liniară. 


1.1 MATRICE 


1.1.1 RANGUL UNEI MATRICE 


Prin matrice înţelegem o aplicaţie A:1x J—K, unde I= 1,2,...,m ; J= Lun , K O 


mulţime oarecare. Ea poate fi reprezentată printr-un tablou de elemente din K, aşezate pe linii şi 


coloane. 


Pentru o matrice de tipul m x n vom folosi notația: 


Ai diz .... Ain 
A= da1 dz ..... dn (1.1.1) 
di Uyy * a amn 


Suprimând din matricea A, m — r linii şi n — r coloane, obținem o matrice pătratică de 


ordinul r. Determinantul acestei matrice se numeşte minor de ordinul r al matricei. 


Definiția 1.1.1. Numim rangul unei matrice A de tipul m x n (m n ) un număr natural r 
cu proprietăţile: 
a) în matrice există cel puţin un minor de ordinul r diferit de zero. 


b) toți minorii de ordinul r + 1 sunt nuli. 


Observaţie. Se demonstrează că dacă toţi minorii de ordinul r + 1 sunt nuli, atunci sunt 


nuli şi toţi minorii de ordin mai mare ca r + 1. 


In continuare vom pune în evidenţă operaţii pe care le vom numi transformări elementare, 
care se bucură de proprietatea că nu modifică rangul unei matrice. Ele au la bază proprietăţile 


determinanţilor. 
1.1.2 TRANSFORMĂRI ELEMENTARE 


Definiţia 1.1.2. Numim transformare elementară aplicată unei matrice una din 
următoarele operaţii: 
(Ty) — înmulţirea unei linii cu un scalar real nenul (az0); 


(T) — adunarea unei linii la o altă linie; 


(13) — schimbarea a două linii între ele. 


Teorema 1.1.1. Transformările elementare nu modifică rangul unei matrice. 

Demonstraţie. Fie matricea A de tipul m x n şi matricele Ti(; =1,3) care se obțin din 
matricea unitate Im. Înmulţind la stânga matricea A cu matricele T, obţinem matricele 4, = TA. 
Presupunem că rangul matricei A este r,.. Vom avea: 

TA, < min Cr, TA) = min (m, r4) =f} (1.1.2) 

Din A; = TA, cum T; sunt ineversabile găsim A =T; "A, şi deci: 


Ta. < min (r,r, ) = min (mr, STA. 1.1.3 
A; T, fA A; A; 


Din (1.1.2) şi (1.1.3) rezultă că r4, = ra, 1=1,3. 


Definiția (1.1.3). Două matrice A şi B care se obțin una din alta prin transformări 


elementare se numesc echivalente ( în privința rangului) şi scriem A ~B. 


1.1.3. FORMA GAUSS — JORDAN A UNEI MATRICE 


Fie A matricea unui sistem de m ecuații algebrice liniare cu n necunoscute (m < n). 


Definiția 1.1.4. Matricea A se spune că are forma Gauss-Jordan dacă conține r (r <m) 


coloane ale matricei unitate de ordinul m. 


T EER i EER OS 00 aaa d a (1.1.4) 
0.. 0...0 ....0 0...0 
0..0....0....0 0...0 


Teorema 1.1.2. Orice matrice nenulă poate fi adusă la forma Gauss-Jordan, printr- 


un număr finit de transformări elementare. 


Demonstraţie. Fie deci matricea A de tipul m x n a unui sistem liniar de m ecuaţii cu n 


necunoscute. 
d ah dij din 
do n dij don 
di din Aij Ain 
äp aya se Agg esr Min 


Fie a, #0. Dacă a, = 0, atunci printr-o transformare de tipul (T3) putem aduce în locul 


lui a, un element a; #0; i=2,m care există deoarece prima coloană are cel puțin un element 


diferit de zero, în caz contrar sistemul nu ar avea n necunoscute. 


Printr-o transformare de tipul (T1) cu a = 23 realizăm în poziția (1,1) un element egal 
dii 


cu unu. Prin m-—l transformări de tipul (T1) şi (T2) în care linia întâi este succesiv inmulțită cu 
numere convenabile: —a,;, —431; —i5---5 — şi adunată de fiecare dată la liniile: a doua, a treia, 
.. a i-a linie, ... a m-a linie, realizăm pe prima coloană toate elementele nule cu excepția 
primului element care este unu. 
Dacă dorim să realizăm pe coloana j toate elementele nule cu excepția elementului din 
poziția (i,j) care să fie egal cu unu procedăm astfel: 


a) ne asigurăm că a; #0 
. l l da 
b) aplicăm o transformare de tipul (Tı) cu a=—, obţinând în poziţia (ij) un element 


egal cu unu. 
c) Aplicăm succesiv m-—l transformări de tipul (Tı) şi (T2), înmulţinâd linia į cu: 


Go td ȘI adunând-o la liniile: întâia, a doua, . .. , a m-a linie obţinem 
pe coloana j toate elementele nule cu excepţia elementului din poziţia (ij) care devine 
unu. Algoritmul se opreşte când pe liniile r+1,m nu mai avem elemente nenule. 


Elementul nenul aij dat se va numi pivot, iar ansamblul transformărilor elementare 
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necesare transformării matricei A într-o matrice cu coloana i din matricea unitate de 
ordinul m în coloana j se va numi pivotajul cu element pivot a;. 

Din forma Gauss-Jordan a matricei A putem observa că rangul matricei este r, acesta 
fiind de ordinul minorului format cu primele r linii şi coloanele /,/,,...j,. Toţi minorii de 
ordinul r+1 sunt nuli. În practică nu este neapărat necesar să obţinem coloanele matricei unitate 
în ordinea lor naturală, numărul lor ne va da rangul. Putem alege pivot orice element diferit de 


zero oriunde s-ar găsi el de preferinţă cel care este egal cu 1. 


Exemplul 1.1.1. Să se aducă la forma Gauss-Jordan şi să se determine rangul matricei: 


O 222 1l Luăm pivot elementul încercuit. Prima linie 
a cert să înmulțită cu 1 se adună la a doua, apoi cu —2 şi o 
2 4 4 3 1 , e À 
adunăm la a treia obținând matricea : 
2 1 Luăm pivot elementul încercuit. Linia a doua o 


1 2 2 
009 3 -1 înmulțim cu 1/3 şi apoi înmulțită cu —2 o adunăm 
0 0 0 


—1 1l la prima obţinând matricea: 


120 0 5/3 Luăm pivot elementul încercuit. Linia a treia o 
0 0 1 1I —1/⁄3 înmulţim cu —1, apoi înmulțită cu —1 o adunăm la 
020; :0: = l a doua şi obținem: 
paz D Rangul matricei este egal cu 3, egal cu 
0 1 0 243 numărul maxim de coloane ale matricei 
0 0 1 -l1 unitate. 
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Exemplul 1.1.2. Să se determine rangul matricei: 


2 1 -l 3 
-1 —2 1 —2 
1 —1 2 -l 
Făcând pivotajele cu elementele pivot încercuite obținem succesiv: 
2 © -l1 3| 2 1 —1 3 
—]1 —2 1 —2|[3 0 —1 4|» 
1 —l 2 -1| 3 0 Ọ© 2 
5 1 0 5 0 10 0 1 00 1 


~®© 0 0 6[=[1 00 1j~lO 1 0 0 
3 0 1 2| 0 O I —Ii| |0 O 1 —l 


Rangul matricei considerate este 3. 


1.1.4 INVERSA UNEI MATRICE 


Fie A o matrice pătrată de ordinul m şi B= A : 1, matricea obținută prin alăturarea 


matricei unitate de ordinul m la matricea A. Se cunoaşte că se numeşte matrice inversă a unei 


matrice pătrate A, o matrice, notate cu A`”, astfel încât: 
AA'=AA=I,, (1.1.6) 
Presupunem pentru moment că există matricea A”. Înmulţind la stânga matricea B cu A” 


obţinem: 
B=A7'B=-|A 'A A Ta =a A 


De la matricea B putem ajunge la matricea B cu ajutorul transformărilor elementare. 
Făcând pivotaje în matricea A în vederea obținerii celor m coloane ale matricei unitate de ordinul 
A : Si ` ; -1 X A 3 S 

m, în dreapta barei vom găsi tocmai matricea A” . Dacă coloanele nu sunt în ordinea lor naturală, 
putem schimba liniile între ele obținând în stânga barei chiar matricea unitate. Dacă nu reuşim să 
obținem coloane ale matricei unitate în număr egal cu ordinul matricei A, atunci matricea nu are 


inversă. 
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Exemplul 1.1.3. Să se găsească inversa matricei: 


1 2 3 
A=|—1 —l —2 
1 1 | 


considerând în matricea B pivotaje cu elemente pivot încercuite avem : 


© 2 3 1 00 [1 2 3 00 0 
Beza! -1 —2 : 0 1 0~Ọ00 1 © 1 1 0x 
1 1 1:00 1| 0 -1 -2 -1 0 1 
1 -1 0 | —2 -3 0 h 00: -1 -1 1 
~o 11 1 1 0|=l0 0 1 0 e - ES 
0 0 0 1 21| 010 l1 2 1 

1 0 0: Aa 1 

=l0 1 0: 1 2 1 

001: 0 -1 -l 


Inversa matricei A este: 


i E, | 
Al 1 2 1 
~=] 


Exemplul 1.1.4. Să se găsească inversa matricei: 


e 29 
A= 2 =] 
| 
Vom avea: 
© -1 2:1 00 [i -12 1 00 
B=|2 2 -1:0 1 0~l0 ® 25 © 1 ül 
3 1 1:001 04 —5 -3 0 1 
0 3/4: 12 14 0 
~l0 1 —5/4 î —1/⁄2 1⁄4 0 


00 0 : -l —1 1 


Matricea nu admite inversă. 
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1.2. SISTEME DE ECUAŢII LINIARE 


Un sistem de ecuaţii algebrice liniare cu m ecuaţii şi n necunoscute 
(m < n) se prezintă sub forma: 


Cm + apm2X +--+ AmnXn = bmo 


unde a,b; (=1,m ; a =1,n) sunt numere reale. 


Dacă există b; 20 sistemul (1.2.1) este neomogen, în caz contrar (b, =0,(v)i =1,m) 


sistemul se numeşte omogen. 


Fie B matricea lărgită a sistemului (1.2.1), de forma: 


ai a2 an i b, 
ec a o Caa 
ami Am2 amn bn 


Să presupunem că rangul matricei A a sistemului este r (r <m). Prin transformări 


elementare putem aduce matricea A la forma Gauss-Jordan. Obţinem astfel matricea B, de 


forma: 
l ag O ae 0 aia T ai 
0 cs ol său 308 Ana 77000 b 
BAO aste Doza E Sua mac Zi be (1.2.3) 
0... 0 0 0 a ba 
0 0 0 0 0 $b, 


Sistemul corespunzător matricei B este: 
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Xi + ETER E TT Fajara Fasai n ==, bi 
Xi F. To aN T.. Han n =b, 


Xp Fă pune Foo ma Sb, (1.2.4) 
0 =b, 
0 =b, 


Deoarece matricele B şi B sunt echivalente, sistemele (1.2.1) şi (1.2.4) vor fi echivalente 


(adică vor avea aceleaşi soluţii). Este suficient pentru aceasta să observăm că transformările 


elementare aplicate liniilor matricei B pentru a fi adusă la forma B au drept consecinţe asupra 


sistemului de ecuaţii: 
(Tı) — înmulţirea unei ecuaţii cu un număr diferit de zero; 


(T2) — adunarea la o ecuaţie a unei alte ecuaţii (eventual înmulțită cu un scalar nenul); 


(T3) — schimbarea a două ecuaţii între ele. 


Evident, aceste operaţii aplicate unui sistem de ecuaţii îl transformă în unul echivalent. 


Teorema 1.2.1. Condiţia necesară şi suficientă ca sistemul (1.2.1) să fie compatibil 


este ca b, =0, i=r+l,m. 


Demonstrație. Să arătăm mai întâi necesitatea condiției. 


Presupunem că sistemul (1.2.1) este compatibil. Aceasta înseamnă că sistemul (1.2.4) 


este compatibil. Din compatibilitatea acestuia din urmă rezultă b, = 0 pentru i =r +1,m. 


Să arătăm acum suficiența condiției. Presupunem că b, =0, i=r+l,m. Rezultă atunci că 


sistemul (1.2.4) este compatibil admițând co” ” soluții și deci sistemul (1.2.1) este 
compatibil. 
În baza acestei teoreme uşor pot fi demonstrate teoremele lui Kronnecker-Capelli şi 


Rouché. Pentru r=m=n se obțin rezultatele cunoscute de la sistemele de n ecuații cu n 


necunoscute. 
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Exemplul 1.2.1. Să se atudieze compatibilitatea sistemului şi în caz afirmativ să se 


găsească soluţia lui. 


It = Xx + X4 = 
de = Xx +2% —2x4, = 
3x4- xn +3% —3x, = 5 
=x +3x, — 5x, +5x4 =—3 


Matricea lărgită a sistemului, după o serie de pivotaje cu elemente pivot încercuite, 


devine: 
© 1—1 1; I 
gaine 2 -2 2l. 
KEI 3a 5 
So e 3) 20] 
1 1—1 1; 1| [1 0 1/2 —1/2: 3/2] 
p- g | NO aa 3 N 
“lo -4 6-6 2 b0 0o o o 
0 4-6 6 -—2| |0 0 0 0; 0] 


Sistemul considerat este echivalent cu sistemul: 


x +1/2x, — 1/2x, = 3/2 
x — 3/2% +3/2x;, =—1/2' 


care este compatibil având soluția: 


x =3/2-1/2æ+1/2ß, 
x, =—1/2+3/2a—3/2B, 


X=0, X,=B, a,peR 


Exemplul 1.2.2. Să se studieze compatibilitatea sistemului: 


x + 2X —3X3 +4x; = 2 


XA 
2x +5x A 


X + X3 
6x3 - 
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Matricea lărgită a sistemului după o serie de pivotaje cu elemente pivot încercuite devine: 


0.28 A 2) [4 2 a 03 2478 


B=|-1 1 1 —1: 3|=]0 © 4 3 5|-l0lI 4/3 1; 5/3 
-3 5 6 1.10) [0 9 12 9:14] 00 oo: ~ 
Cum b, =—170, sistemul este incompatibil. 


1.2.1. EXPLICITAREA SISTEMELOR 
DE ECUAȚII LINIARE 


Fie sistemul (1.2.1). Presupunem că rangul matricei A este m < n. 


Definiţia (1.2.1). Vom spune că sistemul (1.2.1) este explicitat în raport cu un grup de m 
variabile (necunoscute), dacă în matricea A a sistemului coloanele acestor variabile sunt cele m 
coloane ale matricei unitate de ordinul m. 

Deoarece cu cele n variabile ale sistemului putem forma Cs grupuri diferite de câte m 
variabile, rezultă că un sistem liniar, care poate fi explicitat în raport cu cel puţin un grup de m 


variabile, va avea cel mult Cp forme explicite. Acest număr maxim va fi atins dacă explicitarea 


poate avea loc în raport cu oricare dintre cele CF grupuri diferite de câte n variabile şi toate 


formele explicite sunt distincte. 

Într-un sistem liniar explicitat sunt două tipuri de variabile (necunoscute): unele, ai căror 
coeficienţi formează coloanele matricei unitate, altele ai căror coeficienţi formează celelalte 
coloane ale matricei sistemului. Ca de obicei, primele vor fi denumite variabile principale, iar 
celelalte variabile secundare. În teoria programării liniare variabilele principale mai sunt 


denumite bazice, iar variabilele secundare nebazice. 


Definiţia (1.2.2). Numim soluţie de bază a unui sistem de ecuaţii liniare, o soluție 


particulară obținută dintr-o formă explicită, prin egalarea cu zero a variabilelor secundare. 
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Din definiţia dată rezultă că dacă sistemul este compatibil nedeterminat, atunci are cel 
puţin o soluţie de bază şi cel mult Ci soluţii de bază. Evident, sistemele incompatibile şi cele 


compatibil determinate nu au soluţii de bază. 

Dacă o soluţie de bază are exact m componente nenule, atunci soluţia de bază se numeşte 
nedegenerată, iar dacă are mai puţin de m componente nenule se numeşte soluţie de bază 
degenerată. 

Soluţiile de bază ale unui sistem liniar se clasifică şi după un alt criteriu. Dacă toate 
componentele unei soluţii de bază au valori nenegative (20), soluția de bază se numeşte 
admisibilă. În caz contrar, soluţia de bază se numeşte neadmisibilă (are componente strict 
negative). 


Să considerăm că r = m. În acest caz din (1.2.4) obţinem o soluţie de bază: 


xX, =b, xX; Doo see X; =b.. X; =bn, restul până la n egale cu zero. Dorim ca variabila 
principală x, să devină secundară iar variabila secundară x,(jz ji, ...;j,) să devină 


principală. Pentru aceasta fie şi a; 20. Efectuând în matricea B un pivotaj cu element pivot 


a; găsim o altă soluţie de bază: 


II. Ea ae -i b, ci ce 

=b aj Xp 5b aj eee So ne Xyp Dn mi 

ay a dij ij 

şi restul până la n egale cu zero. 
Putem scrie variabilele principale concentrat sub forma: 
_ b, | 
aan „ pentru kz] 
5 y (1.2.5) 

ER =e 
A „ pentru k=] 


Exemplul 1.2.3. Să se determine toate formele explicite şi soluţiile de bază 


corespunzătoare pentru sistemul: 


—X, + 2 — Xa Fa =2 


Matricea lărgită a sistemului este: 
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Pivotajele: 


dacat no ozb tona 


1 
i 2a 12o 1002/1012 


AESA 


S itrle e 
a o ale -1 —1 0 0] [0 1 41 2 
c : 
= “2 KLI IEO a 
TRE 1 1. 2030). [le = 1000 
=i 1002 lo 1012 


k 1 010], 
e) | 


A MN 


= —1 010] |-1 1 —1 0:0 
to- 4 02 10 112 


conduc la formele explicite în raport cu grupurile: 


2 ), X X4), 2, 23), (%2, X4), 03 X4)» 


a) Xz +x, =2 5 Xz +x, =2 
c) Xz Fo =2 d) —X TX — Xa = 0 
x +t =2 X e Paget = 


P x X% +x = 0 
—X, +2x, + Xa =2 


Soluţiile de bază vor fi: 
(S) x, =—2, x =2, x =0, x, =0; soluţie de bază nedegenerată admisibilă 
(5) x, —2, x =0, x =0, x, —2; soluție de bază nedegenerată admisibilă 
(S) x, =0, x =2, x =2, x, =0; soluţie de bază nedegenerată admisibilă 
(S4) x, =0, x, =0, x, =0, x, =2; soluţie de bază degenerată admisibilă 
(S) x, =0, x, =0, x =0, x4 =2; soluţie de bază degenerată admisibilă 
Sistemul nu poate fi explicitat în raport cu grupul de variabile (xx) şi admite doar 


je vie 


patru soluţii de bază distincte din numărul maxim de C = 6 posibilităţi. 
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US2 Spaţii liniare 


Noţiunea de vector şi cea de spaţiu liniar sunt noțiuni centrale în algebra liniară. 
Utilizarea acestora pentru rezolvarea problemelor de optimizare este cât se poate de naturală. In 


prezenta unitate de studiu se prezintă noţiunea de spaţiu liniar cu exemplul cel mai potrivit (adică 


spaţiul R”) dar şi cu alte exemple. De asemenea se aplică elemente din USI (Transformări 
elementare) pentru a stabili relația de dependenţă (independenţă) liniară a unor familii de vectori. 
De asemenea se utilizează sistemele liniare (abordate şi în USI) pentru a stabili concret relaţiile 
de dependență liniară. Se introduce în continuare noțiunea de bază de vectori, precum şi de 
coordonate ale unui vector într-o bază. Schimbarea de coordonate, atunci când se schimbă 
vectorii bazei, este următorul subiect abordat în material. 

Cei care şi-au însușit complet USI nu vor avea dificultăţi pentru parcurgerea 
acestei unități de studiu. Mai întâi trebuie însușite corect noțiunile de bază (vector, bază, 
coordonate). Ulterior trebuie aplicate teoremele de caracterizare a dependenței și independenţei 
(folosind efectiv transformările elementare). Pentru determinarea coordonatelor unui vector se 
vor rezolva (cu transformări elementare) sisteme de ecuații compatibile determinate. 

Însuşirea noţiunilor teoretice presupune un studiu de patru ore. Aplicaţii pentru 
dependenţă şi independență solicită încă două ore. Deteminarea coordonatelor (rezolvarea 
sistemelor liniare) cu aplicaţiile corespunzătoare mai adaugă încă două ore la timpul de studiu. În 
capiutolul acesta se prezintă un algoritm de calcul a noilor coordonate ale unui vector când unul 
dintre vectorii bazei este înlocuit. Însuşirea lui şi aplicarea cu uşurinţă înseamnă încă două ore de 
studiu. În total US2 necesită zece ore de studiu. 


Pentru facilitarea studiului sunt utile din bibliografie pozițiile [5] şi [6]. 


x 
x x 
Înzestrând o mulţime oarecare cu o serie de operații având anumite proprietăți obţinem 
noţinea de structură algebrică. Unele tipuri de structuri algebrice au fost studiate în liceu 
(monoid, grup, inel, corp). Spaţiul liniar se defineşte în acelaşi mod, elementele sale (vectorii) 


generalizând noțiunea de vector liber utilizat în fizică sau matematică. 
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În orice spaţiu liniar se introduce noţiunea de independenţă liniară precum şi cea de 
bază. Vom introduce aceleaşi noţiuni dar numai în cazul particular al spaţiului liniar R” (care 
apare cel mai des în contextul aplicaţiilor economice). 

În studiul problemelor economice se impune problema determinării valorii maxime sau a 
valorii minime a unor funcţii. Un ajutor preţios în acest scop îl constituie formele pătratice care 


generalizează funcția de gradul al doilea. 
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2.1. DEFINIŢIA SPAȚIULUI LINIAR 


Fie A o mulţime oarecare şi K un corp având elementele neutre Ok respectiv Ik. 
Presupunem că pe A se pot defini două operaţii astfel: 

a) (V)x,y € A, (3) un element notat cu x+y € A astfel încât (x,y) —> X +y. 

b) (v)xeaA, (V) XEK, (Ə) un element notat A x € A , astfel încât (A,x) — A x. 

Prima operație, de tip aditiv, este o operație internă iar cea de doua, de tip multiplicativ, 


este o operație externă pentru A. 


Definiție. Mulțimea A formează un spațiu liniar peste corpul K dacă sunt satisfăcute 
următoarele proprietăți: 

Lı) (V)X, yY, ZEA, x+y +z)=(x+y)+z, 

L2) (V) x, YEA, x+y =y+X, 

L) (3)0, EA astfel încât x+0, =x, (YV) XEA, 

La) (Y) x€A, (IÐ —x EA astfel încât x+(—x)=0,, 

Ls) (V) AEK, V) xy EA, Ay) =X +Y, 

Ls) V) AEK, V) pEK, (WxeA,  A+p)x=X py, 

L) VNp EK, (Vx EA, Nx) = LAX) = (Ax, 

Ls) (VX EA, lkXx=X. 


Elementele lui K se numesc scalari, iar cele ale lui A, vectori. Prima operație se numeşte 
adunarea vectorilor iar cea de a doua înmulţirea vectorilor cu scalari din K. Atunci când nu 
este pericol de confuzie vectorul nul 04 se notează simplu 0. De asemenea, elementele neutre ale 
lui K se scriu 0 şi 1 (în imensa majoritate a aplicaţiilor K=R ). 

Consecințe. Din axiome decurg următoarele proprietăţi: 

a) (V)xeA, 0x=0, 

Într-adevăr, din Le) obţinem (A +0) x = Ax +0x, adică Ax = Xx + 0x. 

Axioma L; ) conduce la 0x =0. 

b) V)AEK,  X0=0. 

Plecăm de la Ls): X(x+0) = x +0. Obţinem Ax = A +0. 
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Aceeaşi axiomă L3) arată că: X0=0. 
c) (v)xeA4, (—l)x=—x 

Folosim a), adică 0x =0 şi înlocuim 0 prin 1+(—1). Din Le) obţinem: 
1x+(—Dx = 0. Lg) conduce la x+(—1)x =0, iar L4) arată că (—1)x =0. 
Exemplu: 

Considerăm A=R"=RxRx...R= x=(%,X%,...,X%,)/ Xx; ER şi KR. 
Adunarea elementelor din A revine la adunarea matricelor cu o singură coloană sau cu o 
singură linie. 

(VX S (Orase Y =i ase) din R”, 
x+y =% +Y X HV... Xp PY) ER” 
În mod similar produsul unui scalar din R cu un element al lui R” se defineşte 


prin: XX = NX sta se ta Xp) =X Aiseee AX). 


n 


2.2. DEPENDENȚĂ ŞI INDEPENDENȚĂ 
LINIARĂ 


Fie A un spațiu liniar peste corpul Ķ şi fie x; €A, EK, i= lh. 


Definiție 2.2.1. Se numeşte combinaţie liniară a vectorilor x, cu scalarii X, expresia 


E = NX HAX H.N. (2.2.1) 
Consecința a) din definiția spațiului liniar ne arată că atunci când 


A =0,(v) i= l,h => E =0. Reciproca acesteia nu este întotdeauna adevărată. 


Definiție 2.2.2. Vectorii x; = A, i= 1,h se numesc liniar independenţi dacă: 


At Apt TA 0 > ANA 0. 


In caz contrar, dacă E=0 cu măcar un scalar nenul, vectorii x;se numesc liniari 
dependenţi, adică: 
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= 


Exemplu: 


Fie spaţiul liniar A = R” = M,,(R) şi vectorii: 


t t t $ S 
X = (Xp X25: Xn) X2 = (X1 X225- X2) 3 Xp = (X1 Xp2>:--Xpn) Vectori din A. 
Condiţia E = 0 conduce la sistemul liniar şi omogen: 


N HAX ee Am =0 
A2 FAX ek Apa > 0 (2.2.3) 


Vectorii sunt liniari independenţi dacă sistemul (2.2.3) admite numai soluţia banală şi liniar 
dependenţi în caz contrar. 


Din teoria sistemelor omogene independenţa liniară este caracterizată prin: 
rang M=r,=h (2.2.4) 
unde matricea M € Maa (R) este formată prin componentele vectorilor x,. 
Dacă (2.2.3) are soluţii nebanale, ele se exprimă prin h —rņų parametri reali. Înlocuite în 


(2.2.2) aceste soluţii conduc la una sau mai multe relaţii de dependenţă liniară. 
Pe baza acestui exemplu vom realiza un mare număr de aplicaţii. Proprietăţile familiilor 
de vectori liniar independenţi şi al familiilor de vectori liniar dependenţi sunt date în următoarele 


teoreme: 


Teorema 2.2.1. O mulţime de vectori liniar independenţi nu poate conţine vectorul nul. 


Demonstraţie: Presupunem prin absurd că mulţimea de vectori independenţi 


F= Xp Xa Xio-::oXp ar conţine vectorul nul şi că x, = 0. Combinația liniară : 


E =—0-x, +0-x, +...+1-x; +...-+A,X, =0 cu scalarul A, nenul. Contradicția justifică 


teorema. 


Teorema 2.2.2. Orice submulțime a unei mulţimi de vectori liniar independenţi este formată tot din vectori liniar independenţi. 
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Demonstraţie: Tot prin reducere la absurd presupunem că deşi x,,xp,....X,, sunt liniar 


independenţi, submulţimea X; > Xp eX, o KSA este formată din vectori liniar dependenți. 


Atunci, egalitatea: 
A X FAX Fe FĂ, =0 


LA H bb Ik Ik 
are loc cu măcar un scalar nenul. 
Rezultă ca: Aa X; TA Xi, + PRT + N, x, + Žin 0 j X; — 0 
JE 
cu măcar un scalar nenul ceea ce contrazice ipoteza teoremei. 


Teorema 2.2.3. Într-o mulţime de vectori liniar dependenţi cel puţin un vector se 


scrie ca o combinaţie liniară a celorlalţi. 


Demonstraţie: Definiţia vectorilor liniar dependenţi arată că: 
AXT ADA =0 
cu măcar un scalar nenul. Presupunem că acesta este A;. Cum orice spaţiu liniar este grup 
abelian, regulile de calcul dintr-o astfel de structură permite să scriem: 


AX; = —A X] — Xp occ..c..... —AĂ —Nuăi-u oc... — A,X; 


l 
În corpul K, X, este inversabil iar dacă >! este inversul său 
=i -1 -1 -1 
X; — —N AX T == A; AX -N Nouă iosi —N Xp (2.2.5) 
şi teorema este demonstrată. 
Când corpul de scalari este R inversul X;! al lui >, este 1/X, şi obținem: 


X; =- -3x i X; Siy ee bg (2.2.6) 


i i i i i 


Teoremele prezentate au consecințe importante cu valoare teoretică şi practică. Astfel, 
modificarea componenței unei mulțimi de vectori liniar independenți prin suprimarea unora 
conduc totdeauna la o familie de vectori liniar independenți. Dacă adăugăm noi vectori unei 
astfel de familii, rezultatul nu este previzibil, noua mulțime poate fi dependentă sau 


independentă. Din (2.2.3) rezultă că (2.2.4) nu poate fi realizată dacă h>n , deci numărul maxim 


de vectori liniar independenți din R” este n. 
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Există efectiv mulţimi de vectori liniar independenţi. Una dintre acestea este: 
F= 0 pasate (2.2.7) 
unde: e; = (1,0,0,...,0)}, e, =(0,1,0,...,0},... e, =(0,0,0,..., D". 
Dacă prin transformări elementare modificăm matricea M a vectorilor mulțimii F , 


obținem alte mulțimi formate tot din n vectori liniar independenți. Numărul acestor mulțimi este 


infinit. 


2.3. BAZE DE VECTORI. COORDONATELE UNUI 
VECTOR ÎNTR-O BAZĂ 


Am luat în studiu în paragraful precedent familiile maximale de vectori liniar 


independenți dintr-un spațiu liniar. Rolul acestora va fi precizat în cele ce urmează. 


Definiție. Familia B= u,,u,,...,U, formează o bază în spațiul vectorial V dacă: 
a) U,u,,... U, sunt liniar independenți; 
b) (WueV, mulțimea F= u,uj,u,,..-.u, , este formată din vectori liniar 


dependenți. 


Mulțimea B este aşadar o mulțime maximală de vectori liniar independenți. Când spațiul 
V este R”, această mulțime conține n vectori, deci k = n. 

Spaţiul R? are cel puțin o bază: cea formată din e,e,,...„e, numită bază canonică. 
Aşadar: 

B, = €,€,...sen 

Din aceasta prin procedeul indicat anterior putem obţine o infinitate de baze toate având n 


vectori. Numărul n, al vectorilor oricărei baze din R” se numeşte dimensiunea acestui spaţiu. 
Condiţia b) din definiţia bazei şi una din proprietățile fundamentale ale unei familii de 


vectori liniar dependenţi ne permit să formulăm şi să demonstrăm: 
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Teorema 2.3.1. Orice vector din R” se exprimă în mod unic printr-o combinaţie 


liniară a vectorilor unei baze B. 


Demonstraţie. Teorema 2.2.3 ne arată că cel puţin un vector al mulțimii 


U,U,U,,...,u, în care u,,u,,....u, formează bază în R” se exprimă printr-o combinaţie 


liniară a celorlalți. Acest lucru este valabil şi pentru u, deoarece în caz contrar, din: 
`u AU, +u, Fiu, =0 (2.3.1) 
în care A =0 ar rezulta: 
NU, HAU, +AU =0 (2.3.2) 
şi unul din scalari fiind nenul s-ar contrazice condiţia a) din definiţia bazei. 


Aşadar A #0 şi obținem: 


N » A 
u= u u Zu 2.3.3 
` 1 x 2 ` n ( ) 
A A Naben as 
Notăm: =L m =... u, = gi scriem: 
H ` H2 N Hn X Ş 
U = PU, to et nUn (2.3.4) 


Prin reducere la absurd demonstrăm că scalarii p; din (2.3.4) sunt unic determinaţi. In 
caz contrar, are loc şi relaţia: 


su, (2.3.5) 


Din (2.3.4)-(2.3.5), obţinem: 

(h — v)u; + (p2 —V2)uz +... (Ha — Vadu =0 (2.3.6) 

Cum u,,u,,...,u, sunt liniar independenţi, toți scalarii din (2.3.6) sunt nuli, deci: 
Hi 5 Vi, M2 Matei fie = Vn 


ceea ce este în contradicție cu ipoteza. Teorema este demonstrată. 
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Observaţie. Scalarii p;, (i =1,n) se numesc în mod uzual coordonatele vectorului u în 


baza B= u,,u,...,u, . Unicitatea lor rezultă şi din teoremă dar şi din faptul că (2.3.4) 


n 


conduce la un sistem Cramer (cu soluţie unică). 


Vom nota mulțimea coordonatelor vectorului u în baza B cu up, adică: 


Ug = (His Has- - -Hn )" (2.3.7) 


Un calcul elementar ne arată că în baza canonică B, = €,.€,.....e, coordonatele 


n 


oricărui vector u coincid cu componentele sale, adică: 


t 
u = (uj, Uz5... U) = 4e truse, F:::+use, 


Exemplu: 
Să se arate că vectorii u, = (1,1,1), u, = (1,2,1)t, u, = (1,1,3)! formează o bază în R? şi 


să se determine coordonatele lui u = (6,1,8)! în această bază. 


Relația (2.3.4) devine: 
u = pU; pU + puz 


şi conduce la sistemul: 


u+ H+ H3=6 
u +2 > 
u+ M +3 =8 
Prin transformări elementare de linii în matricea lărgită a sistemului obținem: 


1 1 16) (11 1; 6) {10 1; 11) {10 0}10 
1 2 1i1ij~|0 1 0:—5|=|0 1 0:—5|-|0 1 0:—5 
1 1 3:8] (0 0 2; 2} [o 02; 2] 0 01; 1 


care corespunde sistemului cu r(M)=3 ) deci (u, u2, uz formează baza). 
m =10 
u = —5 
ua =l 


acestea fiind valorile coordonatelor căutate. 
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2.4. SCHIMBAREA COORDONATELOR UNUI 
VECTOR LA SCHIMBAREA BAZEI 


Dacă în locul bazei B= u,,u,...,u, a spaţiului liniar R” considerăm o altă bază 
B' = Vis V>>:-: Va având acelaşi număr de vectori, n, ne aşteptăm ca în relația (2.3.4) 


coordonatele să se schimbe. Dorim să aflăm legea de schimbare a acestora în cazul general şi 


într-un caz particular. 


Să considerăm că în baza B (numită baza veche) vectorul u are coordonatele 
4s PD Qp p deci: 

U = quU +U, ++ aU, (2.4.1) 
şi că acelaşi vector are în baza B’ (numită baza nouă) coordonatele f3,,3,...8, , adică: 

u = Bv, +Bova +... +BaVa (2.4.2) 

Pe de altă parte fiecare din vectorii v; ai bazei noi are în baza veche coordonatele a;;, 


unic determinate 
n 
Vi = AU + au, osul = 2 aju; (2.4.3) 
j=l 


Înlocuim (2.4.3) în (2.4.2) şi comparăm cu (2.4.1). 
Obţinem: 
Qi = aai FaioBa +: Faima i=ln (2.4.4) 
Dacă notăm cu A matricea coordonatelor vectorilor din baza nouă faţă de vectorii din 


baza veche (numită matrice de trecere): 


CU Up (0503: ) şi cu Uy = (BiB2> Ba )' atunci (2.2.4) se scrie matricial : 


up =A'uy, (2.4.5) 
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Cum vectorii lui B’ sunt liniar independenţi, A este inversabilă şi din (2.4.5) putem obţine 
relaţia: 


—1 


uy= A! w (2.4.6) 


Relaţiile (2.4.5) şi (2.4.6) sunt formulele de schimbare a coordonatelor la schimbarea 


bazei. 


Exemplu: 
Se consideră vectorii u, = (1,2,1), u, = (2,1,0), u, =(1,10) din R°: 
a) Să se arate că formează o bază în R? 
b) Să se determine coordonatele lui u = (1,1,—2)' în această bază. 


c) Să se determine coordonatele lui u în baza formată din 


Soluţie: 
a) Matricea formată cu componentele vectorilor u,.,u;,uz are rangul trei, egal cu 


dimensiunea lui R°, aşadar cei trei vectori formează o bază. 
b) Vectorul u se scrie u=au, +au, +ou, care conduce la sistemul 
aq +2a, + o =l 
cu soluția: a, =—2, œ, =—2, œ =7, deci : 


up, = (—2,—2,7)',B= uu, u; 


c) Matricea de trecere de la B la B'= v}, V2, V3 


este : 
—] 1 0 
A=| 1 1 —l 
0 1l 2 
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iar inversa lui A! are forma: 


23/5058 
(A'y!=| 2/5 2/5 —1/5 
1/5-1/5: 215 


Folosind relaţia (2.4.6), obţinem în final: 


uy =(—1,—3,2). 


Observaţie. Dacă vectorii noii baze nu sunt exprimaţi de la început prin combinaţii ale 


vectorilor vechii baze, coordonatele lor (deci elementele matricei de trecere A) trebuie calculate 


separat. În această situaţie este mai simplu să calculăm direct up, folosind definiţia. 


În unele cazuri bazele B şi B’ diferă printr-un singur vector. Formula de calcul a noilor 


coordonate va avea evident o formă mai simplă. Să presupunem că în baza B, vectorul u; s-a 


înlocuit printr-un vector v j care are coordonatele a, Toilette ii adică: 
Vj = Aj + a»ju» +: +aju; + + aaa (2.4.7) 


Ceilalţi vectori ai vechii baze rămânând neschimbaţi, matricea de trecere A devine: 


N O cete OSa O 
0 1l 0 0 
ajj aj aij aj 
0 0 0 1 


După determinarea lui (A!) | şi aplicarea formulei (2.4.6) deducem noile coordonate: 


Q; n 


ij ij 


Formulele (2.4.7), greu de memorat pot fi utilizate prin intermediul unui algoritm de 


calcul bazat pe transformări elementare ale liniilor matricei coordonatelor vectorilor 


u, u, U»,...,U;,...U, şi V j Ilustrăm mai jos modul de organizare a algoritmului. Etapa inițială 


are forma : 
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Baza u; u, u; u, Vj Upg 
u; 1 0 0 0 aij Aj 
u, 0 1 0 0 42; Ap 
u; 0 o0 | 0 a; la 
u, 0 0 0 1 aaj Qi 


Coloanele tabelului conțin coordonatele fiecărui vector în baza de pe prima coloană. Ele 


se presupun cunoscute. Ultima coloană reprezintă coordonatele lui u în baza veche B. 
Etapa finală are forma: 


u, Vj u, Vj Up 
a 
lj 
0 iS 0 0 6i =au aij 
aij ij 
a = aci 
2j B, =a, —a, 
1 H 0 0 la. 
ij 
aij 
1 Q 
1 
0 mii ai 0 | 8; = — 
ij aij 
Q: 
nj n n n 
0 = 1 0 "a; 
a; 


şi conţine coordonatele tuturor vectorilor în noua bază. 


Se observă că linia i din ultima fază se obţine prin împărţirea liniei i din prima fază prin 
a; (care este nenul, altfel matricea de trecere A nu este nesingulară). 
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O linie oarecare k a matricei finale se obţine înmulţind linia i cu —a, după care se adună 
cu linia k a matricei inițiale. Ca efect al acestor operaţiuni, pe ultima coloană se obţin 
coordonatele f;,B.(k zi), pe coloana v; se obţine coloana "i" a matricei unitate I, iar coloana 
lui v; se modifică în mod corespunzător. 

Algoritmul poate fi aplicat înlocuind cauzele cu efectele şi invers. Astfel, putem să 
aplicăm transformările elementare de linii în matricea fazei inițiale pentru a transforma coloana 
lui vj în coloana vectorului u; din faza iniţială. Unii numesc această operaţie “Lema 
substituţiei” şi aplică diverse metode mnemotehnice de utilizare. Noi o considerăm o simplă 


aplicare a transformărilor elementare şi o folosim succesiv înlocuind toţi vectorii bazei iniţiale B 


prin vectorii bazei B’. 


Exemplu. În baza canonică B, = eee se consideră vectorul u = (1,1,—1)'. Să se 
determine coordonatele sale în baza B = yVas6a unde 


v; = (1,2,3), v, =(1,—1,1)*, v, =(0,0,D)'. 


Baza «i èz & Mi V3 u 

e, l 0 0 l l l 

e, 0 1 0 2 -1 1 

e; 0 0 1 3 1 -1 

Mi 1 0 0 1 l l 

e, —2 l 0 0 -3 -1 

& -3 0 1 0 -2 —4 

M 13 13 0 1 0 2/3 
V2 2⁄3 -13 0 0 1 1/3 

& -5/3 —B 1 0 0 —10/3 


Operaţiile — extrem de simple — au urmărit crearea matricei unitate pe coloanele 


vectorilor noii baze B = v,,v»,e 
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În prima etapă am completat tabelul cu coordonatele — coincizând cu componentele — 
tuturor vectorilor în baza canonică B.. S-au aplicat transformări elementare de linii pentru a 
obţine pe coloana lui vı prima coloană a matricei unitate I. În ultima etapă, prin transformări 
elementare de linii se obțin pe coloana lui v> cea de a doua coloană a matricei unitate I. 


Coordonatele lui u în noua bază sunt înscrise pe ultima coloană şi deci: 


2 l 10 
e 


Observaţie. Chiar dacă pare prea complicat în comparaţie cu aplicarea definiţiei în 
determinarea coordonatelor lui u, algoritmul are avantajul principal că în fiecare etapă coloanele 
conţin coordonatele în baza curentă. În plus, ca avantaj secundar este posibilitatea de a-l folosi în 


probleme de programare liniară. El va fi aplicat consecvent în capitolele care urmează. 


34 


US3 Forme pătratice 


Formele pătratice constituie un caz particular de funcţii definite pe un spațiu 
liniar. Ele sunt întâlnite în numeroase capitole de matematici economice (chiar unele modele de 
creștere economică utilizează formele pătratice). 

Dedicăm această unitate de studiu formelor pătratice făcând mai întâi o prezentare 
a noțiunii, după care facem o clasificare (după semn) a acestora. 

Caracterizarea formei se obține cu uşurinţă după aducerea la expresia canonică 
prin una dintre metodele devenite clasice (metoda lui Jacobi şi metoda lui Gauss). Sunt date 
aplicaţii în acest sens, accentuând aplicarea acelorasi transformări elementare de linii în matricea 
asociată formei. 

Studentul trebuie să-şi însușească mai întâi noțiunile fundamentale. Ulterior 
trebuie să înveţe să asocieze unei forme pătratice o matrice simetrică şi reciproc. Ultimul lucru 
este să folosească metoda transformărilor elementare pentru a aplica metoda lui Gauss de 
obținere a expresiei canonice (sau a metodei lui Jacobi pentru a ajunge la acelaşi rezultat). 

Sunt necesare aproximativ şase ore pentru a finaliza studiul şi a fixa, prin 
aplicaţii, cunoştinţele (două ore pentru studiul teoretic, o ora pentru matricea asociată, două ore 
pentru metoda lui Gauss şi o oră pentru metoda lui Jacobi). 


Sursele bibliografice sunt cele de pe poziţiile [3], [4], [5], [6] şi [8]. 


x 


Considerăm în continuare spațiul liniar R” definit peste R precum şi funcția Q: R” x R” — R 


definită prin: 
n n . . T. 
Q(x)= Sax, cu aj=a; (V) ij=l,n. (3.1) 
i= j= 
Această funcție se numeşte pătratică şi generalizează funcțiile de gradul al doilea 
omogene într-o singură variabilă. 
Matricea A € M (R) care are elementele a; ij=l,n se numeşte matricea formei 


pătratice. Prin verificare directă se poate arăta că are loc formula: 


Q(x)=x'Ax. (3.2) 
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Exemplu: Să se scrie matricea corespunzătoare formei pătratice Q:R5—R definită 
prin: 


Q(x) =2x? —x3 +3x2 — XX) +2X2X3. 


Soluție. Folosim (3.1) sau (3.2) elementul a; al matricei A fiind coeficientul lui x;x; în 


expresia acesteia. Cum x;xj =X;X; şi aj =a; în A coeficienții produselor de variabile diferite 


se înjumătăţesc. Aşadar: 


2 —1/2 0 
A =|—1/2 —1 1 
0 1 3 


şi formula (3.2) se poate imediat verifica. 


Observaţie. Matricea A este pătratică şi simetrică conform condiției a; =a, (nu este 


obligatoriu inversabilă). 
Dacă în matricea A numai elementele diagonalei principale pot fi nenule spunem că 


forma pătratică are expresia canonică. Expresia canonică a unei forme pătratice este: 


QA) = Ax? +x? +.. +A, NER. (3.3) 


O formă pătratică Q este nedegenerată dacă expresia canonică are exact n termeni. Când 
numărul lor este mai mic forma se numeşte degenerată. 
Expresia canonică a unei forme pătratice are o importanță deosebită în studiul semnului 


acesteia având în vedere că Q(x) e R. Există forme pătratice care păstrează acelaşi semn pentru 


orice argument x E€ R°. Nu este greu de verificat că: Q(0)=0. Următoarele definiții pot fi 


relevante: 


Definiţie. Forma pătratică Q: R” — R se numeşte pozitiv (negativ) definită dacă: 


Q(x) > 0 (Q(x) <0 (V)xeR", xz0. (3.4) 
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Definiţie. Forma pătratică Q:1R" — R se numeşte pozitiv (negativ) semidefinită dacă: 


Q(x) 20 (Q(x) <0) (V)xeR” (3.5) 


Orice formă pătratică definită este o formă pătratică semidefinită. Reciproca nu este 


adevărată. In cazul în care: 


x, ye R” ai. 0(x)>0 si Q(y)<0 

atunci forma pătratică se numeşte nedefinită ca semn. 
Dacă o formă pătratică este adusă la expresia canonică caracterul de formă pătratică 
pozitiv (negativ) definită (semidefinită) poate fi mult mai simplu apreciat. Astfel, dacă forma 
pătratică Q(x) definită de (3.1) nu este sub forma (3.3), putem obține expresia canonică a 


acesteia: 
QY) = AY MY? te Fan A ER (3.6) 
unde y,,i= ln sunt date de relațiile liniare: 


y Di Fe biz +H bin 
Yı = bă, + bx, + bouXu ; CU b; ER (3.7) 


Yn = DX + bo F ii Dan 


Teoremă: Fie forma pătratică Q(x) definită de (3.1) a cărei expresie canonică 


asociată este dată de (3.6). Dacă: 


a) 4, >0,i= ln atunci Q(x) este pozitiv definită; 


b) 4, 20,i=l,n si (3) 4, =0atunci Q(x) este pozitiv semidefinită; 
c) 4, <0,i= ln atunci Q(x) este negativ definită; 
d) 4, <0,i= Ln si (3) A, = Oatunci Q(x) este negativ semidefinită; 


e) (3)4, >0 si (34; <0 atunci Q(x) este nedefinită ca semn. 


Vom da în continuare două metode de obţinere a expresiei canonice (3.6). 
a) Metoda lui Gauss. 
Se aplică transformări elementare asupra liniilor matricii A corespunzătoare formei 


pătratice, pentru a o aduce la forma (numită superior triunghiulară) de mai jos: 
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ai 12 ain 
0 a% A 
A = 3.8 
0 0 an (SA 
0 0 a! 
Expresia canonică a formei pătratice este: 
1 1 1 
O= a dea (3.9) 
ai 422 Ann 
unde: 
y, = ax, Fax + tanXn 
I~ a22X2 T- F a2nXn (3.10) 
Ya = an Xa 
b) Metoda lui Jacobi 


Fie A,A, A, lanţul minorilor diagonali ai matricii A corespunzătoare formei 
pătratice, daţi de relaţiile: 


A a2 


A Sau » A= az Âp et A (3.11) 


a2 az 
Dacă toți A,70 ; i=1,n atunci expresia canonică este dată de formula: 


Hoa Aia A 
=— + — N m 18 3 3.12 
Q(y) A, yı A, y2 A. Yn ( ) 


Observaţii: 


1) Expresia canonică a unei forme pătratice nu este unică (coeficienţii A, pot avea valori 
diferite). Totuşi în oricare dintre expresiile canonice asociate aceleiaşi forme pătratice numărul 
de coeficienţi 4, pozitivi respectiv negativi este acelaşi; 

2) Metoda lui Jacobi nu poate fi aplicată dacă măcar un minor diagonal al matricei A 
este nul ((2)A,; =0); mai mult prin această metodă nu se pot obţine expresiile variabilelor y, 


definite în (3.7); 
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3) Aplicând metoda lui Jacobi de obţinere a expresiei canonice (3.12), teorema de 


caracterizare a semnului unei forme pătratice conduce la următoarele concluzii: 


a) Dacă A, >0,A, >0,......... „A, >0 atunci Q(x) este pozitiv definită; 
c) Dacă A, <0,A, >0,A, <0,......... atunci Q(x) este negativ definită; 


e) Dacă A, 70, (v)i= ln şi în orice altă combinație de semne decăt cazul a) sau c) 


atunci Q(x) este nedefinită ca semn; 
Evident cazurile b) respectiv d) de semidefinire nu pot fi caracterizate cu metoda lui 


Jacobi 


Exemplu. Să se găsească expresia canonică a formei pătratice: 


Q(x) = 2x? —2x;Xx, +2x2 —2x,X; +3x2. 


Folosim mai întâi metoda lui Gauss. Obţinem succesiv: 
2 —1 0 1 1 -l 1 1 —l 

A=|—1 2 —l|o|—1 2 —1|n|0 3 —2|- 
0 —1 3 0 —1 3 O —1 3 


l 1 —]l 1 1 —l 
0 —l 3 00 7 


Expresia canonică este aşadar: 


| 
90 arta cau 


Prin metoda lui Jacobi obţinem: 


E ni De ali - 40 
A, =2 ^= a A =] 2 -l=3 
= A 
a 3 


Expresia canonică astfel obținută este deci: 


E E 
z)=-z +z +z 
N) a a a TA 


Prin cele două metode am obţinut coeficienți diferiți (deci expresii canonice diferite), dar 


semnele acestora sunt aceleaşi. Forma pătratică este evident pozitiv definită. 
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US4 Programare liniară. Consideraţii teoretice 


Principalele noţiuni asociate programării liniare sunt pe larg analizate în această 
unitate de studiu. După enumerarea (mai mult deducerea) câtorva modele economice remarcabile 
care conduc la probleme de programare liniară, sunt trecute în revistă formele pe care le poate 
lua o astfel de problemă. Se acordă atenție și soluțiilor problemelor de programare liniară şi 
proprietăților acestora. O atenţie deosebită s-a acordat mulțimilor convexe care constituie o 
noțiune fundamentală. 

Principala preocupare a studentului constă în însuşirea corectă şi completă a 
noțiunilor teoretice. Puţine aplicaţii concrete sunt accesibile, unele dintre acestea fiind analizate 
în cadrul activităţilor tutoriale. Una dintre acestea (programarea în două variabile) care se rezolvă 
prin metoda grafică este deosebit de importantă deoarece poate servi la exemplificare 
chestiunilor teoretice introduse. 

Timpul de lucru este apreciat la şase ore: două ore pentru problema de programare 
în sine (cu forme diverse de prezentare), două ore pentru prezentarea noțiunii de mulțime 
convexă și a rezultatelor teoretice asociate și două ore pentru cele trei tipuri de soluţii cu 
proprietățile lor. 

Orice manual de programare liniară conţine chestiunile enumerate aici în diverse 
maniere. Poziţiile [3] — [9] din lista bibliografică folosesc chiar maniera de prezentare adoptată în 


lucrarea de faţă. 


40 


4.1. MODELE ECONOMICE CE CONDUC LA 
O PROBLEMĂ DE PROGRAMARE LINIARĂ 


a) Consumul lunar de materii prime 


Considerăm un proces de producţie în cursul căruia se realizeză m tipuri de produse 
pi = l,m, pentru care se utilizează n tipuri de materii prime M p= In. Dintr-o unitate de 
materii prime M , se pot produce a, unităţi de produs P,. Se doreşte ca numărul de produse de 
tipul P, (realizate zilnic/săptămânal/lunar /etc.) să fie cel puţin în cantităţile b;. Costul unei 
unităţi de materie primă M, este c. Se cere să determinăm consumul lunar de materii prime, 


încât să se realizeze producţia minimă planificată, iar cheltuielile legate de materiile prime să fie 
minime. 
In continuare vom elabora modelul matematic al acestui fenomen economic. Notăm cu 


x; j=l,n cantităţile de materie primă M, ce urmează a fi utilizate. Vom scrie sub formă 


algebrică restricţiile economice impuse mai sus. Aceste restricţii sunt de trei categorii: restricții 
legate de planificarea producţiei, restricţii impuse de sensul economic al variabilelor şi restricţia 


de minimizare a cheltuielilor. Aşa cum a;x, reprezintă cantitatea de produse de tipul P, ce se 


obține dintr-o cantitate x, de materie primă M;, deducem că: 
n . 
aX Apă F... F An Xn = Da, ; i=l,n 


reprezintă cantitatea totală de produse de tipul P,, ce se obține când materiile prime M, se 


utilizează în cantitățile x,,j=1,n . Dar această cantitate nu poate fi mai mică decât b, şi deci 


obținem inecuațiile liniare de forma: 
Saab, >b, i=l,m (4.1) 
j=l 


Cum x, reprezintă mărimi economice (cantități de materie primă) rezultă obligatoriu 


impunerea condițiilor (de nenegativitate): 


x20, j=l, (4.2) 
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Cheltuielile totale obţinute ca sumă a cheltuielilor unitare vor fi date de funcția: 
n T h 
E Do ii X= XXX}, ER, (4.3) 
J= 


care se cere a fi minimizată. 

Din punct de vedere matematic problema poate fi formulată astfel: să se găsească în 
mulţimea soluţiilor sistemului de inecuaţii liniare (4.1) care satisfac condiţiile de nenegativitate 
(4.2), o soluţie pentru care funcţia f ia valoarea minimă. Problema are sens dacă sistemul (4.1) 
are cel puţin o soluţie. Inecuaţiile (4.1) le numim restricțiile problemei, (4.2) condiţiile de 


nenegativitate, iar f poartă denumirea de funcţie obiectiv. 


b) Folosirea optimă a resurselor 


Într-un proces de producţie în care se pot realiza n tipuri de articole Aj, j=l,n sunt 


utilizate m resurse R;, 1>1,m, care sunt limitate de cantitățile b;. Pentru producerea unei unități 


de articol de tipul A; se consumă o cantitate a, de resursă de tipul R;. Profiturile unitare pentru 


articolele A; presupunem că sunt c;, j=l,n. Să se determine cantitatea x, de articol A, j=l,n 


care urmează a se realiza în procesul de producție considerat, cu resurse limitate, astfel încât 
procesul să aibă drept rezultat profit maxim. 


Limitările resurselor conduc la inecuațiile: 


Sax, <b, i=l, (4.4) 
j=l 


Necunoscutele (cantitățile) x,;, j = 1,1, trebuie să satisfacă condiţiile de nenegativitate: 
x; >0, j=l,n (4.5) 
Profitul total este dat de funcţia liniară 
n 
F(x) = 3 cx; (4.6) 
j=l 
care trebuie să ia valoarea maximă. 


sa uite y ; T n sala 
Ne propunem să găsim (măcar) o soluţie xo=— x, X2,- X, ER (numită în cele ce 


n 


urmează soluţie optimă) a sistemului de restricții economice (4.4) care să aibă componentele 
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nenegative, deci să verifice condiţiile (4.5) şi în care funcția obiectiv (profit) să ia cea mai mare 
valoare posibilă. Evident, în modelul matematic (corect) al unei probleme economice reale, 
sistemul de inecuaţii liniare (4.4) are o infinitate de soluţii care verifică condiţiile de 


nenegativitate (4.5). Dorim să găsim acea soluţie xo pentru care are loc: 


(min) f(x)= f (x0) 


b) Problema amestecului optim 


Să presupunem că pentru realizarea unui amestec se pot întrebuința n materii prime 
M, j= Ln, disponibile în candități nelimitate, ce conțin m substanțe S,, i = l,m. O unitate de 
materie primă M;, conține cantitatea a; din substanță S;. Amestecul trebuie să conțină 
substanța S; în exact cantitatea b;. Costurile unitare ale materiilor prime M, sunt c; unități 
băneşti. Să de determine cantitățile x;, j =],n de materie primă M, necesare realizării 


amestecului cu compoziția prescrisă şi la un cost total minim. 
Restricțiile sunt şi acum de trei feluri: restricțiile de compoziție a amestecului, restricțiile 
datorate sensului economic al variabilelor şi restricția de minimizare a costului total. 
Traducerea algebrică a primelor două categorii este: 
Sax: =b, i=l,m (4.7) 
x; >0, j=l,n (4.8) 


Costul total al materiilor prime necesare pentru realizarea unei unități de amestec este dat 


de funcția: 


f(x)= 3 ex; (4.9) 
= 


care se cere a fi minimizată. 
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4.2. FORMULAREA GENERALĂ A UNEI PROBLEME 
DE PROGRAMARE LINIARĂ 


Se constată că modelele matematice particulare, ataşate problemelor economice 
prezentate anterior prezintă o serie de asemănări, fapt ce permite înglobarea lor într-un model 
general. Să considerăm un sistem care poate conţine ecuaţii, inecuaţii de un semn sau altul, adică 


un sistem de forma: 


At FAX Fm Sb 
Aa +X FAX teamă Sb (4.10) 


d +X Fan% F- -a mnXn s la 


numit sistemul de restricții (economice) al problemei, care poate fi scris condensat (utilizând 


indicele de sumare) şi sub forma: 
Zax; Sb, i=l, (4.11) 


(semnul $ înseamnă semnul =, ori <, ori >, ori <, ori >). Condițiile de tipul: 


x; >0, j=l,n (4.12) 
sunt numite condiții de nenegativitate şi sunt necesare datorită caracterului economic real al 


variabilelor (d.p.d.v. strict matematic putem rezolva şi problemele în care aceste condiții nu 


apar). Funcția liniară: 
f(x) = > c;x; (4.13) 
j=l 


care trebuie să o optimizăm (minimizăm sau maximizăm) se numeşte funcție obiectiv (funcție 
cost sau funcție profit conform tipului de problemă studiat). 

Problema generală a programării liniare se formulează astfel: în mulțimea soluțiilor 
sistemului (4.11) ce satisfac condițiile (4.12) să se găsească o soluție (numită soluție optimă) ce 


dă valoarea minimă sau maximă funcției definite de (4.13). 
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4.3. FORME ALE UNEI PROBLEME 
DE PROGRAMARE LINIARĂ 


Modelul matematic al unei probleme de programare liniară poate fi scris sub diverse 
forme potrivit scopurilor particulare urmărite. Aceste forme au doar o semnificaţie strict 
matematică şi uşurează respectiv simplifică demonstrațiile anumitor rezultate matematice 
(propoziţii sau teoreme) absolut necesare în fundamentarea modului (algoritmului) de rezolvare a 


problemelor de programare liniară. 


a) Forma matriceală 


. A . A i a, T 
Considerând matricele: A = (a; ),_ E M mns X5 (X ase) ER”, 


i=l,m 


j=l,n 
b= (b, b»... b,) ER” şi C=(c:6,....,C,)E R”, putem scrie problema de programare liniară 
sub forma: 


(min/max) f(x) = C- x 


Ax Sb i (4.14) 
x>0 


b) Forma vectorială 


Dacă introducem vectorii: 


T T 


T 
A = djasa oB = osoase EI E // AL, S ARD 


Po=b. 


m m2 


atunci problema poate fi scrisă sub forma: 


(min/ max) f(x) = C-x 
XA +h, +.. tă Sh (4.15) 
x>0 


c) forma canonică 


O problemă de programare liniară este sub formă canonică dacă ea se scrie: 
(min) f(x) = Cx 
Ax >b 


2 (4.16) 
x>0 
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sau 
(max) f(x) = Cx 
Ax <b (4.17) 
x>0 
d) forma standard 
O problemă de programare liniară este sub formă standard dacă restricțiile ei sunt sub 
formă de egalități. O problemă de programare liniară poate fi scrisă sub forma standard 


matricială sau vectorială, adică: 


(min/max) f(x)=cx (min/ max) f(x) = C-x 
Ax =b sau A + x a e = B (4.18) 
x >0 x >0 


Observaţie. Formele (4.16), (4.17), (4.18) ale unei probleme de programare liniară sunt 


echivalente dacă ţinem seama de următoarele: 


1. Sensul unei inegalităţi se schimbă prin înmulţirea cu —1; 


2. O inecuaţie de forma ax +a +... +aX, Sb;, poate fi scrisă ca o ecuaţie 


in n — 12 


č ENS ERTE oa eit 
AX + Apo +... Fă, +X; =b;, prin adăugarea a unei noi variabile (variabilă de 
compensare) x; > 0; 


3. O inecuație de forma: ax, +a +.--+4;,X, >b;, poate fi scrisă ca o ecuație 


in n — 1? 


a X FAX +... Fă, —X; =b, , scăzând o variabilă de compensare x; > 0 ; 

4. O ecuație poate fi scrisă ca două inecuații de semn contrar ; 

5. Problemele de maxim pentru funcția f(x) pot fi reduse la o problemă de minim 
pentru funcția f (x) =— f(x); cu alte cuvinte în cazul unei functii (forme) liniare de 
tipul (4.13) are loc egalitatea : 

(max) fœ) = —(min)(—f (00) = Amin) f (x) , pentru (V)x eR" 
Având în vedere cele de mai sus, în continuare ne vom ocupa numai de probleme aflate 


sub forma standard cu cerința de minim pentru funcția obiectiv. 
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4.4. MULȚIMI CONVEXE 


Fie spaţiul vectorial R” şi M o submulțime nevidăasa(MCR", M = Ø). 


Definiţia 4.1. Mulțimea M CR" se numeşte convexă dacă (1)x, x» EM şi (We 0,1 


avem. 


x = x, +(1—1)x, EM (4.19) 


Definiţia 4.2. Un punct (vector) xEM se numeşte punct extrem sau vârf al mulţimii 
convexe M, dacă nu există x, z X, din M şi nu există A € (0,1) astfel încât 
x = Ax, +(1—A)x, (4.20) 


în caz contrar, punctul se numeşte punct interior. 


Definiţia 4.3. Numim combinaţie liniară convexă a vectorilor x; EM, 1>1,k o expresie 


de forma: 


k k 
AX, cu À El0,l] XA =1 (4.21) 


i=l i=1 


Observație. Se poate demonstra (relativ uşor) că orice punct interior se poate scrie 


(exprima) întotdeauna ca o combinaţie liniară convexă de punctele extreme ale mulțimii 


convexe; evident, conform definiției 4.2, un vârf (punct extrem) al unei mulțimi convexe nu 


poate fi scris ca o combinație liniar convexă de puncte interioare mulțimii. 


4.5. SOLUȚII ALE UNEI PROBLEME 
DE PROGRAMARE LINIARĂ 


Fie o problemă de programare liniară aflată sub forma standard matricială: 
(min) f(x)=cx 


Ax =b (4.22) 
x>0 
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Definiţia 4.4. Un vector x € R" ale cărui componente x p» J= ln verifică sistemul de 


restricții şi condițiile de nenegativitate se numeşte soluția admisibilă a unei probleme de 


programare liniară. 


Mulțimea soluțiilor admisibile o vom nota cu S. 
Teorema 4.1. Mulțimea soluțiilor admisibile S, este o mulțime convexă. 


Demonstraţie. Fie x, =x, ES, două soluții admisibile. 
Atunci avem: 
Ax, =b, x >0 ; Ax, =b, x, >0. Să considerăm un vector x de forma: 
x = Ax, +(1—À)x,, cu à € 0,1. 
Să arătăm că x €S, . Într-adevăr vom putea scrie: 
Ax = A 13, +(1—A)x, =1Ax H1 —A)Ax, = àb+(1 —À)b = b 
şi deci verifică sistemul de restricții. 
Cum A€ 0,1, iar X;, X, >0, rezultă Ax, +(1—A)x, =x>0 şi deci x este soluție 


admisibilă a problemei de programare liniară. 


Definiția 4.5. O soluție admisibilă care are cel mult m componente strict pozitive, iar 
vectorii formaţi cu coloanele matricei A, corespunzătoare acestor componente sunt liniar 


independenți, se numeşte soluție admisibilă de bază. 


Dacă soluția admisibilă de bază are exact m componente strict pozitive, spunem că este 
nedegenerată, iar dacă are mai puține, degenerată. Mulțimea soluțiilor admisibile de bază o vom 


nota cu S şi este o mulțime finită (card $,, < C” ) deci evident, nu este o mulțime convexă. 


Exemplul 4.1. 
Fie problema de programare: 


(min) f(x) = 2x, +x, + 2x4 


x + 2X +x = 2 
2X% =X, +2x +x, =4 
x; >0,j =1,4 
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În baza definiţiilor (4.4) şi (4.5) constatăm că: 


x, = 2/3, 1/3, 2/3, 5/3 1 este soluţie admisibilă 
x= 0, 1,0, 5 T este soluţie admisibilă de bază nedegenerată 
x= 0, 0, 2, 0 T este soluţie admisibilă de bază degenerată 


T i în atit i 
x4= 1, 0, 1, 0 este soluție admisibilă, dar nu este de bază deoarece 


(vectorii Pı şi P3 nu sunt liniar independenți) 


Teorema 4.2. Orice soluție admisibilă de bază a unei probleme de programare 


liniară, este un punct extrem al mulțimii soluțiilor admisibile S, şi reciproc, orice punct 


extrem al mulțimii soluțiilor admisibile este o soluție de bază a problemei de programare 


liniară, cu alte cuvinte are loc: 


xESp p © xes, este punct extrem (4.23) 


Definiția 4.6. O soluție admisibilă x se numeşte soluție optimă, dacă realizează minimul 


sau maximul funcției obiectiv f, mulțimea soluțiilor optime o vom nota cu So 


Se poate demonstra că mulţimea soluţiilor optime este o mulţime convexă. Evident, 


teoretic Sọ poate să fie vidă (problema de programare aferentă nu are soluţii optime), să aibă o 


infinitate de elemente (în acest caz spunem că problema de programare aferentă are o infinitate 
de soluţii optime). Dar cazul cel mai des întâlnit şi cel care apare în imensa majoritate a 
problemelor economice reale este cel în care S, are numai un element (problema de programare 


aferentă are soluţie optimă unică). 


Teorema 4.3. Dacă o problemă de programare are o soluţie optimă finită, atunci 


există un punct extrem al mulţimii soluţiilor admisibile S, în care funcţia obiectiv ia 


aceeaşi valoare, adică: 


(min) f(x) = fo e R(finiD > (Axe Sai. f= fo (4.24) 
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USS Metoda Simplex 


În unitatea de studiu precedentă s-a prezentat modelul de bază al unei probleme 
de programare liniară. În unitatea prezentă dăm o metodă de rezolvare a acesteia: metoda 
Simplex. Aici se vede importanţa unităţilor de studiu anterioare (USI şi US2) deoarece se 
lucrează cu vectori (elemente ale spațiului vectorial al restricțiilor) şi folosim (pentru 
determinarea componentelor soluţiei) metoda transformărilor elementare. De mare utilitate este 
algoritmul prezentat în US2 care se aplică aici fără modificare. 

Se dau principiile algoritmului Simplex (criteriile de optim, de intrare, de iesire 
etc.) şi se descriu etapele algoritmului. Se prezintă apoi cazul în care soluţia inițială trebuie 
construită (metoda celor doua faze). 

Studenţii trebuie să rețină mai întâi criteriile, apoi au de reţinut forma tabelului pe 
care se aplică algoritmul Simplex, cu precizarea clară a modului de completare a fiecărui element 
al acestuia. Ulterior ei trebuie să deprindă ordinea în care se aplică cele trei criterii fundamentale 
şi calculele care se fac pentru fiecare. Ultimul lucru îl constituie efectuarea calculeleor 
(transformările elementare) pentru obținerea unei soluții mai bune. După fixarea temeinică a 
algoritmului se poate trece la metoda celor două faze. 

Sunt necesare zece ore pentru parcurgerea acestui material (două ore pentru 
parcurgerea principiilor, o oră pentru descrierea algoritmului, patru ore pentru aplicații ale 
algoritmului, trei ore pentru metoda celor două faze). 

Din lista bibliografică sunt utile poziţiile [3] — [6], [8], [9]. În [7] este prezentată o 


variantă puţin modificată dar interesantă. 
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5.1. PRINCIPIILE METODEI SIMPLEX 


Să considerăm problema de programare liniară sub forma standard matricială: 


(min) f(x)=cx (*) 


Ax=b (%*) (5.1) 
sau vectorială 
(min) f(x) = C-x (%) 
yA +ah +.. thh E) (5.2) 


Raţionând în mod simplist, pentru a rezolva o problemă de programare liniară (adusă la 
forma standard), ar trebui să parcurgem următorii paşi: 
1. Să determinăm toate soluţiile sistemului liniar compatibil nedeterminat (**), lucru relativ 
uşor şi cunoscut încă din liceu; 


2. Să eliminăm (cum!?) din infinitatea de soluţii a lui (**) acele soluții (tot o infinitate) care 


nu verifică condiţiile de nenegativitate (***); 
3. Din infinitatea de soluţii admisibile rămase (adică soluţiile sistemului (**) care satisfac 


condiţiile (***)) să o găsim (cum!?) pe aceea sau acelea care fac funcţia obiectiv să ia 


cea mai mică valoare, adică care satisfac(e) şi condiţia (*). 

Rezultatele indicate în capitolul precedent ne oferă un mod simplu şi elegant de 
determinare a soluţiei optime a unei probleme de programare liniară. Astfel, teoremele 3.1.2 şi 
3.1.3 ne indică că trebuie să căutăm soluţiile optime ale unei probleme de programare liniară, în 
mulţimea soluţiilor admisibile de bază (S,, este o mulţime finită având cel mult CF elemente). 
Cu alte cuvinte, va trebui să parcurgem următorii paşi: 


1. Determinăm (eventual folosind transformări elementare) toate soluţiile de bază ale 


sistemului (**); sunt cel mult C” soluţii de bază, deci un număr finit (şi nu infinit); 
2. Eliminând soluţiile cu componente negative (neadmisibile), care nu satisfac condiţiile 
(***), obținem mulțimea S „a soluțiilor de bază admisibile, printre care se află şi soluţi- 


a/ile optim-ă/e căutat-ă/e; 
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3. Calculăm valoarea funcției obiectiv în toate soluțiile admisibile de bază aflate la pct. 2). 
Soluția optimă va fi cea în care funcția obiectiv f(x) ia valoarea minimă. 

Totuşi, chiar şi în cazul problemelor în care numărul restricţiilor (ecuaţiilor)m respectiv al 
variabilelor (necunoscutelor) n este relativ mic, numărul de calcule necesar aplicării metodei de 
mai sus este extrem de mare, ceea ce o face aproape imposibil de aplicat în realitate. De 
exemplu, pentru m=10, n=40 numărul de soluții de bază al sistemului (**%), care trebuie 
determinate, este de Cp = 846.660.528 soluţii (!!). 

Vom prezenta în continuare o metodă (algoritm) numită metoda (algoritmul) SIMPLEX, 
elaborată de matematicianul german G.B. Dantzig în anul 1906, algoritm care are avantajul unui 
volum extrem de redus de calcule şi utilizarea numai a transformărilor elementare ca aparat 
matematic. 

Să presupunem fără a restrânge generalitatea că matricea lărgită a sistemului de restricţii 


este de forma: 


Î 20-a URC aa. E 
O zare: 10 la e aa ROS 
2m+1 2n e (5.3) 
OA E E EEE E E 
Vectorii Pı, P2, . . ., Pm formați cu primele m coloane ale matricei A, fiind liniar 


independenți, formează o bază în R™. În acest caz un vector oarecare P;, j=m+1,n se exprimă 


în mod unic cu ajutorul vectorilor bazei sub forma: 
m 
= 

Să notăm cu Cp vectorul linie format cu primele m componente ale vectorului C, adică 


coeficienţii variabilelor x,,x,,...,x, din funcţia obiectiv ce corespund vectorilor care au intrat în 


bază. 


Introducem mărimile: 


z-c;=z;-CBP, = d j=l,n (5.5) 


ij > 


Dacă presupunem acum că f;, i=1,m sunt pozitive şi dacă avem în vedere (5.3) deducem 


că vectorul x de componente, 
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X, = BX; > Pose -Xn Bms Xm > 0---» Xp>0, este o soluție admisibilă de bază a 


problemei de programare considerată. 


Teorema 5.1. (Criteriul de optim) Dacă diferenţele z, —c; S 0, j= m+l,n, atunci 
problema de programare liniară are optim finit şi x este soluţia optimă căutată. 

Observaţii: 

a) se constată uşor că z; -cj =0 pentru j= l,m 

b) În problemele de programare, în care avem cerința de maxim pentru funcția obiectiv, 
criteriul de optimalitate va fi: z jc; 20, j=m+l1,n 

c) Din demonstrația teoremei rezultă că în cazul în care există z; — c; >0 în probleme de 


minim pentru funcția obiectiv pot exista soluții admisibile pentru care funcția obiectiv să ia o 
valoare mai mică decât pentru X, deci X nu este optimă. 


In cazul în care X nu este soluție optimă vom căta să obținem din vechea soluție de bază 
admisibilă X o nouă soluție de bază admisibilă y “mai bună” decât vechea soluție X, în sensul 
că valoarea funcției obiectiv f(x)este mai “bună” (adică mai mică) decât în vechea soluție ( 

fO)< fŒ). 
Dantzig demonstrează că, acest lucru se poate obține înlocuind un vector P; (corespunzător 


unei variabile principale a soluției X )din bază, cu un vector Pj (corespunzător unei variabile 


secundare a soluţiei X) din afara bazei. Să notăm cu yk , k=1,2,...,1—1,Jj, i +1,...,m 


componenetele noii soluţii de bază y. 


Din (1.2.5) înlocuind b, cu By, â; cu a Xp CU Yg, Xj cu B; obţinem: 


ij? j 


B i (5.6) 
y= , ksi 
(94 
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Vectorul P; (care iese din bază) şi vectorul Pj (care va intra în bază) se determină pe baza 
următoarelor două teoreme: 
Teorema 5.2. (Criteriul de ieşire din bază). Condiţiile necesare pentru înlocuirea 


vectorului P; cu vectorul P; sunt: 


1 


asf snina =f, oii ay >0, (5.7) 
Q; 


Să vedem acum ce condiții trebuie să îndeplinească un vector P; pentru a intra în bază în 
locul vectorului P; pentru a obține o soluție "mai bună". 

Teorema 5.3. (Criteriul de intrare în bază) Condiția pe care trebuie să o 
îndeplinească un vector Pj, j=m+1,n pentru a intra în bază în locul vectorului P; în 


vederea obținerii unei soluții "mai bune" este: 


zj—cj>0 (5.8) 


5.2. ETAPELE ALGORITMULUI SIMPLEX PRIMAL 


Rezultatele enunțate în capitolul anterior se concretizează în următorul mod de lucru, 
cunoscut în literatura de specialitate sub denumirea de algoritmul Simplex primal. Acesta 
constă dintr-o succesiune de iterații, care îmbunătățesc treptat o soluție admisibilă de bază 
inițială. Pe baza considerațiilor din aliniatul 5.1. putem pune în evidenţă etapele algoritmului care 
ne vor conduce la soluția optimă a problemei dacă aceasta există. 

1. Se determină o soluție admisibilă de bază, fie aceasta X (oare e atât de simplu!?); 


2. Se construieşte tabelul Simplex corespunzător acestei soluții (este redat mai jos); 
3. Se calculează diferențele zj —c;, j=1,n aflate pe ultima linie a tabelului inițial. 
Aplicăm criteriul de optim. Putem avea următoarele patru situații: 
a) dacă z,-c,<0, j= m+1,n atunci soluția admisibilă de bază X este optimă (şi 
unică!); 
b) dacă z, =c, <0, j= m+l,n şi (3)z ;=€; = iar vectorul P, corespunzător are şi 


componente strict pozitive, atunci soluția admisibilă de bază X este optimă (dar 
nu este unică; problema admite o infinitate de soluţii optime cu aceeaşi valoare 


a funcţiei obiectiv); 
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c) Dacă există indici j, j=m+l1,n pentru care z; —c; 20, iar toate componentele 


vectorului Pj corespunzător sunt mai mici sau egale cu zero (a; <0), atunci 
problema de programare are optim infinit; 

d) Dacă (3)J;j= m+1,n pentru care z,—c; >0 şi nu toţi a; <0, atunci se trece la 
etapa următoare (soluţia nu este optimă). 

Se aplică criteriul de intrare, determinând vectorul P; care intră în bază ca fiind acel 

vector pentru care: 


Z; —C: = max (Z, —c 
ic rax y k) 


Se aplică criteriul de ieşire. Părăseşte baza vectorul P; pentru care 
^i — min Bă Să <> 0, = min!6, >0) 


Se face schimbarea de bază obţinând o nouă soluţie admisibilă de bază y , pentru care 


etapele algoritmului se reiau. Datele problemei de programare liniară şi elementele ce 
au intervenit în teoremele 5.1, 5.2, 5.3 respectiv în algoritmul prezentat mai sus, le 


sintetizăm în următorul tabel Simplex: 


C Co .. Goo Cm .. G Ch 

B Cp Po p 9 _fo 
Po P; Pr ip Pà m 

1 ` 

P, c |6 1O00 e ICE IE A a An 0, 

P» c2 B» 0 1 w e 0) e Oj ze An 6, 

P; ci | 6&6 10 0 Lo 0 a; a, |? 

Pm cm Ba 10 0 0O > 1 CR Am | On 

Z;-C; z 0 0 0 0 Z;-C; Zn “Cn 
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Observaţii: 


a) Valoarea funcției obiectiv pentru soluția admisibilă de bază este dată de: 
f(x) = cip; = CP 
i=l 
b) Diferenţele z; —c, se calculează cu relațiile: 


i) Dacă soluția nu este optimă (există z; —c, > 0) alegem diferența pozitivă care este cea 


mai mare, determinând în felul acesta vectorul P; care intră în bază. 

ii) Dacă există mai multe astfel de diferențe, atunci poate intra în bază oricare dintre 
vectorii ce corespund acestora. 

iii) Dacă pe coloana vectorului Pj care trebuie să intre în bază toate elementele sunt mai 
mici sau egale cu zero, atunci problema admite optim infinit. 

iv) Dacă există şi elemente strict mai mari ca zero considerăm toate rapoartele dintre 
componentele lui Po şi componentele strict pozitive ale vectorului P;, rapoarte de forma 


9 = b; _ componentele lui P) 


L 


>0 cu 4;>0. Cel mai mic raport 6, >0 va indica vectorul care 


4; componentele lui P, 


părăseşte baza. Dacă vectorul P; părăseşte baza atunci noua soluţie admisibilă de bază se obţine 


cu ajutorul unui pivotaj cu elementul pivot oj. 


Exemplu 5.1. Pentru fabricarea a 3 tipuri de produse P,, j=1,3 se utilizeată 2 tipuri de 


resurse Rı şi R2. Consumurile specifice, cantitățile de resurse şi profiturile unitare sunt date în 


tabelul: 


P Limitări 
P; P2 P; resurse 
Rı l l 2 30 
R2 2 3 l 57 
Profit 4 l 5 
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Să se determine cantităţile de produse P), j=1,3 ce urmează a fi realizate pentru a obține 


un profit maxim. 
Soluție. Notăm cu x,x,,x, cantitățile de produse de tipul Pı, P2, P3 ce urmează a fi 
realizate. Modelul matematic al problemei este: 
(max) f(x) = 4x, +x, +5x% 
X + x, +2x <30 


2x +3x, + e <57 
x; 20, j=1,3 


Aducem problema la forma standard folosind variabile de compensare x, şi x; şi 
considerăm problema de minim pentru funcția —f . 
(min) (— f&)) = —4x, — x, —5% 
x+ x% kr 2 +x; = 30 


2x +3x, + x% +x; =57 
x; >0,j=1,3, x4 20, x; >0 


Vom avea problema: 


O soluție admisibilă de bază este x= 0, 0, 0, 30,57 E 


Tabelul simplex se prezintă astfel: 


Zi] -1171 0 E: ga | i 


În prima etapă au intrat în bază vectorii x4 şi x5. Pe linia diferențelor z; —c, observăm 


că există diferenţe strict pozitive, deci conform cazului 3) de la criteriul de optim soluția iniţială 
nu este optimă. 


Cea mai mare diferență pozitivă este z-c> =5, şi conform criteriului de intrare, 
vectorul P; va intra în bază. 


Avem rapoartele 6, = Lo PSR: 15, 0, = Ths 2 = 57. Cel mai mic raport este 
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6, =15, ceea ce înseamnă (conform criteriului de ieşire)că vectorul P? părăseşte baza fiind 
înlocuit de vectorul P3 . Un pivotaj cu element pivot 2 (încercuit) conduce la o nouă soluție 
admisibilă de bază x= 0,0, 15,0, 42 T care nu este optimă (există diferenţe z ;,—c,>0 în 


tabelul corespunzător noii soluții găsite). 


Continuând cu o nouă iterație cu elementul pivot 3/2 (încercuit), obţinem soluția 
admisibilă de bază x= 28,0, 1,0,0 $ care este optimă şi unică. 

Lăsând la o parte variabilele de compensare introduse, obţinem soluţia optimă a 
problemei iniţiale: x = 28, 0,1". 


Din punct de vedere economic aceasta înseamnă realizarea a 28 unităţi de produs P; şi a 


unei unităţi de produs P; , profitul maxim posibil de realizat fiind de 117 (u.m.). 
Exemplul 5.2 Să se determine soluţiile optime ale problemei: 


(min) f(x) = —3x, — 2x, — 4x; 
x +2x, — x +x, =7, 

x — 4x +x +x =l 

x; > 0; j=1,5 


O soluţie admisibilă de bază iniţială este x = 0, 0,0,7,1 T Tabelul Simplex corespunzător este: 
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a ci las dr pata sd mie 0 
Pi Po RoR P P: 0, 
P, 0 | 7 | De ep 3] o | 8%<0 
P, 0 | 1 -4 © 0 1 | 8= 
ze e 100 [8 2 4 0o 0 
P, 0 8 Do e j 1 | <0 
P |-4] ı Îi A d 10 1 | 8%<0 
Zid | ei Vea  08 0 e red 


; : 7 : Aa Sa Ea la T 
Algoritmul se opreşte, chiar dacă soluţia admisibilă de bază găsită x = 0,0,1,8,0 nu 
este optimă (diferenţa z, —c, =18>0), deoarece vectorul P» care ar trebui să intre în bază are 


toate componentele negative. Problema de programare liniară considerată admite optim infinit 


conform cazului 3).. 
5.3. METODA CELOR DOUĂ FAZE 


Algoritmul Simplex este un algoritm rapid convergent la soluția optimă. Astfel, plecând 
de la o soluţie de bază admisibilă inițială a sistemului în formă standard, cu m ecuaţii şi n 
necunoscute, după cel mult m+n iterații obţinem soluţia optimă. De exemplu, în cazul unui 
sistem cu 10 ecuaţii şi 40 de necunoscute, plecând de la o soluție admisibilă iniţială, în cel mult 
10+40=50 de paşi (iterații) se ajunge la soluția optimă (dacă există). 

Determinarea (eventual prin transformări elementare) a unei soluţii de bază pentru un 
sistem liniar, chiar de dimensiuni mari, este extrem de simplă. Problema care apare, este 


admisibilitatea acesteia. De exemplu, în cazul sistemului prezentat anterior cu Cipo = 846.660.528 


soluţii de bază, se poate întâmpla ca foarte multe dintre acestea, sau chiar toate, să fie neadmisibile. 
Am putea astfel calcula un număr imens ( sute de milioane) de soluţii neadmisibile (deci inutile), 
până obținem acea soluţie admisibilă iniţială. După găsirea acesteia, aplicând algoritmul Simplex 


vom obține soluţia optimă foarte rapid, în cel mult 50 de iterații. 
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Vom indica o metodă mai lesne de aplicat în practică pentru a obţine o soluţie admisibilă de 
bază inițială: metoda celor două faze. 

Să considerăm că problema de programare liniară (5.1) este sub formă standard, termenii 
liberi sunt pozitivi şi că matricea A a sistemului de restricţii nu conţine nici o coloană a matricei 


unitate. Vom adăuga la fiecare restricţie o variabilă nenegativă, numită variabilă artificială. Fiind 


m restricții vom adăuga m variabile artificiale xp, k=n+1, n+m , şi considerăm o nouă problemă 


de programare liniară, numită problema artificială ataşată problemei iniţiale, de forma: 


A zi n+m 2 
(min) g(x )= 2 x 
k=n+1 
Ax +I x?= b (5.9) 
x > 0,x" >0 
Aici I i itate de ordinul LXS ee n i | variabilel 
ICi 1, este matricea unitate de ordinul m şi X = Xpy Xp Xu+m vectorul variabilelor 


artificiale. Se constată că orice soluție admisibilă a problemei (5.9) în care variabilele artificiale 
sunt nule, devine după înlăturarea acestora o soluţie admisibilă de bază a problemei (5.1). 


Metoda celor două faze constă în: 


Faza întâi. Se determină soluţia optimă a problemei (5.9) folosind algoritmul Simplex. 
Soluţia de bază admisibilă iniţială artificială, se determină direct din sistemul artificial, considerând 


variabile principale pe cele artificiale. 


Aplicând problemei (5.19) algoritmul Simplex primal, ajungem la una din situaţiile: 
a) (min) g(x*)=0 şi nici un vector corespunzător variabilelor artificiale nu este în bază. 
b) (min) g(x*)=0 şi în bază rămân şi vectori corespunzători variabilelor artificiale. 
c) (min) g(x")z0. 
d) 
Faza a doua, În situaţia a) soluţia optimă a problemei (5.9), după înlăturarea variabilelor 
artificiale, reprezintă o soluţie de bază admisibilă, inițială, nedegenerată pentru problema (5.1). Cu 
această soluţie inițială, se aplică algoritmul Simplex primal în vederea găsirii soluţiei optime a 


problemei inițiale. Din tabelul final al fazei întâi se elimină elementele ce țin de variabilele 
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artificiale (coloanele vectorilor artficiali), se modifică coeficienţii funcţiei obiectiv şi se calculează 


noile diferenţe z, —c,, obținându-se primul tabel Simplex al fazei a doua. . 


Situaţia b) conduce la o soluţie admisibilă de bază degenerată pentru problema inițială. 
Consideraţiile anterioare rămân valabile, însă nu vor fi eliminate coloanele corespunzătoare 
variabilelor artificiale, pentru care vectorii artificiali corespunzători au rămas în bază. 


În situaţia c) problema iniţială nu admite soluţii admisibile. 
Exemplul 5.3. Să se determine soluţiile optime ale problemei de programare liniară: 


x +3x +2% >3, 
3x +3 + X 24 
x; > 0, j=1,3 
(min) f(x) = 3x, — 4x, +2x% 


Aducem problema la forma standard, ea devenind: 


x +3% +2% — Xa =3, 
3x +3x, + x — x; =4, 

x; >0, j=1,3, x4 > 0, x5 >0 
(min) f(x) = 3x — 4x, +2x 


Vom întrebuința metoda celor două faze. 


În prima fază căutăm soluţia optimă a problemei (artificiale): 


xq 3x93 +2x — x4 +x =3, 
3x +3% + x% =x; +a = 4, 
x; >0, Xis X5, Xe, X3 D0 
(min) g(x*) =x +x%5 


Soluția admisibilă de bază de pornire este x= 0, 0, 0, 0, 0, 3, 4 A A 


6l 


Tabelul simplex are forma: 


P: 0 1/2 l 0 1/2  —1/2 -1⁄2 1⁄2 


Am găsit soluția optimă x= 1/2, 5/6, 0, 0, 0, 0,0 "a problemei artificiale. Soluția 


admisibilă de bază pentru problema inițială va fi: x = 1/2, 5/6, 0, 0, 0 i „iar sistemul de restricții 


cel ce rezultă din ultima etapă a algoritmului simplex pentru problema artificială. 


Tabelul Simplex corespunzător fazei a doua are forma: 
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Problema admite optim infinit deoarece vectorul P; nu are şi coordonate strict pozitive. 


Exemplul 5.4. Să se determine soluţiile optime ale problemei de programare liniară: 
xX; — 2x; +x; + 2Xg =4 
3x, — X4 +X; — 8X6 =3 
X, +X, +X4 — X; — 3x; = 2 
xj>0,j=1,6 
(min)f(x) = —2x, — 3x, — X, +4x, +5x; +8X6 


Prima şi a treia coloană a matricei A sunt coloane ale matricei unitate de ordinul trei şi 


atunci vom adăuga la restricția a doua o variabilă artificială x3. 


Vom căuta în prima fază soluția optimă a problemei: 


xı — 2x4 +X; +2X6 = 4, 
3x, — X4 +X; — 8x; +X} =3, 
Xp X3 HX4 — X; — 3X6 = 2. 
x;>0,j=1,6, x} >0 


(min)g(x*)= x3 


O soluție admisibilă de bază pentru această problemă este: 
x= 4, 0, 2, 0, 0, 0,3" 


Tabelul simplex are forma: 


4 o 0 

p? Ega 30 0 i E g ii 

Pp | 10 0 0 ı 1 so se 14) 
zen e a 0 Ap g E a de 8 — 0 

Pi Aae 0: 0 ea e 2 1 

P [0 0. si 0 -1⁄3 1⁄3 -8⁄3 14 

P |01 o 0 1 43 -43 AB -1/3 
zc;| -]ļ]ol|lo 0 0 0 0 0 4 
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Soluţia admisibilă de bază pentru problema inițială va fi x = 4, 1, 1, 0, 0,0 i 


Tabelul simplex pentru problema inițială (faza a doua) are forma: 


i aia Dă 1 4 5 8 
S Ps P4 O OPs O Ps 
P | 2 | 4 | o Ww d 2 
Pa | A 0 | 0 -13 1⁄3 -83 
Pa et e 0 0 1 43 A43 -13 
zgr | e e o 08 e DU E 


Soluţia optimă a problemei este x, = (4,1,1,0,0,0)', (min)f(x) = —12. 
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US6 Probleme de transport 


Problemele de transport sunt probleme de programare liniară particulare. Ele se 
referă la distribuirea unui produs P aflat în cantităţi date într-un număr de depozite şi solicitat de 
un număr de centre de desfacere, de asemenea în cantități cunoscute. Această problemă nu se 
rezolvă cu algoritmul Simplex din cauza dimensiunilor matricii asociate. 

Algoritmul prezentat are aceleaşi caracteristici cu algoritmulk Simplex: un criteriu 
de optim, un criteriu de intrare şi un criteriu de ieșire. Toate acestea sunt descrise în prezenta 
unitate de studiu. Este descris apoi algoritmul, modul în care se obține o soluție mai bună dintr-o 
soluţie dată. 

Studentul trebuie să înveţe mai întâi să echilibreze problema. Ulterior el trebuie sa 
obțină o soluţie nedegenerată printr-una dintre cele două metode uzuale (metoda diagonalei şi 
metoda costului minim). După verificarea optimalității soluției, urmează construirea unei soluții 
mai bune. 

Calculele nu sunt complicate, dar algoritmul pune în general probleme studenţilor. 
Sunt necesare două ore pentru însuşirea noţiunilor teoretice şi încă patru ore pentru aplicarea 
cursivă şi ușoară a criteriilor. 

Titlurile bibliografice [3] — [9] conţin tratarea subiectului în manera pe care și noi 


am abordat-o. 


Să considerăm că un anumit produs este stocat în m depozite D; în cantitățile a; , i =1,m. 


19 

Produsul este solicitat de n centre de consum C; în cantităţile b; , j= ln. Costul transportului 
unei unități de produs de la depozitul D; la centrul Cj este a; u.m. Se cere a se determina 
cantitățile x; de produs care urmează a fi transportate de la depozite la centrele de consum astfel 


ca disponibilul a fie epuizat, în fiecare depozit, cererea să fie satisfăcută exact, în fiecare centru 


m n 
de consum, iar costul total al transportului produsului să fie minim. Se presupune că Xa; =} b; 
i=l a=l 


„ adică disponibilul din depozite este egal cu cererea totală a centrelor de consum. O astfel de 
problemă este denumită problema de transport echilibrată. 
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Restricţiile problemei sunt din nou de trei tipuri: restricţii asupra disponibilului şi cererii, 


restricţiile datorate sensului economic al variabilelor şi restricţia minimizare a costului. 


Datele modelului economic le prezentăm într-un tabel de forma: 


C 
CGIC ai 
D 
C11| C12 
D, a] 
XI  |XI2 
C21| €22 
D- a2 
X21 [X22 
Cil Ci2 
D; ai 
Xil Xi2 
Cm1 Cm2 
Dm am 
Xml |Xm2 
da; 
j 


Traducerea algebrică a celor trei tipuri de restricții conduce la: 


l mi T 
(minfE) = 32 32 CjXjs X5 Xip X25 mn 
i=l j=l 


j=l 
aN 
x; >0, i=l,m; j=l, 


(6.1) 


Observație. Dacă disponibilul este mai mic decăt cererea sau mai mare, pentru a 


echilibra problema se introduce un depozit fictiv sau un centru fictiv cu cantitatea (fictivă) 


existentă sau cerută astfel încăt problema să devină echilibrată. Costurile unitare de transport 


pentru depozitul fictiv sau centrul fictiv se consideră nule. 
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Teorema 6.1. Problema echilibrată de transport echilibrată (6.1.) admite cel puţin o 
soluție admisibilă, iar o soluţie de bază admisibilă de bază are cel mult m+n-l 


componente diferite de zero. 


m n a.b A 
Demonstrație: Cum Xa, = Xb, =S, considerând x, = T avem: 
i=l j=l 


S E E E 
iao SS S păi 

x = Sici a=- S=., j=], 
d ds Tg dai gS] 


şi cum X, 20, ajungem la concluzia că X(x,) verifică sistemul de restricţii şi condiţiile 


nenegativitate, fiind deci soluţie admisibilă a problemei de transport. Putem observa că din cele 
m+n restricţii doar m+n-l sunt independente. Este suficient pentru aceasta, de exemplu, ca 
din suma primelor m restricţii, să scădem suma următoarelor n—1 şi vom obţine atunci ultima 
restricție. Pe de altă parte ştim că numărul componentelor nenule ale unei soluţii de bază 
admisibile de bază este cel mult egal cu numărul restricţiilor independente şi în felul acesta 


teorema este demonstrată. 


6.1. SOLUȚII ADMISIBILE DE BAZĂ PENTRU 
O PROBLEMĂ DE TRANSPORT 


Soluţiile optime pentru o problemă de transport vor fi căutate în mulţimea soluţiilor 
admisibile de bază şi din acest motiv vom pune în evidență două metode prin care putem găsi 


aceste soluţii. 


a) Metoda diagonalei. Se începe cu determinarea componentei corespunzătoare căsuţei 


din stânga sus a tabelului deci x; 
Vom lua x, = min(a,b,) şi vom atribui valoarea zero tuturor variabilelor de pe aceiaşi linie sau 
coloană cu x în funcţie de următoarele situaţii: 


1. dacă a, <b,, vom lua x =4, x; =0, j=2,n 
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2. dacă a, > b,, vom lua x, =b, x =0, i=2,m 

3. dacă a, =b,, vom lua x, =a, =4, iar xp sau x, =0, 
şi restul componentelor de pe linia şi coloana lui x, egale cu zero. În ultima situaţie obţinem 
soluţii degenerate. 

Modificăm apoi ai şi bı înlocuindu-i cu af =a; — xX, b =b — x. Procedeul se repetă 
la pasul următor pentru calculele rămase necompletate care formează un tablou cu m — 1 linii şi 
n coloane ; m linii şi n—1 coloane; sau m—1 linii şi n-—l coloane după cum primul pas s-a 
desfăşurat în situaţiile 1), 2), respectiv 3). 

Prima componentă ce va fi determinată acum va ocupa aceeaşi poziţie (colţul stânga — 
sus), în noul tabel. Procesul se termină în cel mult m +n -—l paşi, la fiecare pas determinându-se 
complet o linie (situaţia 1), o coloană (situaţia 2), sau o linie şi o coloană situaţia 3. La fiecare 
pas se determină o singură componentă nenulă. Componentele nebazice nu se completează în 
tabel, căsuţele corespunzătoare rămânând libere pentru a nu se confunda componentele nebazice 


cu eventualele componente bazice nule. 


b) Metoda costului minim. Această metodă este întru-totul similară celei a diagonalei, 


cu singura diferenţă că se va determina mai întâi componenta x;; pentru care c; = min Cy , luând . 
x, = min(a;,b;) . 


Exemplul 6.1. Să se determine o soluţie admisibilă de bază pentru problema de transport 


ele: C [œ [Ca T. 
D ^ 


Dı 35 
D, 2 2 4 5 40 
D; 3 4 2 3 20 


3 
b | 50| 10 | 25 io. ji 


68 


a) Metoda diagonalei. Problema este echilibrată. Vom căuta o soluţie cu 


m+n—1l=—3+4—1=6 


35,0 


40,25, 15, 
0 


20, 10, 0 


Pasul întâi: deoarece  min(a,,b,)>a=35, punem conform situației 1) 
E g A A . x 7 
x% = 35, Xp =X%3 = Xa =0. Modificăm pe a şi b înlocuindu-i cu a, =aj—xX4=0, 


Pasul al doilea: se pleacă de la tabelul format prin înlăturarea primei linii,: deoarece 
min(a,,b/)>b/=15, punem (conform 2)) x, —15, x, =0 (x, componentă nebazică). 

Pasul al treilea: a rămas de completat tabelul format din liniile a doua, a treia şi ultimile 
trei coloane. Cum min(a,,b,)>b,=10 punem x, =10, x =0 (x, componentă nebazică) 

Pasul al patrulea: a rămas de completat tabelul format din ultimile două linii şi coloane. 
Cum min(a),b.)>a”=15, punem (conform 1)) x =15, x, =0(x,, componentă nebazică). 
Modificăm pe a, şi b, înlocuindu-i cu a= a} — x =0, 


Pasul al cincilea: a rămas de completat tabelul format din ultima linie şi ultimile două 


coloane. Cum min(az,b))F10 punem x =—10. Modificăm pe a, şi b, înlocuindu-l cu 


l o — N =: 
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Pasul al şaselea: obligatoriu vom pune x, = 10. 
S-a obţinut o soluţie admisibilă de bază: 


X 1 = 35, x72 =0,x3 = 0, a = 0, x = 15, x =10, x3 = 15,34 = 0, x = 0,X3 =0, 


costul total al transportului fiind 300 u.m. 


b) Metoda costului minim. Vom considera din nou tabelul simplificat: 


25 110 35, 10,0 
2 2 4 5 
30.110 40, 30, 0, 
3 4 2 3 
20 20, 0 
10 25 10 
0 0 0 0 


Pasul întâi: minc;j=c = C3 = Co = Co = C33 =2. Putem determina ori care dintre 


componentele x2, X3, X213 X223X33. Luăm x„=min(a,,b,)=10, xp =X; =0 (componente 


. a ; - 7 
nebazice). Înlocuim pe a, şi b, cu aj =a, — x =30, b, =b, — xX» =0. 


r N x EA Po a 2 é a 1 R — E 
Pasul al doilea: min c;j> c= Co C33= 2. Luăm x,„=min(a3,b,)=30, x3 = x34 =0 


: a . . Li 
(componente nebazice). Înlocuim pe a; şi b, cu af =a} — x =0, b =b —x =20. 
Pasul al treilea: min c= C43= C33= 2. 
Luăm x = min(a,,b) = 25, x, =0 (componentă nebazică). Inlocuim a, şi b, cu 


7 7 


Pasul al patrulea: min c;;> Cu4= 034 C33= 3. 

Luăm x = min(a3,b4) =20, x, =0 (componentă bazică), x}, =0 (componentă 
nebazică). Înlocuim a; şi b! cu a/=— a, — x =0, b!"=b!— x; =0. 

Pasul al cincilea: obligatoriu x4 = 10. Am obținut pentru problemă o soluţie 
admisibilă de bază degenerată (x, = 0 componentă bazică). Ea este indicată în tabelul de mai 


sus.Costul unui asemenea transport este 220 u.m. Metoda costului minim ne conduce de cele mai 
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multe ori la o soluţie admisibilă de bază mai bună decât metoda diagonalei, în sensul că 


realizează o valoare a cheltuielilor de transport mai mică. 


6.2. SOLUȚII OPTIME PENTRU O 
PROBLEMĂ DE TRANSPORT 


Soluţiile optime (soluţii admisibile ce realizează minimul funcţiei obiectiv) ale unei 
probleme de transport le vom căuta în mulțimea soluţiilor admisibile de bază nedegenerate. Să 


considerăm că am determinat pentru problemă o soluţie admisibilă de bază nedegenerată: 
x=(x; >! . Pornind de la aceasta să vedem cum putem obţine o soluție admisibilă. 
Convenim ca o pereche de indici (i, j) să o numim celulă. Când vom spune celulă liberă 


înţelegem că ea se referă la o componentă nebazică (secundară), adică nu se găseşte printre cele 
m+n-1. Celule ocupate vor fi cele care se referă la componentele bazice (principale). Să 


considerăm o celulă liberă (i, j). 
Definiţia 6.1: Un şir de celule care începe cu celula liberă (i,j) şi se termină cu 


aceasta, conţinând în rest numai celule ocupate şi anume câte două din aceiaşi linie şi 


coloană se numeşte ciclul celulei libere (, j). 
Cel mai simplu ciclu al unei celule libere (i,j) se prezintă sub forma de mai jos unde în 


dreptul celulelor am scris şi componentele corespunzătoare ale soluţiei. 


Ai 


G6, ) > (ij) 


Gi Je f) 
Să atribuim con Xi j ïX; (componenta 1 X; j ), o valoare 0>0.Vom obține atunci 


un vector care are m+n componente diferite de zero, acesta trebuind să verifice sistemul de 


restricții. 
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Acestea sunt verificate dacă vom scădea şi aduna succesiv valoarea O la componentele 


care au intrat în ciclul celulei libere (i, j). Indicăm pe acelaşi ciclul acest lucru. 


ra) x. —0 


i 
(1) Gui) 


(iJe (iJi) 


Xij O 0 NT 4 


Pentru a fi verificate condițiile de nenegativitate va trebui ca 0 să fie luat astfel încât 
Xij —0 >0 şi x; —9>0. 
Odată îndeplinite aceste condiţii putem spune că am obţinut o nouă soluție admisibilă 


x(x,) unde: 


i ai 0, xi = Xij =, Xij DĂ], T (6.2) 


Xij =X; 0, x, = Xœ pentru k =i}; j= j, ji 
Dorind acuma ca noua soluție admisibilă să fie şi admisibilă de bază este suficient să 
considerăm: 
0 = min(x, Xj) (6.3) 
În acest fel se determină care componentă bazică devine nebazică. Relația (6.3) reprezintă 
forma criteriului de ieşire din algoritmul pentru găsirea soluțiilor optime a unei probleme de 
transport. 


Să vedem acum care sunt condițiile ca soluția admisibilă de bază x să fie soluție optimă. 


Să evaluăm diferența f x -f x . Vom avea: 


m n 
fx-fx =E Ecke Xhe — Xe = cd +e 0 -e0 +e; 
— E eee 


sau 
f x —f x =08; (6.4) 


unde 
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dia (6.5) 


id 1J 


8, = —C; +, =c 


care reprezintă suma algebrică cu semn alternat a costurilor celulelor corespunzătoare ce au intrat 
în ciclu şi care începe cu —c;. 

Din (6.5) se observă că dacă pentru toate celulele libere (i, j) vom avea 5, <0, şi aşa cum 
0 >0, soluţia admisibilă de bază x(x,) considerată reprezintă soluţia optimă a problemei de 
transport. Rezultă în felul acesta că criteriul de optimalitate este 

ô; <0 (6.6) 

pentru toate celule libere (i, j). 

Tot din (6.4) se constată că dacă există 5, > 0, soluţia x poate fi îmbunătăţită. Pentru 


aceasta fie: 
Ojj > Max Öke (6.7) 


Ske 20 
Atunci componentei nebazice x;îi vom atribui o valoare 0 >0. Relația (6.7) reprezintă 
forma criteriului de intrare în algoritm. Mărimile 5, joacă rolul diferențelor z; —c; din problemele 


de programare de care ne-am ocupat. 


Sintetizând cele de mai sus în practică procedăm în felul următor: 


a) Se determină o soluţie admisibilă de bază nedegenerată prin una din cele două metode 
prezentate. 


b) Pentru toate celulele libere se calculează cantităţile 5, . 

- dacă toate cantităţile ô sunt mai mici sau egale cu zero, atunci soluţia este optimă; 

- dacă există 6, >0 se trece la etapa următoare. 

c) Se aplică criteriul de intrare. Vom atribui o valoare 0 >0 componentei x, (nebazică) 


cea pentru care avem îndeplinită relația (6.7). 


d) Se consideră ciclul celulei x, , se introduce 0 în ciclu şi se aplică criteriul de ieşire. 


ij? 
Se determină valoarea lui 9 după relaţia (6.3) indicându-se deci care componentă bazică devine 
nebazică. 


e) Etapele algoritmului se reiau până se găseşte soluţia optimă a problemei. 
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Observaţii: 

a) Prezenţa semnului de egal în criteriul de optim şi de intrare în bază, indică existența mai 
multor soluţii optime ; 

b) Pe parcursul aplicării algoritmului pot să apară soluţii admisibile degenerate care pot 
conduce la un fenomen numit fenomen de ciclaj. Pentru evitarea apariției soluţiilor 
admisibile degenerate se întrebuințează aşa-numita metodă a perturbării. Aceasta 


constă în înlocuirea cantităților a, cu a; +£, şi cantităţii b, cu b, +me (unde reprezintă 


o cantitate foarte mică, strict pozitivă). La finalul algoritmului, după determinarea soluţiei 
optime a problemei peturbate, se face e=0, obținându-se soluţia optimă a problemei 
inițiale. 

c) Dacă la aplicarea criteriului de intrare sunt mai multe mărimi ô; pozitive maxime egale 


putem lua oricare dintre acestea. 


Exemplul 6.2. Să se determine soluțiile optime ale problemei de transport din exemplul 6 
Pentru a fi siguri că nu apar soluții admisibile de bază degenerate aplicăm de la început 
metoda perturbării, adică considerăm: 
a,=35+e, a,=40+e, a,=20+e, b, = 50, b, =10 b} =25, b} =10,+3e 
O soluţie admisibilă de bază determinată prin metoda costului minim din tabel este dată în 


tabelul: 


C 
N C C | G 
4 2 2 
Di 10 125 
D, 2 2 3 
40+e 
D, 3 4 J 
10—e 
bi so 10 | 25 


Să găsim cantităţile 5, pentru celulele libere. Considerăm ciclul celulei libere (1,1) şi în 


dreptul celulelor scriem costurile. 
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(3, ha 8,4) 
3 3 


Vom avea: ô=—4+3—3+3=-1 


Analog: 


SE e ei 
SARIA 2E 
ei E 
Di) ae 2) 
—2+2—3+3=0 


a n n n o 
Ww” N X Ww” N 
MII 


Deoarece toate cantităţile ©; sunt mai mici sau egale cu zero am găsit soluţia optimă a 
problemei. Făcând s=0 găsim: 

Xx =0, x2 =10, x, = 25, x4 =0, x1 = 40, x» =0, x =0, Xa =0 

X3, =10, x; =0, x33 =0, x34 =0, (min)f(x)=210 u.b. 
Deoarece 5, =0 , algoritmul poate continua atribuind componentei x, (nebazică), o valoare 


0 > 0. Considerăm ciclul celulei (3,3) şi introducem pe 6 în ciclu pentru a satisface restricțiile. 


25—98 e +0 


(1,3) —————————————> (1,4) 


(3,3 je (3,4) 


Aplicăm criteri 0 ieşire luând 0=m 10+28—0 )=10e. 
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Componentele noii soluţii sunt date în tabelul: 


Făcând £=0 găsim soluţia optimă. 


x =10, x32 20, x 210, x, =0, (min)f(x) =210 u.b. 


Exemplul 6.3 Să se determine soluţiile optime ale problemei de transport: 


Problema nu este echilibrată. Vom adăuga un depozit cu costuri zero. 


Vom avea deci: 


C2 C3 aj 
2 3 50 
4 l 25 
75 
20 30 90 


D Ci C2 C3 ai 
D: 2 50 
D2 } 25 
D3 X 15 

90 
bj | 40 | 10 | 25 [og 


76 


O soluţie admisibilă de bază determinată prin metoda costului minim din tabel este 
prezentată în: 


Componentele soluţiei au fost determinate în ordinea: 
Xa = 5 dop 295 4 29. ăia = 20, 

X3 =5, X2 => 0, X = 0, x33 =0, x33 =0 

Pentru celulele libere calculăm cantitățile $,;. 


—3 +1—3 +1 = —4 
—4 +2 —3 +1 = —4 
õp =—0+2—1+0=1 
833 =—0+3—1+0=2 


[2] 
N 
N 
III 


Soluţia nu este optimă. Avem că 6, = max ôjj=2. Considerăm ciclul celulei (3,3) 


6; >0 


componentei x, (nebazică) îi atribuim o valoare 0 > 0 şi introducând pe 0 în circuit avem: 


25+0(30) 5—0(0) 
(1 > (1,3) 


15—0(10) (5) 


După criteriul de ieşire luăm 0 = min(1,15) = 5. Noua soluţie este prezentată în tabelul: 
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SI AE AE IE 
De a “|. 50 
D, Ie a 
Dig aler aas 
b; | 40 | 20 | 30 [SS 


Pentru celulele libere vom avea: 


õ13 

21 
ôn 
ò 


32 


—3+1—-0+0=-2 
—3+1—0+0=-—2 
—4+2—1+0—0+1=-—2 
—0+2—-1+0=1 


Soluția nu este optimă. Atribuim componentei x,» o valoare 0 >Q. 


30 +0(40) 


20—610) 


(1, D) ——————> (1,2) 


(3, L je—— 8,2) 


10—80) 


Noua soluție este prezentată în tabelul: 


010) 
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Această soluţie este optimă deoarece avem: 
63 = — 1l, = —36p = —3,6, =—l. 
Soluţia problemei iniţiale va fi: 
x, = 40, x2 = 10, x3 =0, x1 =0, xx =0, x = 25 . (min)f(x) = 8Su.b. 


Se observă că cerințele centrelor C2 şi C3 nu vor fi satisfăcute în totalitate. 
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US7 Serii de numere reale 


Problema sumelor infinite (a seriilor) este adesea întâlnită în asigurări şi finanțe. 
Rezolvarea ei presupune a stabili dacă o astfel de sumă există şi în caz afirmativ, care este 
valoarea ei. Nu este o problemă ușor de rezolvat pentru ca nu toţi studenţii cunosc bazele analizei 
matematice (mai precis teoria şirurilor de numere reale). Prin activităţile tutoriale completăm 
acest lucru, după care definim noţiunile fundamentale (serie convergentă, serie divergentă, suma 
unei serii convergente, criteriu de convergență etc.). Se insistă asupra faptului că suma unei serii 
convergente se poate aproxima, după ce (printr-un criteriu dintr-o listă destul de bogată) am 
stabilit natura ei. 

Studenţii trebuie să înţeleagă mai întâi ce înseamnă serie convergentă analizând 
exemplele clasice: seria geometrică, armonică, armonică generaliată, telescopică. Urmează apoi 
reținerea criteriilor care se aplică oricărei serii (criteriul general de convergenţă, criteriul general 
de divergență). După această etapă urmează criteriile pentru fiecare tip de serie (criterii pentru 
serii cu termeni pozitivi, criterii pentru serii cu termeni alternaţi, criterii pentru serii cu termeni 
oarecare). 

Deşi este vorba de un singur subiect teoretic, sunt necesare aproximativ 8 ore 
pentru însuşirea lui: două ore pentru noțiunile teoretice generale, două ore pentru seriile 
remarcabile, două ore pentru seriile cu termeni pozitivi şi două ore pentru seriile cu termeni 
oarecare şi pentru seriile alternate. 

Poziţiile [1] şi [3]-[6] din lista bibliografică pot fi de un real ajutor, dar pot fi utile 


si alte tratate de analiză matematică. 


x x 


Este bine cunoscut faptul că, în prezent, fenomenele economice devin din ce în ce mai complexe 
şi interdependente; mai nou, tendința de globalizare a economiei (şi nu numai) duce la o 
concurență acerbă pe mai toate piețele şi în mai toate domeniile de activitate economică, 
fenomen care se simte din ce în ce mai pregnant şi în economia românească. 

Aceste situații, cer din partea specialiştilor în economie, cunoştinţe cât mai precise în 


vederea observării, înţelegerii, modelării şi rezolvării pe baze riguros ştiinţifice a problemelor 
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complexe cu care se confruntă. Această abordare nu mai este astăzi posibilă, fără recurgerea la 
un aparat matematic cât mai adecvat. 

Modul mai puţin riguros, chiar empiric, de luare a deciziilor economice doar pe baza unei 
simple experienţe, este total depăşit, astfel încât cei care nu se vor adapta noilor condiţii şi 
cerinţe economice, vor fi eliminaţi încet, încet de pe piaţă. 

În acest capitol se prezintă, într-o manieră accesibilă studenţilor economişti, câteva 
elemente de bază din teoria seriilor numerice şi a seriilor de funcţii (cazul particular al seriilor de 
puteri). Aplicațiile practice ale acestei teorii se întâlnesc în multe domenii de cercetare din 
matematică, fizică, chimie, biologie etc. În particular, aplicaţiile cu caracter economic ale teoriei 
seriilor sunt multiple şi permit modelarea şi rezolvarea unor probleme economice cu caracter atât 
determinist cât şi aleator (probabilistic). 

Elemente şi aplicaţii ale teoriei seriilor sunt regăsite în Teoria probabilităților şi Statistica 
economică, Micro- şi Macroeconomie, Gestiunea stocurilor şi a fenomenelor de aşteptare etc. 

Pentru parcurgerea şi înţelegerea noţiunilor prezentate în acest capitol sunt necesare 


noţiuni de teoria limitelor şi a şirurilor studiate la Analiza Matematică din clasa a XI-a. 
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7.1 SERII NUMERICE 


7.1.1. ȘIRURI DE NUMERE REALE 


Definiţia 7.1.1: Numim şir de numere reale, o aplicaţie: 


f:N>R 
f(n)=a, €R 


Vom nota şirul cu simbolul: (a, )ney; ân se numeşte termenul general al șirului. 


Definiția 7.1.2: Spunem că şirul (a, )pen este convergent către un număr a €R, dacă, 


orice vecinătate a lui a conține toți termenii șirului cu excepția unui 
număr finit de termeni. 


Vom nota cu a, —a sau lima, =a. 


N—OO 


Numărul "a" se numeşte limita şirului şi dacă există este unică. Dacă şirul (a) ey nu este 


convergent, atunci acesta se numeşte şir divergent. Un şir divergent sau nu are limită sau 
limita sa este +00 sau —co. 


Exemple: 

n n? | PIN za ui l 
a) a= ; ba = 3 Ca S|I+—| sunt şiruri convergente (cu limitele 0, — 
n +1 2n +3 n 2 

respectiv e=2,7.11...; 

— n? +] 
b) a, =n ; b, = aE =Nl—n? sunt şiruri divergente cu limitele +00, +00, 

n 


respectiv —00; 
c) a, =(D",b, =(—1)" Inm? +n +1) ‚c, =sin(n +1) sunt şiruri divergente care nu au 


limită. 
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Teorema 7.1.1 (de caracterizare a şirurilor convergente) 
Şirul a, —a dacă şi numai dacă: (v)e >0.(2)n. EN astfel încât a, —a| <€, pentru 


(n >n,- 


Teorema 7.1.2. (Cauchy) 


Şirul (a, „en este convergent dacă şi numai dacă: (v)e >0, (E)n, EN astfel încât : 


au —a,| <e, pentru (Y)n >n, şi (Y)k EN”. 


Observații: 

i) Cele două teoreme, ca şi definiția, nu sunt folosite în mod curent în aplicații practice; cu 
ajutorul lor se pot demonstra rezultate a căror aplicare este mult mai simplă (de exemplu, 
teorema de convergenţă a şirurilor monotone şi mărginite, teorema cleştelui etc.). 

ii) Aşa cum se observă, aplicarea teoremei 7.1.2. nu presupune cunoaşterea apriorică a 


limitei a, în enunţul acesteia intervenind numai termenul general a, respectiv a,x; 


acest lucru este deosebit de util în aplicațiile care conțin şiruri recurente sau în care 


limita şirului este dificil de determinat. 


7.1.2 NOȚIUNEA DE SERIE NUMERICĂ. 
CONVERGENȚĂ. 
PROPRIETĂȚI GENERALE ALE SERIILOR. 


După cum se ştie, oricărei mulțimi finite fa; az, zi san) de numere reale i se poate asocia 
un număr real s=a, +a, +...+a, numit suma acestora. Se pune întrebarea firească: dacă 


mulțimea conţine o infinitate de numere reale a,,a»,...,a (cu alte cuvinte un şir de 


poe 


numere reale (a, ue ) l-am putea asocia un număr SER care să generalizeze conceptul de 


neN 
sumă finită? 
Fie şirul : (a,) 


ir ap apare day at 
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Vom asocia acestuia, un nou şir de numere reale, (S, )uen : So» Si; S25- --» Sp, --- definit 


astfel: 


Ve) 
deo | 
II 
o 
O 


S; =aọ +a 
(7.1.1) 


n 
S, =aọ +a; +a, +...+a, = Xa, 
k=0 
pe care îl vom numi şirul sumelor parțiale ataşat şirului (a, )nen- 


Definiţia 7.1.3. Numim serie numerică, de termen general a, expresia 


00 def 
Ža, = lim S, =a Pay +a, na +... (7.1.2) 


n=0 n—>0 


00 
Žan ) sau 
n=l 


Observaţie: [n unele cazuri numerotarea termenilor seriei poate începe de la n = ( 


o0 
de la n = nọ oarecare fixat | X 3 (vezi şi proprietățile seriilor numerice). 


n=ng 


OO 
Definiția 7.1.4.: Spunem că seria >. a, este o serie: 


n=0 


a) convergentă (C ), dacă şirul (S en este şir convergent; 


b) divergentă (D) , dacă şirul (S, ue este şir divergent. 


Cu alte cuvinte problema convergenţei unei serii numerice este echivalentă cu problema 


convergenţei unui şir de numere reale şi anume, şirul sumelor parțiale (S, ue - 


Observaţii: 


OO (8,9) 
i) Dacă S, — S, vom spune că seria >. a,(C) şi are suma S (notăm: Xa, =S); 
n=0 n=0 


n=0 


OO e9] 
ii) Dacă S, —> oo, spunem că seria } a, (D) şi are suma +œ notăm > a, = +c0o0); 
n=0 
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OO 
iii) Dacă (S )aey nu are limită, spunem că > a„(D) şi nu are sumă. 
n=0 


In studiul seriilor numerice, trebuie cercetat (ca şi în cazul şirurilor) două aspecte: 


- natura seriei (dacă este convergentă sau divergentă) 


- suma seriei (în caz de convergenţă) 


Dacă pentru şiruri stabilirea convergenţei şi determinarea limitei acestora este în cele mai 
multe cazuri una şi aceeaşi chestiune, în cazul seriilor situaţia se schimbă radical; de cele mai 
multe ori este relativ uşor să stabileşti convergenţa seriei, dar foarte greu, ba câte odată chiar 
imposibil de determinat suma acesteia. 

De obicei, cu metodele specifice analizei numerice, se determină o aproximaţie S’ a 
sumei S a seriei (S ~S). 

Pentru a se cerceta natura unei serii, nu se poate folosi întotdeauna definiţia, ci adesea se 
utilizează condiţii echivalente cu aceasta sau numai condiţii suficiente de convergență numite 


“criterii de convergenţă”. Vom prezenta unul dintre aceste criterii în cele ce urmează. 


Teorema 7.1.3. (criteriul general de convergenţă a lui Cauchy) 


Seria D a, este (C) dacă şi numai dacă: (V)e > 0,(3)n, EN astfel încât: 
n=l 


au +a,p Fana] <£, pentru (v)n >n., (Vk EN. 


Această teoremă se obține ca o consecință a Teoremei 7.1.2, considerând în loc de şirul 


(a, en » Şirul sumelor parțiale ($, ) Într-adevăr, conform (7.1.1) avem: 


neN * 


Sark = Sa 


= |(a, Faaa Fany Fetan) ACU fas hdi) 


any Fe: Pânza 
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7.1.3. EXEMPLE REMARCABILE DE SERII 


1) Seria geometrică 5 aq; cu apq ER 
=0 


Termenul general al acestei serii are forma a, =aq" (*) şi depinde de doi parametri 


(constante reale) a, q. Evident, este logic ca natura seriei şi suma acesteia în caz de convergenţă 


să depindă de aceşti parametri. Astfel, dacă: 


1) a = 0, atunci a, =0, (Y)q €R, deci seria este convergentă şi are suma S=0. 
ii) a z 0, avem: 
5) > fi 
S, Sao +a; +a, t...a, =ar+ragrag 1+...+ag 


=a(l+q +q" +...+q") 
deci S, reprezintă suma unei progresii aritmetice cu rația q (de aici şi numele 


seriei) şi conform rezultatelor de la acest capitol (vezi manualul de Algebră, clasa 


a X-a) vom avea: 


E3 n+l 
d d = jel 

S, = 1—q 
a(n+1);  q=l 


Trecând la limită în relaţia de mai sus, obţinem: 


a 


; ge(—Ll) 
1—q 
lim S, =œ ; qe[l,00) (depinde de semnul lui a) 
(4) ; qE€(~%,-—1] 


cu alte cuvinte, pentru: 
De e ăi a a 
q€ —1,l => $, aq =—— [seria converge cu suma S=——— 


a20 1—q l—q) (11.3) 
qE€(—œ,1]U[, cœ) => 5 aq"(D) cu suma S = zoo sau nu există 
n=0 
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Observaţie: De fapt, această serie cuprinde o infinitate de serii particulare (pentru a, q 


având valori fixate) a căror natură o cunoaştem; foarte important este că, în caz de convergenţă, 


ştim şi suma seriei. 


Exemple: 
e 1 (1) fă l. cei 
a) > za => |>| „este serie geometrică cu a =1, q = 3 şi conform (7.1.3) rezultă că: 
n=0 n=0 
œ% l 1 1 l 
=1+-4 una Fae DAC 
e De ph O) 


b) 5S asi = 5 5] „este serie geometrică cu a=l şi q=- şi conform (7.1.3) 


rezultă că: 
a (—1)” l 1 (—D" 3 
d = LT +. = (0) 
dp PIB a 
ER S [ BA)I e SIRE l 
0) D z7 = 05. ala este serie geometrică cu a= 5, 1>2 şi conform (7.1.3) seria 
n=0 n=0 
SS gn 
diverge |în acest caz X` z7 = +00]. 
n=0 


d) F e J2 i =5 (—1) =o i este serie geometrică cu a=-—l, q=—N2 şi deci 
=0 =0 


seria este divergentă; mai mult suma nu există. 


n=l 


~ l 
2) Seria armonică generalizată aia € R| 
n 


Evident, ca şi în şirul seriei geometrice, natura (şi în caz de convergență suma) seriei 
armonice generalizate, va depinde de valoarea parametrului a. Aplicând, Teorema 4.2 se poate 
arăta că, seria: 


(C) a; a>l 


pe =] (7.1.4) 


1 n 
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Observaţii: 


i) In caz de convergență (a > 1), nu există o formulă generală de stabilire a sumei S în 


funcţie de parametrul a, asa cum s-a stabilit în cazul seriei geometrice; 


ii) Evident, în cazul în care seria diverge (a <1), suma seriei este S= -+00 (ca serie cu 


termeni pozitivi). 


Exemple: 


a) 


b) 


c) 


forma 


œ% | œ | RU | 1 
> — (D),deoarece a=1 —=1+- +- +...+— +... = +00) 
n n= N 2 Ca n 


se al e l 1 1l 1 

— (C),doarece a=2>1 i SE E Ma a E Bai A, 
z n? o 5 n? 2 7 n? | 
1 


3/2 
n 


œ 1 
= 3 — (D), deoarece a=2 SI 
n=l 3 


n3 


K | 
Seria >D — se numeşte seria armonică (simplă) sau seria lui Riemann, deoarece 


n=1 n 
dau tă a ue Sf as i 32 1 1 
termenul general an este media armonică a vecinilor săi, adică satisface relația — = + 
an an-ı an+ 
z . |2 1 
3) Seria telescopică | 5 
n=1 n(n +1) 
z ; DR y de 1 
Această serie este convergentă. Într-adevăr scriind termenul general a, = (i) sub 
n(n + 
1 1 . roii 3 
(=) a, sa obținem termenul general al şirului sumelor parțiale (S ax de 
n 


forma: 


S, =a +a, +...+a 


Lo 1 1 1 
n=- + + +... + ———— = 
b2 23 Aa n(n+1) 


1 1 1 1 1 1 1 
in = (după simplificări 
; i J+ i): 2 +] (dup p ) 
l 
n+1 
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AC Sa ae ta 8 e dia h- T AEE 


nı n(n+l) n-o n-o n+ nl n(n +1) 
COP 1 ; iii 1 ; 
Dacă notăm cu a, =—, atunci termenul general al seriei a, = ———, se poate scrie 
n n(n +1) 
sub forma: 
aa = On = (7.1.5) 


Seriile numerice a căror termen general an au forma de mai sus, se numesc serii 


telescopice şi se poate observa uşor ca seria Sa, (C) > (a, an (0). 


n=l 
Mai mult seria va avea suma S=q;—-@, unde a, este primul termen al şirului 


(0, en» lar a= lim a. 
n—% 


Vom enunţa în continuare câteva din proprietăţile generale ale seriilor numerice, care se 


pot demonstra relativ uşor, ţinând cont de definiţia seriei şi a operaţiilor cu şiruri convergente. 


7.1.4. PROPRIETĂȚI GENERALE ALE 
SERIILOR NUMERICE 


00 
(Pı) Dacă unei serii > a 
n=l 


i se adaugă sau i se suprimă un număr finit de termeni, 


n? 


natura acesteia nu se schimbă. 


Observaţie. Evident, în caz de convergenţă, suma S a unei serii se modifică 
corespunzător prin adăugarea (sau scăderea) a sumei termenilor finiţi care se adaugă (sau se 


suprimă). 


00 
(P2) Sacă unei serii convergente > a,, i se asociază termeni în grupe finite, cu 
n=l 


păstrarea ordinii lor, atunci se obține tot o serie convergentă şi cu aceeaşi sumă. 
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Observaţii: 


i) Proprietatea de mai sus, poate fi scrisă şi sub forma: 


00 
(P;) Dacă Sa, =a +a, +...+a, +...=S, atunci şi: 


(a +a, +... +a) + ana tean) PC aa Piata te ag) HSS 


ii) Asociind termenii unei serii divergente, în grupe finite cu păstrarea ordinii se 


pot obţine serii convergente. 


Exemplu: 


Fie seria 5 (D° =1—1+1—1+...+(—D° +... (termenul general fiind a, =(—1)") 


n=0 
Atunci: 


deci: S», =1 Şi Sau =0 ; evident şirul sumelor parţiale (S, )aen diverge deoarece are două 


subşiruri care converg la limite diferite. Rezultă că seria X (—1)” (D) şi nu are sumă. 
n=0 
Dacă asociem termenii în grupe de câte doi (sau de câte patru, şase, etc.) obţinem: 
(—1)+(—-D+(-—D)+...(—D+...=0+0+...+0+...=0 


deci noua serie obţinută converge şi are suma S = 0. 


00 00 
(P3) Dacă seriile > a, şi > b, sunt convergente cu sumele S; respectiv S2, atunci seria 
n=l n=l 


> (a, +b,) converge şi are suma Sı + S2. 


n=l 


Cu alte cuvinte dacă: 


Š a, =S, È b, =S, = È (a, +b,)=S; +S, (7.1.6) 
n=l n=l n=l 


00 00 
v * . .. . . v . A 
(P4) Dacă NER , atunci seriile X, a, şi > Aa, au aceeaşi natură; mai mult în caz de 
n=l n=l 


00 (9.9) 
convergenţă, dacă > a, =S atunci > Aa, = AS. 
n=l n=l 
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OO 
Teorema 7.1.4. Dacă seria ` a, (C), atunci lim a, =0 
n=1 N—OoO 


Demonstraţie: 
. * . 
Fie S, =a; +a, +...+a,,(v)n e N . Deoarece seria converge, avem: 


© a, = lim S, =S () 
n=1 n—0 
Dar a, =S, —Sa trecând la limită în această relație obținem: 
lim a, = lim (S, — Sp) = lim S, — lim Spa eR S=0 
n—>0 n—>0 n—>0 n—>0 
q.e.d. 


Observaţie: Condiţia din teorema de mai sus (a, — 0) este o condiţie necesară dar nu şi 
suficientă pentru convergența unei serii; deci, dacă termenul general a, > 0 acest lucru nu 


implică obligatoriu că seria > a, (C). Ca exemple se pot da cazurile particulare ale seriei 


n=l 


L (D) chiar dacă a, = i — 0. 


armonice. Astfel > 
n 


n= N 
O formulare echivalentă a teoremei 7.1.4. des folosită în aplicaţii practice, este 


Teorema 7.1.5. Criteriul general de divergență (C.G.D) 


OO OO 
Fie seria > a, a, ER.Dacă lim a, #0 atunci seria > a, (D). 
n=l N—00 n=1 


Exemple: 
a) > —— z D), deoarece lim a, = lim n _l 
n=l 2n n—>0 n—% 2n + 1 2 


b) S a (D), deoarece: 
n 


lim a, = lim (=) = lim (1+) =e% 2,71...+0 
n 


n= n= n n= 
2 2 
e l+n-—n l+n—n 
c) > ——— (D), deoarece lim a, = lim ——— = —oco #0 
na 2n n-o n-o 2n 


9] 


In concluzie, dacă: 


i) lima,z0>Voa, (D) 


n=l 
ii) lima, =0 = Xa, poate fi (C) sau (D) (trebuie făcută o analiză suplimentară în 
n—oo n=] 


acest caz) 


7.2 SERII CU TERMENI POZITIVI 


Definiţia 7.2.1. Numim serie cu termeni pozitivi , o serie de forma: 


Sa, cua, >0, (DneN. 
n=l 


Teorema 7.2.1. (de caracterizare a convergenţei seriilor cu termeni pozitivi) 


00 
Seria cu termeni pozitivi Xa 
n=l 


as A, Z0 converge dacă şi numai dacă şirul sumelor 


parţiale (S, en este mărginit. 


Demonstraţie: 
Avem 


Saa —S, = (a +a, +-.-+anu)— (a +a, +--- anu) any Z0deci Sau 25, adică 


n+l 


şirul (S, )aen este monoton crescător. 


Dacă (S, nem este şi marginit (superior), rezultă că este convergent (conform teoremei de 


convergenţă a şirurilor monotone şi mărginite) deci Xa, (C). 
n=l 


Reciproc, dacă Za, (C) => (Sa en (C). Dar un şir convergent este obligatoriu mărginit, 


n=l 
deci teorema este demonstrată. 


q.e.d. 


92 


Observaţie: În aplicaţiile practice, Teorema 7.2.1. nu este folosită foarte des deoarece, în 
general, este relativ complicat a se demonstra mărginirea şirului (S, )aen ; întradevăr, trebuie 


demonstrată relaţia: 


<M, (WneN cu M>Oconstantă care nu depinde de 


|S,|=|a, +a, +...+a, 


neN.. 
Mult mai utile sunt următoarele rezultate (numite criterii de convergenţă), care pot fi 


clasificate în două categorii: criterii de comparație şi criterii cu limită. 


ID) Criteriul de comparaţie de specia I (cu inegalităţi simple) 


Fie seriile: Sa, şi X b, care satisfac relaţiile 0<a, <b,, (WneN*. 


n=l n=l 


Atunci: 


i) Dacă Sb, (O > Sa (© 
n=l n=l 


ii) Dacă Sa, D) > Sb. (D) 
n=l 


n=l 


œo 1 | 
Exemplu: Fie seria >. 2; Termenul general a, = z satisface: 
n=1 N n 
EIRA . not 
pal a ella LA er apt 
n-:n-n....n nnn n nn n 
Dar: 
OO 09 2 S] 1 , ER À 2 1 
2 b, =} — =2} — (O) (ca serie armonică generalizată > — cu a=2>1) 
n=l n= N n= N n= N 


; E eram esa Să À 
Atunci conform punctului i) din criteriu > —— (C). 


n=l 
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II) Criteriul de comparaţie de specia a III-a (cu limită) 


O (9,9) 
Fie seriile Ya, şi > b, cu a, b, >0. Dacă există: 


n=l n=l 


lim > =k € (0,00) (7.2.0) 


n>a 


atunci seriile au aceeaşi natură. 


A n e 1 
Exemplu: Fie seria >, cu a, = > 
n=1 2N+3 2n+3 


9.9) OO 
Luăm ca serie de comparaţie, seria armonică (simplă) > b, =} „ despre care ştim 


1 
n=l n=l n 
că diverge. Deoarece: 


1 


E EE R a E EL ra 
n= b, n= 1 n—o 2n + 3 2 


= k=- € (0,00) deci conform criteriului, 


n 
. p9 . . Ss) . . 
seria a, are aceeaşi natură cu seria >. b, , deci diverge. 
n=l n=l 
Observații: 


i) Criteriul de comparație la limită, tratează şi cazurile când k=0 sau k = -+oo, dar s-a 
renunțat la prezentarea lor pentru simplificarea expunerii. 
ii) Aplicarea criteriilor de comparație, nu este în general foarte simplă, deoarece trebuie 


satisfăcute simultan următoarele două condiţii: 


a) alegerea seriei de comparație (de exemplu, Sb, ), funcție de seria dată (de 


n=l 


OO 
exemplu, a, a cărei convergență vrem să o stabilim); 
n=l 


O 
b) trebuie cunoscută aprioric natura seriei >. b, (convergentă sau divergentă 
n=l 


Prima condiție este influențată foarte mult de abilitățile matematice ale rezolvitorului, iar 


cea de a doua de “experiența” sa. Pentru această din urmă condiție, seriile geometrică şi 
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armonică pot furniza o infinitate de serii particulare (a căror natură o cunoaştem) posibil de 
folosit ca serii de comparaţie. 


iii) In schimb, următoarele criterii (numite criterii cu limită) se aplică mult mai uşor, 


DO 
formularea acestora depinzând exclusiv de termenul general an al seriei 3a a 


n=l 


cărei natură vrem să o determinăm. 


III. Criteriul raportului (D” Alembert) 


: e 08 ; . a : s 
Fie seria Xa ,, a, >0 şi: lim == =) € 0,00 . Atunci, dacă: 


n=l n-o 24 


DA<l> Xa (O 
n=l 
i)\>l=> 3 a, (D) 
n=l 
iii) = 1; caz de nedeterminare (nu se poate preciza natura seriei) 
Exemple: 


n! n! i 
— ; cu a, = ră respectiv a, 


D 5 


n=1 N 


_ (+l)! ; 
= (na . Atunci 


| n l. n 
1 = lim 22 = lim (n+D! D im n!(n +1) n 


n>% a, n= (m1) n! i sa (n +1)" -(n+1) n! p 


(după simplificări)= lim | x | Sa <1, deci seria 2 Re converge. 
n=œ În+l e n n’ 
sul n a n+l 
) 2 ao > CU aa Se ŞI Anu POEZII 
Atunci: 


n 2 
E e Pieter a tal DE SIE 


n>œ a 


i n? +1 
- z =2 lim RA 
z n=% n +2n +2 2 n-o n +2n +2 


=2.] = 2 > [ldeci, seria 


—;— „ diverge. 
n=1 n +1 
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3) È — ; cu a, =— ŞI ap = —. Avem 

n=1 N 
lim 1+ = lim —— =1, deci conform acestui criteriu nu putem preciza natura seriei 
n>% a, n> n +1 


l l la a a : z l 
armonice > — (despre care ştim însă că este divergentă, vezi formula (7.4)). 
n=1 N 


IV. Criteriul rădăcinii (Cauchy) 


o0 
Fie seria > a 
n=l 


a 2>0şi lim ya, =A e€ [0,9%]. Atunci, dacă: 
n—>0 


n? 


panler O 
n=l 


i) X>1> Sa, (D) 


n=l 
II)X =1= caz de nedeterminare (nu se poate preciza natura seriei) 


Exemple: 

ie ad 2 
DL Die cu pan >0. Atunci: 

n=l 2 2 

2 n 2 2 
. i i l şi ata E 

A= lim ya, = lime 2 — lim (ln) = — <1 deci, seria > = converge. 

n—09 n—o 2 n—o 2 2 n=1 2N 


n 


2 2 
58 1 n i 
D L = | ; cu a |) > 0. Atunci: 


l . : TOR: 
Wa = I =], deci nu putem preciza cu acest criteriu, 


e mat a LAN > ; iz ORS 
natura seriei >, =F (despre care ştim că este convergentă, ca seria armonică generalizată, 


n=l 


pa cua=2>1). 
n=1 N 
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Observaţii: 


i) 


În unele din exemplele de mai sus am folosit limitele fundamentale: 


limYn=1 şi lim 1/P, (n) =1 (7.22) 


n= 


unde: P, (n)=a,n“ +a, n“ +...+an? +an-+ay polinom de grad k cu 


coeficienți reali şi a, > 0 


ii) 


iii) 


v) 


Se observă că aplicarea acestor criterii, constă doar în calculul unei limite; totuşi dacă 
termenul general al seriei a,, are o formă mai complicată determinarea lui A poate fi 
o problemă deloc uşoară. Se poate apela în aceste cazuri la criteriile de comparaţie 
sau la alte criterii de o formă mai complicată care nu sunt prezentate aici (se poate 
consulta în acest scop, orice curs de analiză matematică care tratează convergența 


seriilor numerice). 


Datorită relaţiei: lim #4 =)= lim yan (7.2.3) 


n-o a, n—oo 
rezultă că dacă un criteriu ne conduce la cazul de excepție A =1, nu se încearcă 
aplicarea celuilalt criteriu, deoarece se va obține aceeaşi concluzie; obişnuit aplicarea 
criteriului raportului sau rădăcinii se face ținând cont de forma termenului general a, 
(astfel încât calculul limitei să fie cât mai simplu). 
Se observă, că există situații — chiar pentru serii foarte simple — în care aceste criterii 
nu pot preciza natura seriilor (à = 1). De obicei în aceste cazuri se aplică criteriile de 


comparație respectiv criteriul general de divergență (C.G.D), sau următorul criteriu 


(a lui Raabe-Duhamel). 


V) Criteriul lui Raabe — Duhamel 


Fie seria >_a, „a, Z0 şi lim | Se e] 


= p € R. Atunci, dacă: 


n=] n—co 


a 


n+l 


i) p>1> $a, (© 
n=l 


i) <11> 5 a, (D) 
n=l 


iii) =] ; caz de nedeterminare (nu se poate preciza natura seriei) 
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Exemple: 
œ 1l 1 i tou mecla ata 

D X — ; u a=- Z0. Am văzut cum criteriul rădăcinii ne conduce la cazul de 
n=] N n 
nedeterminare A =1. Vom aplica criteriul lui Raabe — Duhamel pentru a “scăpa” din 


această situație. 


Avem: 
2 2 2 
p= lim d tim n e4- iin COA e |- 
n= an+ n= n n= n 


. 2n++l 
= lim 


n—>oo n 


= 2 >1, deci seria converge. 


DL ; cu a, = >0. Se verifică imediat că atât criteriul raportului cât şi 
n=1 2n+l 2n+l 
cel al rădăcinii ne conduc la cazul A =1. Dacă aplicăm criteriul lui Raabe — Duhamel 
obținem: 
p= lim n| —1|= lim [zt -1 lim | A - 
n—00 any n—oo 2n +1 n=>œ\| 2n +1 


deci (în acest caz) nici criteriul lui Raabe — Duhamel nu ne “scoate” din cazul de 


nedeterminare. 


Dacă aplicăm în schimb criteriul de comparație la limită, considerând seria armonică : 


œ% œ l a n 1 
b. =5— (D), obținem k = lim — = lim ———— = —. 
2 = Z os n>œb, n-o 2n+l 2 
l prea ; a , . 2 1 
Deci cele două serii au aceeaşi natură, cu alte cuvinte seria >. ER (D). 
n=l n 


Obeservaţii: 


i) Criteriul lui Raabe — Duhamel nu se aplică (în mod obişnuit) decât atunci când 
criteriul raportului sau rădăcinii ne conduc la cazul de nedeterminare A =1. În cele 
mai multe dintre situaţii, acest criteriu ne scoate din cazul de nedeterminare; în caz 
contrar se va încerca aplicarea criteriilor de comparaţie. 

ii) Totuşi, există situaţii unde datorită formei termenului general al seriei aplicarea 
directă a criteriului lui Raabe — Duhamel aduce o simplificare a modului de calcul a 
limitei. 
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7.3 SERII ALTERNATE 


9,9; 
Definiţia 7.3.1. O serie $ u, , U, E R se numeşte serie alternată (cu termeni alternați) 
n=l 


dacă, termenii săi satisfac relațiile: 


uU <0 , (Vn eEN* (7.3.1) 


Observații: 
i) Din (7.3.1) rezultă că termenii seriei alternează ca semn; atunci termenul general u, 
poate fi scris sub una din următoarele forme: 


u, =(— Dra a, > 0 (dacă primul termen al seriei este pozitiv) 


n? 
respectiv: 
u, =(—1)a a, >0 (dacă primul termen al seriei este negativ) 


n 2 
ii) Deci orice serie alternată poate avea una din următoarele două forme: 
es n+l n+l 
S(—D a, =a a, +a; —a +... +(—D ` a, +... 


n=l (7.32) 
S (Da, =—a +a, —as ka... +D’ an +... 
n=l 


Exemple: 


oo 1 1 1 1 1 
a piia PRE e (0 | i pad area Sa 
) 2 ) 3 x (—1) z 


by CP psd ED o 4 CE ii 
n=l 


În studiul convergenței seriilor alternate nu pot fi folosite criteriile de la serii cu termeni 


pozitivi. Vom prezenta aici un singur criteriu specific seriilor cu termeni alternanţi. 
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Criteriul lui Leibniz: 


(9,9) 
Fie seria alternată > (—1)'a, cua, >0. Dacă şirul (a, un satisface condiţiile: 


n=l 
D a, —0 


ii)(a, nen este monoton descrescător (a, > anp, (n € N*) 


E 
atunci seria X` (—l1)'a, converge. 
n=l 


Observații: 
i) Criteriul lui Leibniz oferă numai condiții suficiente de convergență nu şi necesare; 
cu alte cuvinte dacă condițiile nu sunt satisfăcute nu rezultă că seria diverge, ci doar 


faptul că criteriul nu poate fi aplicat. 


ii) Totuşi, dacă lim a, #0, atunci ` (—1)"a, (D), dar datorită aplicării criteriului general de 


n=o0 n=l 
. m . x n 
divergență lim u„ = lim (—1)* ap #0 . 
N—O n= 
iii) In cazul în care a, — 0, dar şirul (a, „ex nu este monoton descrescător criteriul lui Leibniz 


(9.9) 
nu este satisfăcut, şi nu putem afirma nimic despre natura seriei X (—1)"a, . 


n=l 
Exemple: 
~% ny l 1 1 1 ny l l Ma P 
a) 5 GCD .—=1 | H... +(— D) H... (seria armonică alternată) 
= n 2 3 4 n 


l 1 À P i gi a aa n, da e 
Avem şirul a, =—— 0 şi este monoton descrescător, deci conform criteriului lui 
n 


Leibniz seria converge. 


a, n 1 Poi a EEE EOE, 
b GCD’ —— ; cu a, = — — deci prima condiție a criteriului lui Leibniz nu 
n= 3n +1 3n+1 3 
este satisfăcută, acesta neputând fi aplicat. Totuşi, deoarece: 


; ; n ras _ PI E 
lim u, = lim(—1)" —— nu există, rezultă conform C.G.D. că seria diverge. 
= 3n+1 


N—O n> 


„2 i 2 
o9 sn n sn n 

& EG „cu a, = >0 
n=l vn NN 


i 1 A T i . 
Deoarece 0<sin?n <1 > 0< a, <->. Trecând la limită (teorema cleştelui) obținem 


Jn 


lim a, = 0, deci prima condiție din criteriul lui Leibniz este satisfăcută. Se poate arăta relativ 


n> 


uşor că şirul (a, )aey nu este monoton, deci condiția a doua nu este satisfăcută şi criteriul lui 


Leibniz nu poate fi aplicat. 
Pentru a putea totuşi preciza natura acestei serii, se pot aplica criterii de la serii cu 


termeni oarecare (neprezentate aici) sau folosi noţiunile din următorul subcapitol. 


7.4 SERII ABSOLUT CONVERGENTE. 
SERII SEMI — CONVERGENTE. 


(9,9) 
Evident, în cazul seriilor cu termeni oarecare e aa’ dn € r), nu se pot aplica în studiul 


n=l 
naturii acestora, criteriile formulate pentru seriile cu termeni pozitivi sau pentru seriile cu 
termeni altenanți, enunțate în secțiunile precedente 7.2 respectiv 7.3. 

Totuşi, în secțiunea 7.1 au fost enunțate două rezultate care pot fi aplicate tuturor seriilor 
numerice, şi anume: Criteriul general de convergență a lui Cauchy, respectiv Criteriul general de 
divergență. Aplicarea lor, în special pentru stabilirea convergenţei seriilor, este aşa cum am mai 
spus, de cele mai multe ori extrem de dificilă. De aceea, matematicienii au căutat şi au stabilit, 
noi criterii, mai simple de aplicat în practică, cum ar fi: Criteriul lui Abel, Criteriul lui Dirichlet, 
etc. 

Chiar şi aşa, aplicarea lor necesită abilităţi matematice mai deosebite, precum şi o 
oarecare practică în lucrul cu serii. Deoarece lucrarea de față se adresează în mod special 
studenţilor economişti, s-a considerat utilă renunţarea prezentării acestor noi criterii; vom arăta 
totuşi că, în anumite situaţii, studiul convergenţei / divergenţei seriilor cu termeni oarecare se 


reduce la studiul convergenţei / divergenţei seriilor cu termeni pozitivi. 
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(9,9) 
Definiţia 7.4.1. Fie seria cu termeni oarecare Xa, ,, a, ER. Spunem că: 
n=l 


OO O 
a) Seria J a, ,, este absolut convergentă (A.C.), dacă seria moduleleor > lanl este 
n=l n=l 


convergentă (ca serie cu termeni pozitivi). 


b) Seria 5 a, , este semi - convergentă (S.C.), dacă ` a, (C)dar, 5 lanl (D). 
n=l n=l n=l 


(9,9) 
Teorema 7.4.1. Dacă seria > a,, a, ER este absolut convergentă, atunci şi seria 
n=l 


Sa, este convergentă. 
n=l 
Demonstrație. Deoarece, seria Xa, (A.C.) > Xa] (©. Aplicând Teorema 7.1.3. 
n=l n=l 


seriei modulelor, obținem: 


(v)e >0, (n. EN aî. 


anu] ans] = antan <E n n, 
şi (V)k E N* (5) 

Din proprietățile funcției modul, rezultă că: 
ngeia Elaa e a 


(e9) 
Din (%) +(*%) rezultă că Va, (C), conform criteriului general de convergență a lui 
n=l 


Cauchy. 
q.e.d. 


Observaţie: 


Reciproca teoremei de mai sus nu este adevărată; cu alte cuvinte, convergenţa seriei simple 


nu implică convergenţa seriei modulelor Sa (© $ San O A 
n=l n=l 
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; e e 1 ; . E fi 
Exemplu: Fie seria Ya, = $X (—1)°—. Evident este o serie alternată şi converge deoarece 


n=l 


œ œ | 
satisface criteriul lui Leibniz. Totuşi, seria modulelor, X lanl = > —, diverge ca serie armonică 
n 


n=l n=l 


simplă. 


A . . 9 . . . . . 
1. Întrucât seria valorilor absolute, X, lan „ este o serie cu termeni pozitivi, pentru studiul 


n=l 


naturii acesteia se vor folosi criteriile enunțate în secțiunea 4.3 


A . E y v . . v C9 
2. In cazul în care seria modulelor X lanl (D), am văzut că despre seria simplă >. a, , nu se 


n=l n=l 


poate afirma nimic. Totuşi, se demonstrează că în cazul în care divergența seriei 


(9,9) 
modulelor > a, | a fost stabilită cu Criteriul raportului sau cu Criteriul rădăcinii, atunci 
n=l 


O 
şi seria (simplă) J, a, diverge. 
n=l 


OO 
3. Pentru seriile cu termeni pozitivi 5 ap an >0], noțiunile de convergență şi 


n=l 
7 , (9.9) (9,9) 
convergenţă absolută coincid (deoarece 3-a, => lanl ). 
n=l n=l 
4. În anumite situații, este utilă şi următoarea formulare echivalentă a Teoremei 4.7. “Dacă 


(e0) O 
seria X a, a, E R diverge, atunci şi seria X |a,| diverge”. 


n=l n=l 


A . A . .. .. . e 
In concluzie, în studiul naturii seriilor cu termeni oarecare, Ya, a, ER, se poate 
n=l 


aplica următorul: 
Mod de lucru: 
. . . . . oo . 
i) Se ataşează seriei cu termeni oarecare, seria modulelor > lanl (care este serie cu 
n=l 


termeni pozitivi); 
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(2,9) 
il) Se studiază natura seriei >. la, | folosind criteriile enunțate pentru serii cu termeni 
n=l 
pozitivi; 
iii) Dacă: 


a) > la (O => > a, (©) 


(9,9) 
b) Xa] (D) (divergenţa fiind stabilită cu Criteriul raportului sau Criteriul 
n=l 
. .. Sat, 
rădăcinii => > a, (D) 
n=l 
(9.9) 
c) 5 la,| (D) (divergenţa fiind stabilită cu alte criterii) = despre natura seriei 
n=l 
| A . . . .. 
X a, nu se poate spune nimic. În acest caz trebuie folosite alte criterii (Abel, 
n=l 
Dirichlet etc.). 
Exemple: 


n a TI 
-sin” a 
cu a€ 0,27 parametru real. 


DA = 
n=l Jn i 
Evident seria este o serie cu termeni oarecare (datorită lui sin'a) şi nu cu 


termeni alternanți. 


[e0] OO Pa i sin” a] 
Ataşăm seria © |a,|= 5 R şi aplicăm criteriul raportului. Atunci, obținem: 
n 


n=l n=l 
n+l . n+l / 
; a] __2 sin a] vn , n j. 
A= lim = lim ———= .—— = lim 2 —— -jsinal 
noo Ap n=co n+l 2" Isin” al n=œ \n+1 
şi deci: 
A = 2lsinal 


Atunci, pentru: 


i) A=2ļsina|<1, seria X |a,| (©) = $ a, (©) (conf.T.7.1) 
n=l n=l 
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ii) X2lsina| >1, seria 5 a.| ©) => 3 a, (D) (conf.Obs.3) 
n=l n=l 


iii) XE 2jsina|=1, nu putem preciza natura seriei Dlan | şi deci a seriei inițiale Sa, : 


n=l 
Rezovând inecuațiile trigonometrice de mai sus, avem: 


g U Aan U IS 
6 3 6 6 


2" sin” a 


= (a) Dă (A.C) 


i) pentru a € 


m 2m) (7m lln sin” a 
ii entru a E€|—,— |U D) —— = (D 
D p z aE j Sp Eta o 
PS n 2n 78 lln A yi oœ . 2 
iii) pentru «elz, e |> nu putem preciza natura seriilor > |a,| şi X a, . 
n=l n=l 


Dar, în cazul iii) cele două serii (simplă şi a modulelor) au expresiile: 


a) Sa =% (- ID NA (C), ca serie alternantă care satisface criteriul lui Leibniz 
=l 


; 1 i 
(şirul —= — 0, monoton descrescător). 


n 
b) Sk, => 


giga SER A] 
(D), (ca serie armonică > —— cu a = 5 <1) 


1 
ari bă 
=INN n=l ar n=1 N 
n . 2 2" sin” a 
Deci, în acest caz, seria 5 (— 1)" =— = 
n=l Jn 


2) s] E] , cu aeR\ |i] parametrul real. 


n 


este (S.C.) 


Ataşăm seria modulelor: = la,|= 5f a ii i == şi aplicăm criteriul rădăcinii: 
ni N — ză 
; A nir E la n+l |l-a l-a 
`= lim pla = lim n . = lim —-. = 
n= Pnl n—09 | n = n>% n = E 


Atunci, pentru: 


D) A= E ASE la,| (C) „deci şi seria > a, (C) (conf. T.4.7) 
— ză n=l n=l 

i) A= — >13 Sla, (D), deci şi seria A (D) (conf. Obs.3) 
— ză n=l n=l 
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1l—a 
1—2a 


iii) >=] 


=] = nu putem preciza natura seriilor > a, | şi > an 
n=l n=l 


Atunci, rezolvând inecuațiile de mai sus, concluzionăm că: 


i) pentru XE —œ,0 U P = 5 ENA este (A.C.); 
3 anal n 1—2a 


A 2) [1] sny (1-a Y 
ii) pentru vefot s=] f = este (D); 


n=l n = ză 


iii) pentru XE [E 


| => nu putem preciza natura seriei. 

Dar, pentru: 

co œ (n+1 Ý 

a) A=0=> a, Z}. |— | „care (D) conform C.G.D. 
ıl n 


n=l n= 


Într-adevăr 


lim a, = lim =] = lim h+] ze 
n—oo n= n n—>oo n 
n+1 j 

n > 


Dais Sa Sen | 
3 n=l n=l 


care (D) conform C.G.D., deoarece a, 0. 
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US8 Serii de puteri 


Seriile de puteri sunt un caz particular de sume infinite. Ele sunt utilizate pentru 
aproximarea funcțiilor transcendente. 

Capitolul de faţă abordează modul în care poate fi studiată o serie de puteri: 
determinarea razei de convergență, studiul convergenței în extremitățile intervalului de bază, 
determinarea sumei, proprietățile sumei. 

Studentul trebuie să rețină metoda de aflare a razei de convergență și să realizeze 
că în capetele intervalului de bază obţinem o serie numerică de tipul celor studiate în US7. 
Pentru partea teoretică două ore sunt suficiente. Pentru aplicații sunt necesare măcar patru ore: 
trei pentru mulțimea de convergenţă şi una pentru studiul in capetele intervalului. 


Utile din lista bibliografică se dovedesc a fi lucrarile de la poziţiile [3], [6] și [9]. 


x 
x * 

Aşa cum s-a observat din secţiunile precedente, noţiunea de “serie numerică” apare ca o 
generalizare (extindere) firească a noţiunii de “şir numeric”. 

În mod natural, aceste două noţiuni au căpătat în Analiza matematică noi extinderi, din ce 
în ce mai complexe, lucru de altfel deseori întâlnit nu numai în matematică ci şi în multe alte 
domenii ale ştiinţei. 

Astfel, pornind de la noţiunea de şir numeric: 

(a ne doo apa: da ++- (*) 
se introduce noţiunea de “şir de funcţii”: 

f(x) Ei hE) fî (3), (**) 
unde funcțiile f, (x) sunt definite pe acelaşi domeniu DC R. 


Exemple: 


fi :R—R 


a) f(x) nene? CU termenul general: | e atoian 


adică, şirul de funcții are forma: 


; 4 6s In s 
xsinx, 2x sin2x, 3x° sin3x,... nx" sin nx,... 
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f, :(0,00)— R 


b) fi) ao unde termenul general este: | fe) a E e 
deci, şirul de funcţii are forma: 


2Inx, 5ln2x, 101n3x,...,(n2 +Dln(nx)... 


Ținând seama că noţiunea de “funcţie “ este mult mai complexă decât cea de “număr”, 


este evident că noţiuni precum convergență sau divergență au căpătat pentru şirurile de funcții 
valenţe noi. Astfel, un şir de funcții f, (x) ney Poate fi: simplu (punctual) convergent, 
uniform convergent, slab convergent, tare convergent etc. 

Este natural ca în caz de convergentă, şirul de funcţii să aibă limita o funcție f(x) şi nu 


un număr. De asemenea, ţinând cont că termenii şirului sunt funcţii, care pot avea diverse 
proprietăţi (funcţie continue, derivabile, integrabile etc.), este natural să cercetăm în ce condiţii 


= 


acestea se “transferă” la funcţia limită f(x). 


In mod cu totul analog noţiunea de “serie numerică”: 


(2,9) 
> a, =a ta ta, boi Han tH (FEE) 


n 
n=0 


capătă o primă generalizare şi anume cea de “serie de funcţii”: 


ăi MAAE AES AOE ES ACE A em) 


unde termenii seriei sunt funcții definite pe acelaşi domeniu DCR. 


Exemple: 


a) 2 (n? + Darcte(x" +1) = arctgl + 2arctg(x +1) +5arctg(x? +1) +...unde termenul 
n=0 
fi :R—>R 


general al seriei este: | £ (x) = (n? + Daretg(x" +1) 


n! Il! 2! 3! 


b) = z t= 
nı Î+xX" 1l1+x 1+x 1+x 
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£:01—R 


unde termenul general are forma: Anm n! 
> 1+x°” 


Analog, în caz de convergență (simplă, uniformă, etc.) seria va converge la o “funcție 
sumă” s(x). Din nou este interesant şi important de studiat, în ce condiții proprietățile termenilor 
seriei de funcţii (integrabilitate, derivabilitate, continuitate etc.) se transferă asupra funcţiei sumă 
s(x). 

Un caz particular de serii de funcții, dar foarte important şi des intâlnit în aplicațiile 
economice, este cel al seriilor de puteri, unde termenii seriei sunt funcții putere. Astfel, avem: 


Definiția 8.1: Numim serie de puteri centrată în x, , o serie de funcții de forma: 
X a (X-X0) = ag ta(X-X0) FoR +... +a, (xx) +... (8.1) 
n=0 


unde (a, )ney este un şir numeric (Numerele a, ER se numesc coeficienții seriei de puteri). 
Observaţii: 
i) un caz particular, dar foarte des întâlnit în aplicaţii, este cel al seriilor de puteri 


centrate în 0 (pe scurt, serii de puteri) deci când x = 0. Vom obține atunci seria: 


Sasa aaa? En Ta Eau, (8.2) 
n=0 
il) Deoarece prin substituţia: 
YA (8.3) 


o serie de puteri centrată în x, (Şi variabila x) se reduce întotdeauna la o serie de puteri centrată 


în O (şi variabila y), adică avem: 


5 a(x x) = 2 a y” (8.4) 


ne vom ocupa numai în continuare numai de seriile de puteri centrate în 0, deci de forma (8.2). 


(9,9) 
Definiţia 8.2: Spunem că, seria de puteri Y a,x” este: 
n=0 


OO 
a) convergentă (punctual) în x =X, ER, dacă seria numerică ataşată, Y a,Xp este 
n=0 


convergentă; 
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b) convergentă (punctual) pe domeniul D) C R , dacă converge în orice punct x € Dg 


(D, se numeşte mulțime de convergență a seriei). 


(9,9) 
Teorema 8.1. Mulțimea de convergenţă a unei serii de puteri, y_a,x", este 
n=0 


întotdeauna nevidă. 


(0,9) 
Demonstraţie. Fie seria de puteri >» a,x” cu (a, en un şir de numere reale oarecare şi 
n=0 


D, mulțimea sa de convergenţă. 


(9,9) 
Pentru x =0 seria > a,x” ay, deci converge şi are suma S= ay. Rezultă că 0eD,, 


n=0 


deci D #0. q.e.d. 


Teorema 8.2. (lui Abel) 


(9,9) 
Fie seria de puteri, Y_a,x", a, ER. Atunci, (Dr € 0,00 (“r” se numeşte rază de 


n=0 
convergenţă), astfel încât: 
i) pentru x € (—r,r) = z a,x” converge (C) 
= (8.5) 


ii) pentru x € (—00,—r)U(r, +00) = a,X" diverge (D) 
n=0 


Observaţii: 


i) Teorema lui Abel nu precizează natura seriei în capetele intervalului de convergență: 


(-r, r). In exemplele concrete, se studiază separat natura seriei în cele două capete, 


adică seriile numerice: 


œ n. X i 
X Ar ŞI X a (r) . 
n=0 n=0 
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OO 
ii) Dacă raza de convergenţă r =0, atunci evident seria de puteri > ax" converge 


n=0 
numai în x = 0 şi diverge în rest (pentru x € R*). 
iii) Dacă raza de convergență r=-+oo, atunci mulțimea de convergenţă a seriei 
Do = R, deci seria nu diverge pentru nici o valoare reală. 


iv) Aşadar, folosind teorema lui Abel, dacă aflăm raza de convergenţă r, ştim mulțimile 


œo 
pe care seria Y_a,x" converge sau diverge (evident cu excepția valorilor x = r ). 
r=0 


Următoarele două teoreme ne oferă două formule pentru aflarea razei de convergență 


r. 


Teoremă 8.3 (Teorema raportului). 


oo a 
Fie seria de puteri, Y_a,x" cu a, ER şi lim | = = pe 0,œ . Atunci raza de 
n=0 n— oo an 
convergență r , este dată de relațiile: 

l 
— ;pe(0,00) 
P 
0  ;p=+ 


Teorema 8.4 (Teorema rădăcinii) 


(9,9) 
Fie seria de puteri, X a,x” cu a, ER, şi lim Ya, | =p€ 0,00 . Atunci, raza de 
n=0 n= 


convergență r, este dată de relațiile: 
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r=i+oo ;p=0 (8.7) 


Observaţie: Din rezultatele enunțate anterior, se observă că în determinarea mulțimii 


(intervalului) de convergenţă respectiv a mulţimii (reuniune de intervale) de divergență pentru 


(2,9) 
seria de puteri > ax", trebuie aplicat următorul: 


n=0 


Algoritm de lucru: 


1) Se calculează valoarea “p” aplicând una din formulele: 


azi sau p= lim ja, | 
| n—00 


an 


p= lim 
2) Să determină raza de convergență “r” cu relațiile (8.6) sau (8.7) (care de fapt coincid); 
3) Se scriu intervalele de convergență / divergență ale seriei 3 aX”; 
n=0 


4) Se studiază natura seriei în cele două capete x = +r. 


Las . w A 09 . 
Observație: In cazul în care seria este centrată în xy | > a, =] se reduce la seria 
n=0 


centrată în 0 | D any] folosind substituția (8.3) şi apoi se aplică algoritm de lucru, obținându-se 
n=0 


mulțimea de (C)/(D) pentru variabila y. Folosind aceeaşi substituție se scriu mulțimile de (C)/(D) 


în variabila x. 


Exemple: 


a) D CŽ; Avem: a, =(-1) a deci je, 
n= n n n 


Folosind teorema rădăcinii, obținem: 
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p= tim dar = lim e = fim ae =i 


(am folosit aici limita fundamentală: lim a/n = 1) 


N—O 


1 | 
deci conform (8.7), avem: r = — z deci r =1. Atunci, conform Teoremei 8.2. (Abel) 
p 
obținem: 


i) pentru x €(l, SFE Îj SO 


n=l 


ii) pentru x € (—œ,—1)U (1, o0) = 5 (—D a (D) 
n=l n 


z œ RE na : 
iii) pentru x = +1 > 5 (—1)" — nu-i putem preciza natura. 
n 


n=l 


Dar, pentru: 
. 29 n 9 n 1 . $ . $: . . .. 
x =1 obținem: a,x =} (—1) — (C) (conform Criteriului lui Leibniz de la serii 
n=l n=l n 


alternante) 


iar, pentru: 


, œ i 3 1 1 
x=-—l1, obținem: ` a,x" =>) — (D)(ca serie armonică S — cua=l) 
n=l n= 2 n n= n“ 


i si (C), dacă xe(—1,1] 
În concluzie, 5 (—1)" A <7 
n=l n (D), dacă x E€(—œ,—1]U (1, +020) 


o R 1 f 
2) SC; avem a, = (1Y .— , deci |a,| e 
n=0 n! n! n! 


Folosind teorema raportului, obținem: 


| 
n: —0 


Cao fa] (n+)! noni 


şi deci conform (4.6.6) avem r = +oo. Deci, seria D (=i L (C) pentru (YV) x ER. 
n=0 n: 
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œ 1 ; i z 
DE = (x +1)", este o serie de puteri centrată în xy =—1. 
n=0 N: 
oa œ% RS 1 
Notăm y =x +1 şi obţinem X Ha yua = aar = la] 
n=0 n:2" n:2" 
Atunci: 
MOET I R 2.49 1 
= lim Part = tim e ca dia a e 


= . == 1m == 
no [a] n=% (nr D221 n+1 2a>œn?+2n+1 2 


} 7 1 : 
Atunci raza de convergență r = — = 2, şi avem: 
p 


| oo 1 i 
d pentru ye(-2,2)>5 Ey (0) 
n=1 n2 
A eS 1 
ii) pentru y E(—œ,—2)U(,œ)> 5 =k (D) 
n—A N 


iii) pentru y = +2 =? (nu putem preciza natura seriei) 


Revenind la variabila x(x = y —1) , avem: 


n+l Hp (©) 
n2” 


i) pentru x € (—3,1) = 5: 
n=l 


HI + D) 


ii) pentru x € (—~œ,—3) U (l, œo) => 2 z an 
n=l N: 


iii) pentru x = —3, x =1 =>? (nu putem preciza natura seriei) 


Dar, pentru: 


x=—3, obținem D a 8+ = 5 Tan (D) (conform C.G.D, deoarece 
n=l n=l n 


termenul general (—1)" iul. — 0) 
n 


iar pentru: 
. SS, nano & n+l 
x = 1, obținem > a (%3 +D" => —— (D) (conform C.G.D, deoarece termenul 
n=l n=l n 


general sa 


—1z0). 


114 


US9 Funcţii de mai multe variabile 


Funcțiile de mai multe variabile constituie generalizarea firească a funcțiilor de o 
variabilă reală. Studiul lor se justifică prin prezenţa acestora în majoritatea modelelor economice. 
În prezenta unitate de studiu se extind treptat noţiuni familiare studenților: domeniu de definiţie, 
grafic, limită, derivată, diferenţială, punct de extrem. Am presupus cunoscute doar noţiunile 
fundamentale de analiză matematică şi am prezentat aproape totul în material. Credem că astfel 
ajutăm studenţii proveniţi din licee care au în programă mai puţină matematică să îşi însușească 
şi să aplice aceste metode. Noţiunea centrală este cea de punct de extrem cu care capitolul se 
încheie. 

Studenţii trebuie în principiu să înveţe să privească funcțiile de mai multe 
variabile ca o modalitate mai riguroasă de a descrie realitatea. Ei trebuie să înveţe să calculeze 
derivatele parţiale (de diferite ordine) ale unei astfel de funcții, să determine diferențialele (de 
ordinul I şi II) şi să le aprecieze rolul în studiul funcţiilor, să determine şi să caracterizeze 
punctele de extrem. 

Este un capitol amplu, cu multe noțiuni noi şi cu nenumărate aplicaţii. Pentru 
acest capitol zece ore pot fi suficiente: două ore pentru partea teoretică, două ore pentru 
deprinderea calculului derivatelor parțiale de ordinul I şi a diferenţialei, două ore pentru 
derivatele parțiale de ordin superior şi a diferenţialei a 2-a, patru ore pentru aplicaţiile la punctele 
de extrem. 

Orice manual universitar de analiză matematică are aceste noțiuni incluse. Pentru 


doritori trimiterile bibliografice [1], [3] şi [6] pot fi de asemenea folositoare. 


x 


x x 


Fenomenele vieții reale sunt descrise cel mai adesea de funcții care depind nu de o singură 
variabilă ci de mai multe variabile independente. Acest lucru este valabil şi în viața economică. 
În acest capitol se introduc elementele de bază ale teoriei acestui tip de funcţii. Pentru a 
depăşi dificultățile extinderii conceptelor fundamentale (continuitate, limită, derivabilitate etc.) 
de la funcţii de o variabilă la o funcţie de mai multe variabile am abordat numai cazul funcțiilor 


de două variabile, dezvoltarea ulterioară făcându-se în aceeaşi manieră. 
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9.1. ŞIRURI DE PUNCTE ÎN R’ 


Să considerăm spaţiul vectorial R" şi să considerăm un vector oarecare 


X = (Xp X35.::X,) al său. Pentru a simplifica scrierea vom renunţa la simbolul de transpunere şi 
în cele ce urmează vom scrie pur şi simplu x = (X,, X3s--:5X). Acestui element al spaţiului R” 


îi asociem un punct P care are coordonatele x,,x,...,X, , tot aşa cum unui vector x =(xX,,X5) 
din R° îi asociem un punct P cu coordonatele x,,x, într-un sistem de axe ortogonale. 


Asociem unei perechi de astfel de puncte P, Po P(x „x... „x, )Py(x0,x9,...,X0) un 


număr nenegativ numit distanţă şi definit prin: 


02 032 02 
IERI (a, 70) Hr 0 a a (911.1) 
Evident, distanţa este nulă dacă şi numai dacă P şi Po coincid. Convenim că însuşi 
vectorul x să se numească punct din R”. În acest mod, distanţa d se asociază unei perechi de 


vectori din R”. 


Definiţia 9.1.1. Se numeşte şir de puncte din R” funcţia care asociază oricărui 
număr natural k un punct unic x“ din R’. 
k k) 


k k 
k—x = (X Xaa Xa 


Şirurile (i hn (hea (0 ei se numesc şiruri de coordonate. 


Definiţia 9.1.2. Şirul Jan converge la punctul x’ dacă: 


lim d(x“, x°) =0 (9.1.2) 


În spiritul definiției cunoscute a limitei unui şir de numere reale aceasta înseamnă: 
(Y) e >0,3 k(£) aî.(Y)k >k(2) d(3*, x°) <€ (9.1.3) 
Folosind (9.1.1) deducem că: 


Ne ee AA e e E AEE p S D a 
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pentru k > k(e) adică 


xE —x9|<e (Y) k>k(e), (Y)i=1,n (9.1.4) 


sântem conduşi astfel la: 


Teorema 9.1.1. Şirul x“ este convergent la x" dacă şirurile de coordonate x: sunt 


convergente la x? pentru orice i=1,n. 


Aan E . Şirurile de coordonate sunt 
k 3k 


> 


Exemplu. Considerăm în R° şirul x“ =l 


i ar x$ = ss şi pentru k — oo au limitele 0 şi respectiv K, Aceasta înseamnă că 
k 3 3 


x“ are limita x? = l5) 


9.2. FUNCŢII DE DOUĂ VARIABILE 


În cele ce urmează vom considera spaţiul R? ale cărui puncte x au doar două coordonate 
x = (x,,X). Pentru ca scrierea să fie mai simplă vom înlocui perechea de coordonate (x,,x) ) 
prin perechea (x,y). 

Fie D o mulţime de puncte (x,y) din R? şi I o mulţime din R. 

Definiţia 9.2.1. Se numeşte funcţie de două variabile tripleta formată din D.I şi o 


lege de corespondenţă care asociază oricărui punct (x,y) € D, un singur element z € 1. 


Exprimăm această corespondenţă prin z = f(x,y). 


Observaţie. Când D este din R” iar ICR în mod similar se defineşte o funcţie de n 
variabile. 
Dacă luăm un sistem cartezian de coordonate Oxyz şi reprezentăm mulţimea punctelor 


P(x,y,z) unde z=f(x,y) obținem imaginea geometrică a funcţiei sau graficul său. De data 


aceasta graficul e o porţiune dintr-o suprafaţă. 
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Z 


Fig. 5.1. 


Exemplu: Să se reprezinte graficul funcţiei z = x? +y7.. 


Soluţie. Ridicăm la pătrat egalitatea şi obţinem z” = x? +y” . Dacă facem aici x = 0 


(intersectăm graficul cu planul yOz obţinem z = +y adică două drepte în acest plan. În 


mod similar intersectăm graficul cu planul xOz şi obţinem dreptele x = +z. 
Graficul este un con infinit cu vârful în originea sistemului de coordonate. 


Fig. 9.2 


118 


9.3. LIMITĂ ŞI CONTINUITATE 


Să considerăm o funcţie f:DCR?—R şi (x9,yg) un punct de acumulare al lui D. 


Aceasta înseamnă că există cel puţin un şir de puncte (x,.Yi)en din D care converge la 


(0 Yo). 


Definiţia 9.3.1. Funcția f(x,y) are în (Xọ,Yyọ) limita L, dacă oricare ar fi şirul 


(Xk Yk)ken €D, (Zk Yk)ken — (Xo Yo); (Sk Yk) 7 (Xo: Y0) rezultă că f(X k Yk) >l. 


Vom scrie: lim f(x,y)=£ sau lim f(x,y)=£. 
(x.y)—>(X0-Y0) X—X0 
Y—Yo 

Această limită se numeşte limită globală deoarece x şi y tind simultan şi independent 

către Xp respectiv yo. Dacă xọ şi y tind succesiv la xọ respectiv yo adică dacă 


L= lim (lim f(x,y)), (5 = lim ( lim f(x,y)),, se obțin limitele iterate sau parţiale. Are loc: 
Y—Y0 X—X0 


X—X0 Y—Y0 
Teorema 9.3.1. Dacă limita globală există atunci şi limitele parţiale există şi are loc 
egalitatea: 
A A =k (9.3.1) 
Reciproca nu este adevărată. 
Demonstraţie. Pentru a justifica rezultatul este suficient să considerăm mai întâi şiruri 


generale (x, Yi en —(Xo>Yo) şi apoi să le particularizăm în şiruri de formă particulară 
(Sko Yo en ~ (X05Y0) ŞI (Xo> Yr en ~ (X0» Yo). Ceea ce e adevărat pentru general se păstrează 


şi pentru particular. 
Faptul că reciproca e falsă se poate proba prin următorul exemplu. 


Să se calculeze limitele iterate ale funcției f(x,y) = e 
Asy 


în punctul (0,0). 


A = imi Sim tim 
x—0| y—0 x+y x—0 x x—0 


l, = lim (iim im = lim(—1) =] 
y—0|x—0 x +y y—0 y y—0 


Cele două limite nefiind egale nu există limita globală a funcției în origine. 
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Observaţii: 

19. Dacă (, = l, nu rezultă că limita globală / există ci doar că dacă aceasta există este 
egală cu valoarea lor comună. 

2%. S-ar putea ca în punctul (xo,yo) funcţia f să nu fie definită adică să nu se poată calcula 
f (x0,yo) cum este în exemplul precedent. 


Definiţia 9.3.2. Funcţia f (x.y) este continuă în (x0,y0) dacă: 


lim f(x,y) = f(0,Y0) 
X—X9 


Y—Yo 
Ca şi în cazul limitei globale acest tip de continuitate se numeşte continuitatea 
globală. Se poate şi aici pune problema continuității în raport cu fiecare argument când 
celălalt este fixat, adică continuitatea parţială. 


Exemplu. Să se studieze continuitatea globală şi cea parţială pentru funcţia: 


f(x,y) = x? 
0 (x,y) = (0,0) 


Xy 
DI (x,y) = (0,0), 
+y 


Calculăm limitele parțiale: 
H= lim f(x,0) = lim 0 =0 
lL, = Ie A 0 =0 
Se observă că /, =f(0,0)=/, =0, deci funcția este continuă în raport cu fiecare 
argument. 
Limita globală nu există pentru că luând şiruri (Xk:Yk)ken (0.0) în care 
Yk = MĂ, mz0 obținem: 


X, MX m 
PN ital < k — 
fak Yk) = 2 Ja je 2? 
x; MAX l+m 


deci f(x: Y) = a care depinde de m. Aşadar continuitatea globală nu există. 
+ 


2 
m 
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9.4. DERIVATE PARȚIALE ȘI DIFERENȚIALE 


9.4.1. DERIVATE PARȚIALE DE ORDINUL ÎNTÂI 


Să considerăm funcția f: D—R şi (xo, yo) un punct din interiorul mulțimii D. 


Considerăm de asemenea un punct (x, y) din D diferit de (xo, yo). 


Definiţia 9.4.1.1. Se numeşte derivată parțială a funcției f în raport cu variabila x în 
punctul (xo, yo). 
na f(x, yo) —f(x0- Y0) (9.4.1) 


X—X0 X— Xo 


dacă limita există. Dacă aceasta este şi finită, spunem că funcția este derivabilă parțial în raport 


cu x în punctul considerat. 
ss . . . Of 7 A . . . 
Notăm valoarea limitei cu zi (0) Yo) sau f; (x9Y0). In mod similar se introduce 
x 
derivata parțială şi derivabilitatea parţială în raport cu y în (xo, yo): 


A Yo) = f(xo Yo) = m Aie T n 
Observaţie. Derivata parțială nu aduce nici o noutate față de cazul funcțiilor de o singură 
variabilă. Practic, când se derivează parţial o funcţie de n variabile în raport cu una, toate 
celelelte variabile se consideră constante (parametri). 
Dacă f (x, y) este derivabilă parțial în raport cu x sau/şi y în toate punctele unei 
submulțimi D; a a lui D, spunem că f este derivabilă parțial pe D.. 
În acest caz pentru orice punct (X0>Y9) ED; stabilim  corespondenţele 
(Xo Yo) — Z ayo) (Xo Y0) > Soyo) care definesc două funcții noi numite derivatele 


parţiale ale lui f în raport cu x respectiv y în D.. Ele se vor nota: 


SL y) sau f/ (x,y) respectiv pula y) sau f(x,y). 
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Observaţie. Calculul derivatelor parțiale CL), e y)se face folosind regulile 
x 


obişnuite de calcul pentru derivata sumei, produsului, câtului, funcțiilor compuse etc. 


Exemplu. Să se calculeze CAE y), ol x y) dacă f(x, y) = xy e”. 
Ox Oy 
x „Of vtr 
Soluţii. Calculăm ai y) considerându-l pe y parametru constant. 
x 
of = I Xy Xy; / xy xy 2\ Xy 
a =(xy)e ” HYE”), = ye” +xy(ye ”) =(y +xy")e 


Calculăm E (x,y) considerându-l pe x parametrul constant: 


te Y) EYE HEYNE), = xe Haye) = (x + yo” 


Observații : 


1°. În exemplul de mai sus, simetria în raport cu variabilele x şi y permitea calculul celei 


de-a doua derivate parțiale schimbând între ele variabilele x cu y în a y). 
x 


20. Toate consideraţiile anterioare se pot extinde de la funcții de două variabile la funcţii 


de n variabile fără nici un fel de dificultate. 


9.4.2. DERIVATE PARȚIALE DE ORDIN SUPERIOR 


Să considerăm din nou f:DCR? — R, derivabilă parţial în raport cu x şi y pe D. Cele 
două derivate parţiale fý (x,y) = 8, (x,y) şi f(x,y) =82(x,y) pot fi la rândul lor derivabile pe 


o submulțime D; a lui D. Apare natural să considerăm derivatele lor parțiale drept derivate 


parţiale de ordinul al doilea ale funcției f: 
og O (of of j 
——— a = —— | — P = —— 2 = fix A 5 
3 (x,y) 3 k (x »)| P (x,y) (x,y) 


A T E RR 
ðy (x,y) = 2an] ðyðx (x,y) = fy (3, y), 
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of 
OxOy 


Dea aya DE 
ðy (x, y) a əx [a (x, v) 


(x,y) =f (X,Y); 


ô ð (əf af ; 
E y) = Ey = (Ze J = op y) z fay, y) - 


Pentru funcţia de două variabile f(x,y) obţinem astfel patru derivate de ordinul al doilea. 
Dintre acestea se remarcă f(x,y) şi fy (x,y) numite derivate parţiale mixte de ordinul al 
doilea. 


Natural se pune întrebarea dacă între cele două derivate există legături. Exemplul următor 


sugerează că răspunsul este afirmativ. 


Exemplu. Să se determine derivatele parţiale de ordinul al doilea ale funcţiei 
f:DCR—R, f(x,y) =x°y+2xy +1. 


Soluție. Obținem succesiv: 


of of _ 2 
— =2xy +2y ; — =x +2x; 
xy + 2y zx 


ôx 

2 2 

Zi (E): A E RE, 
ôx? Exo &xây âxloy 


of 2 (2) of (2) 
= —| —|=2x+2 ; == |=0 
ôyôx ðy\ əx ôy” ôy \ əy 


Se observă că cele două derivate parțiale mixte f,,fy, sunt egale. Nu orice funcție f 


(x,y) are această proprietate. Fără demonstrație enunțăm următoarea teoremă: 


Teorema 9.4.2 (Criteriul lui Schwarz) 


Dacă în domeniul D, c D funcţia f (x.y) satisface condiţiile: 


a) are derivatele 5 ale ` 
0X 
2 8? 
b) derivatele mixte de ordinul al doilea ; sunt continue global, atunci 
Ox0y OyOx 


derivatele parţiale mixte de ordinul al doilea sunt egale. 
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Observaţie. Funcţia f (x,y) din exemplul anterior satisface condiţiile criteriului lui 
Schwarz. 


Fiecare derivată parțială de ordinul al doilea poate fi la rândul său derivată parțial în 


raport cu cele două variabile. Se obţin astfel derivate parţiale de ordinul al treilea şi procesul 
poate continua. 


9.4.3. DERIVATELE PARȚIALE ALE 
FUNCŢIILOR COMPUSE 


Să luăm funcţia f :D c R? — R şi să presupunem că la rândul lor, argumentele x şi y ale 
lui f sunt funcţii de variabilele u,v,w: 


f(x, y) = f(x (u, v, w), y(u, v, w)) = ọ (u, v, w) (9.4.2) 
Ne interesează derivatele parțiale ale funcției ọ în raport cu cele trei variabile u,v,w. Un 


calcul nu foarte complicat care foloseşte definiția derivatei parțiale conduce la formula de calcul 
a celor trei derivate parțiale: 


ôp_ ôf ôx ôf ôy 
Qu ôx ôu ôy ôu’ 


Op _ of Ox of Oy 
ðv ôx ðv Oy 


Op ðf Ox əf ðy 
ðw Ox ðw ðy ðw` 


Exemplu. Să se afle derivatele parțiale ale funcţiei f(x,y) = xy? + x7y , cunoscând că x = 
utv,y=u-V. 
Soluție. Notăm (u,v) = f(x (u,v), y (u, v)) 


op ôf ôx əf ôy 5 > 3 3 
= —— x = (y +2xy) 1+ (2xy +x) l=y +4xy +y", 
a ooa a (y^ +2xy)-1+(2xy +x"): 1=y" +4xy +y 


Op of Ox of Oy 2 2 2 2 
— = —-— +: = (yl +2xy)-1+(2xy + )-l=y x. 
A a (y +2xy)-1+(2xy +x*^) 1=y" -x 
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Observaţie. Formula (9.4.3.1) se poate generaliza când f are n variabile iar fiecare este 
funcţie de m variabile. Evident aplicarea formulei este posibilă când toate funcţiile implicate sunt 
derivabile parțial. 

9.4.4. DIFERENŢIALELE FUNCŢIILOR 
DE MAI MULTE VARIABILE 


Fie f:DCR?—R, (x0,y0) şi (x,y) două puncte din D. Să notăm h = x - xo, 


k = y - yo. Numim creşterea totală a funcţiei în punctul (x0,yo) cantitatea Af (x9, Yo) definită prin: 


Af (xo: Yo) = f(x,y) — fo» Y0) (9.4.3) 


Definiţia 9.4.2 Funcţia f este diferenţiabilă în (x0,yo) dacă există două constante A şi B 


independente de h şi k şi două funcţii A(X, Yo, h,K), B(xo- Y0; h, K) astfel încât: 


O(X9> Yo „0,0) == 0 > B0, Y0-0,0) = 0 E) lim O(X9, Yo- h,k) = lim Bo, Yo- h,k) => 0 
h—0 h—0 


k—0 k—0 
Af (Xo; Yo) = Ah + Bh + ha(x o, Y0: h,K) + k- B(x0; Y0; h, K) (9.4.4) 
Se poate demonstra următoarea teoremă: 
Teorema 9.4.3 Dacă f (x,y) este diferențiabilă în (Xo,yo) atunci ea este derivabilă 


parțial în raport cu x şi y în acest punct şi Cru Yo) =A, bo Yo) =B. 
x 


Definiţia 9.4.3 Expresia Ah +Bk se numeşte diferenţiala funcției f în punctul (Xoya) şi 
se notează: 
df (x0,yo0) = Ah + Bk. (9.4.5) 
Dacă funcţia f (x,y) este diferențiabilă în (x0,yo) atunci diferenţiala sa df(xo,yo) e dată de 


of of 
df (x0, Y0) = ao Yo)h + ay (%0 Yok (9.4.6) 


Cum dx = h, dy = k obţinem diferenţiala de ordinul întâi a funcţiei f (x,y) în forma 


finală: 
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of of 
df (x9, Yo) = ao Yo )dx + & (So: Yo)dy 


Exemplu. Să se determine diferenţiala de ordinul întâi a funcției f(x,y) = xy? +2x°y în 
punctul (1,1). 
Soluție: f! (x,y) =y +6x°y fi(LD)=7 
' 3 ' 
f (x,y) =2xy + 2x (1) =4 
Obţinem: df(1,1) = 7dx +4dy. 


Observaţie. Diferenţiala de ordinul întâi ale funcției f (x,y) serveşte la aproximarea 


creşterii totale a funcției f, Af (x9, Yọ) conform formulei (9.4.3). 


Dacă f (x,y) este diferenţiabilă în toate punctele (x9, Yọ) unui domeniu D; CD, putem 


defini o corespondenţă: 
(Xo Yo) > df(xo- Yo) (V) (x0- yo) ED. 
Care se numeşte diferenţiala lui f şi se scrie df(x,y). 
Ea poate fi la rândul său diferenţiabilă. Diferenţiala sa se numeşte diferenţiala de ordinul 
al doilea a lui f şi se notează d”f(x, y). Un calcul elementar conduce la formula d7f(x,y). 
í Of ai 


Of O 
d’f(xo, Yo) = z Eo Yo)h? ag 2 (Xo: Yo)hk T o o Yo) k? 
Ox OxOy Oy 


Exemplu. Să se determine d"f(,1) pentru funcţia: 


f(x,y) =xy + 2x0y 


Soluţie. 
Of 3 Of 2.32) 
=y+4x — =x +6x 
E y y öy y 
2 2 2 
a Ya t =1+12xy , ÊT L 12x?y 
ox OxOy Oy 
2 2 2 
ara : Ei (,1)=13 , ien şi în final 
Ox OxOy Oy 
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d2F(,D) = 4h? + 26hk +12k? . 


Observaţii: 

1°. Procesul poate continua obţinând diferenţiale de ordin superior lui 2. 

20. Diferenţiala de ordinul al doilea este după cum se poate vedea din exemplul de mai 
sus o formă pătratică de variabile h,k. Când variabilele sunt x,,xp,...,X, această diferențială 
devine o formă pătratică de n variabile. Ea joacă rolul derivatei a doua la funcţii de o singură 


variabilă adică serveşte studiul convexităţii unei funcții de n variabile. 


9.5. EXTREMELE FUNCȚIILOR DE 
MAI MULTE VARIABILE 


9.5.1. FORMULA LUI TAYLOR PENTRU FUNCŢII 
DE MAI MULTE VARIABILE 


Se cunoaşte formula lui Taylor pentru funcţii de o variabilă care are derivate până la 


ordinul k+1 într-o vecinătate a punctului xo. 


2 
EO = fo) EP, EDP, 


(x —Xo)k i E (kH) 


0 fiind un punct intermediar lui Xo şi xX. 

Această formulă permite aproximarea valorii f în x cu ajutorul valorilor acesteia şi a 
primelor sale k derivate în punctul xo, dacă ultimul termen (restul de ordin k) este neglijat. 

Căutăm o exitindere a formulei la cazul funcțiilor de două variabile după care, natural, 
extinderea se poate face la n variabile. 

Fie aşadar, f : D C R° —> R şi (xo,yo) un punct interior din D. 


Presupunem că (x,y) este un punct oarecare din D şi notăm h =x —x9, k =y-—yo. 


127 


Teorema 9.5.1. Dacă f (x,y) are toate derivatele parţiale până la ordinul m în (x0.yo) 
iar derivatele parţiale de ordin m+1 există într-o vecinătate V a lui (x0,yo) atunci oricare ar 


fi punctul (x,y) € V , are loc formula: 


fy) ii a aa i 


(2) 
A (xo yoh Crono) e (9.5.1) 
ə (m+) 
riley ayk FA+Ra 
unde : 
Gi) 
Z aoyo Oe o le f reprezintă puterea formală “i” a expresiei diferenţiale 


din paranteză aplicată funcției f iar R, numit restul de ordin m în formula lui Taylor este 
exprimat prin: 


(m+) 


R = S o +h, yo +OK) + (they +OK f; OEOD (9.52) 
X 


9.5.2. TEOREMA LUI FERMAT PENTRU FUNCŢII 
DE MAI MULTE VARIABILE 


In teoria diferențială a funcțiilor de o variabilă, teorema lui Fermat afirmă că într-un 


punct de extrem local din interiorul domeniului de derivabilitate al funcției, D; , derivata f ! se 


anulează. Încercăm aici o extindere la funcții de două variabile care se poate apoi transpune la 


funcţii de n variabile. 


Să considerăm f:DCR?—R şi (Xo Yo) ED. 


Definiţia 9.5.1 Punctul (x), Yọ) se numeşte maxim (minim) local pentru f(x,y) dacă 


există o vecinătate V a acestuia astfel încât (V) (x,y) e V are loc: 


f(x,y) < fo Y0) € y) 2 f(xoY0)) (9.5.3) 
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Teorema lui Fermat poate constitui o metodă prin care se determină punctele de extrem 


sau posibilele puncte de extrem. 


Teorema 9.5.2 Dacă f(x,y) are derivate parţiale finite în raport cu x şi y în punctul 


de extrem local (x), yo) din interiorul lui D, atunci: 


f f 
Z koyo) =: Soyo) =0 (9.5.4) 


Demonstrație. Considerăm (x)=f(x, yọ) şi Yy) =f(39,y). Cum (Xo; Yo) este în 
interiorul lui D, (vezi fig. 5.3) xo se află în interiorul lui (a,b) iar yo în interiorul lui (c,d). Se poate 
aplica teorema lui Fermat pentru (x) cu xo punct de extrem local şi ẹ(y) cu yo punct de 
extrem local. 


Obţinem: 


f 
1) e 


dp of 
— = > = 0, — 
dis Xo) J (Xo-Y0) dy J 


ceea ce demonstrează teorema. 


Fig. 5.3 
Soluțiile sistemului (9.5.4) se numesc puncte critice sau puncte staționare. Reciproca 


nefiind adevărată, punctele critice rămân posibile puncte de extrem. Condițiile (9.5.4) pentru 


punctele (X9, Yo) sunt doar condiții necesare de extrem. 
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Exemplu. Să se determine posibilele puncte de extrem ale funcției: 
f(x,y) =x? +y —3xy +15 


Soluție: Calculăm 2A „Avem: 


x Oy 


A = 3x2 —3y, > =2y—3x 


Rezolvăm sistemul oE =0, —=0 
Ox 


3x? —3y =0 
2y—3x =0 


Obţinem soluţiile (0,0) şi E i a7) care sunt posibile puncte de extrem. 


9.5.3. CONDIȚII SUFICIENTE DE EXTREM 
(DETERMINAREA EXTREMELOR) 


Am văzut în cele de mai sus că un punct de extrem este punct staționar dar din 
teorema lui Fermat nu putem obține şi reciproca. Căutăm aşadar condiții suficiente de 
extrem prin care să putem aprecia dacă un punct staționar este sau nu este punct de 
extrem. 


Să remarcăm mai întâi că formula lui Taylor (9.5.1.2) se poate scrie şi sub forma: 
1 1 lra 
f(x,y) = f(xo> Yo) a, Yo) Eri EES f(x0 Yo) Rm (9.5.5) 


Dacă m = 2 atunci 


f(x,y) =f (x0: Y0) 4 Laf | et FR, (9.5.6) 


Să presupunem acum că (Xo,Yọ) este un extrem local. În virtutea condiţiilor (9.5.2.2) 


rezultă că : 


Of of 
df (x0, Yo) = o Yo)h ta Yok =0. 
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În plus dacă x —Xp=h şi y—y =k sunt suficient de mici, Rz poate fi neglijat (acolo h 
şi k apar la puterea a treia). Rezultă că pentru (x,y) într-o vecinătate suficient de mică a lui 
(X0Y0): 

sgn f(x,y) —f(x0.Yy0) =sen d’fF(xo, Yo) 


obținem, astfel: 


Teorema 9.5.3. Dacă forma pătratică d"f(x9,yg) este pozitiv definită în punctul critic 


(X0>Yg) atunci acest punct este minim local. Dacă în acelaşi punct (sf (X9>Y9) este negativ 
definită, punctul este maxim local. 
Condiţiile în care o formă pătratică este pozitiv (negativ) definită au fost prezentate 


anterior. Pentru prezenta formă pătratică, matricea este: 


Of of 
— (Xo: Yo) —— (Xo: Yo) 
Ox? OxOy 
H=| A (9.5.7) 
O Of 
azoy eW g e 


care se numeşte matrice hessiana şi se poate generaliza la funcții de n variabile. Suntem 


conduşi astfel la următoarea teoremă: 


Teorema 9.5.4 Fie P (xo, yọ)punct critic pentru funcția f(x,y) şi A, ^, minorii 


diagonali ai matricii Hessiene în acesta. Atunci: 


1. Dacă A >0 , A, >0 atunci P este punct de minim local pentru f(x,y); 
2. Dacă A, <0 , A, >0 atunci P este punct de maxim local pentru f(x,y); 
3. Dacă A >0, A, <0 sau A <0, A, <0 atunci P nu este punct de extrem local 


pentru f(x,y); 


4. Dacă A =0 sau A, =0 atunci nu putem preciza natura punctului critic P. 


Exemple : 


1°. Să se determine extremele funcției: 


f(x,y) =x? +y —3xy +15 
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Soluţie. Determinăm punctele de extrem local printre punctele critice adică printre 


soluţiile sistemului: 


f f ; YI 
Să = 0, O. = 0. Am văzut anterior că sunt două puncte critice (0,0), 22). 


ox Oy 


Fiecare poate să fie sau nu punct de extrem. 


Oi suc ai Pf 


ii E Ee i 2 
0x’ OxOy o 


Matricea hessiană: 


6x —3 
H(x, y) = e 4 


O calculăm în fiecare punct critic: 
0 —3| . 
H(0,0) = E i şi A, =0, A, =—9 
Punctul (0,0) nu este punct de extrem. 


2 Pe Mae E 
(2.2) 2 şi A=, A, =18 


=3 2 


Punctul E , a este un punct de minim. 


2°. Să se determine extremele funcţiei: 
SAAN ENEA. 
fF(x1;X2,X3) =X] +X3 +X3 —X1X3 
Punctele staționare (critice) sunt soluții ale sistemului 


CETE =0, DE cag =0, A E =0 
9x, x, Ôx; 


Obţinem un singur punct critic (0,0,0). Acesta este şi singurul posibil punct de extrem. 


E ME N E A 
ôx? ” @? ° a? ” ôxðy, ” x ðx 
AORN e Ana ate AE ue 

ðx;ðx, Oxi, OO xð, 
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Matricea hessiana devine: 


2 0 -l1 
H(xX2X3)=| 0 2 0 
-1 0 2 


şi are lanțul minorilor principali: 
A2 A= Ag =6 


Rezultă că (0,0,0) este un punct de minim local. 
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